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Introdução 
às Equações 
Diferenciais de 
Primeira Ordem 


Neste primeiro capítulo vamos discutir os principais conceitos associados às 
equações diferenciais, suas origens e evolução, analisar alguns exemplos prá- 
ticos e apresentar a montagem das equações diferenciais e seus diferentes 
tipos e os meios para identificá-los. 

H7 OBJETIVOS 


■ Refletir sobre os conceitos básicos de equações diferenciais. 

■ Entender o processo de construção de uma equação diferencial. 

■ Associar as equações diferenciais aos conteúdos de derivadas e integrais. 

■ Saber classificar as equações diferenciais. 
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1.1 Conceitos Básicos 


Historicamente, as equações diferenciais começaram a ser estudadas por Isaac 
Newton e Gottfried Wilhelm Leibniz no século XVII. Os irmãos Bernoulli tam- 
bém estudaram equações diferenciais. Além deles, Leonhard Euler, Joseph- 
Louis Lagrange e Pierre-Simon de Laplace foram importantes na criação de 
métodos para resolução de equações diferenciais assim como no estudo de suas 
aplicações. Outros matemáticos como Bessel, Legendre, Hermite, Chebyshev e 
Hankel também colaboraram, como colocado em BOYCE e DiPRIMA (2014) 
para o desenvolvimento da área, levando seus nomes em diversos métodos e 
equações em homenagem a eles. Além deles, Jean-Baptiste Hoseph Fourier, 
que também foi aluno de Laplace e Lagrange fez estudos importantíssimos na 
área de equações diferenciais que possuem aplicações tecnológicas importan- 
tes nos dias de hoje nas áreas de telecomunicações e em processamentos de 
imagens digitais. 

Ainda segundo BOYCE e DiPRIMA (2014), a evolução nos estudos das equa- 
ções diferenciais se confunde com a evolução do próprio Cálculo Diferencial 
e Integral, já que as duas áreas guardam uma relação muito próxima entre si. 
Em especial, a maior impulsão em ambas está na articulação das teorias com 
as aplicações práticas, ou seja, na interdisciplinaridade da Matemática com ou- 
tras áreas do conhecimento, como a Física, Química, Biologia, Economia, entre 
outras. 

Pensando por este lado, é natural associar os estudos das equações diferen- 
ciais a Matemática Aplicada, devido à amplitude de problemas práticos que po- 
dem ser caracterizados através de equações diferenciais. 

Mas, do que se trata essa caracterização matemática? 

O termo caracterização matemática está relacionado com a arte da constru- 
ção de modelos matemáticos, ou simplesmente modelagem. O objetivo da mo- 
delagem é representar um problema qualquer por meio de equações, no caso, 
equações diferenciais. Na construção dos modelos, como bem explicado em 
BOYCE e DiPRIMA (2014), é importante que sejam identificadas, com correção, 
as variáveis envolvidas e determinar corretamente qual a dependente e qual a in- 
dependente, bem como, atribuir unidades de medidas adequadas e selecionar o 
princípio ou lei matemática que representa a relação funcional entre as variáveis. 

Nesta disciplina, abordaremos tais modelos matemáticos em situações 
onde o princípio matemático é dado por uma equação diferencial. 
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Sim, mas uma questão essencial neste início de conversa é: 


0 que é uma equaçao diferencial? 


Para responder a essa pergunta, convidamos você a refletir sobre outra 
questão: do que se trata o conceito de diferencial no contexto do Cálculo? 

Se você pensou em derivadas, está com a razão. O Cálculo Diferencial está 
associado com derivadas e as taxas de variações de funções matemáticas. 

Você estudou diversas situações práticas onde a taxa de variação era um fe- 
nômeno de interesse, tão importante quanto avaliar a própria função. Através 
das técnicas de derivadas, era possível avaliar a taxa de variação e verificar o que 
ocorria com a função para um dado ponto do domínio. 

Neste sentido, BOYCE e DiPRIMA (2014) se refere às equações diferenciais 
como as relações matemáticas que utilizam a taxa de variação de uma função 
no equacionamento da relação, ou seja, uma equação que além das variáveis e 
da forma funcional, também possui uma derivada como componente. Assim, 
uma equação diferencial nada mais é do que uma equação que envolve variá- 
veis e derivadas dessas variáveis. 

Neste contexto, você pode argumentar: quais situações práticas podem ne- 
cessitar analisar suas taxas de variação e derivadas, juntamente com as variáveis? 

Bem, apenas a título de apresentação, podemos considerar questões rela- 
cionadas com o movimento de fluidos, a variação da corrente elétrica em circui- 
tos, o comportamento da temperatura na superfície de objetos, a propagação 
de fenômenos sísmicos e também o comportamento relacionado ao cresci- 
mento populacional. 

Continuando a analogia com o estudo das funções e equações, dependente 
do tipo de equação ou função, você se via obrigado a utilizar um rol de técnicas 
apropriadas para a manipulação da mesma. Por exemplo, em se tratando de 
uma equação exponencial, utilizaria propriedades de exponenciação e de loga- 
ritmos, caso fosse uma equação trigonométrica, as relações do ciclo trigonomé- 
trico e outras mais. 

Da mesma forma, podemos ter diferentes tipos de equações diferenciais, 
mas como sabemos que essas equações possuem derivadas como componente, 
é natural supor que a construção, representação, resolução e interpretação das 
equações diferenciais exigirão domínio sobre as técnicas de derivação, e por 
consequência, domínio sobre as técnicas de integração. 
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As derivadas ou taxas de variação representadas nas equações diferenciais 
podem ser parciais ou não, de primeira ordem ou de ordem superior. Nesta dis- 
ciplina trataremos de todos os tipos e classificações principais, focando prin- 
cipalmente as aplicações dentro da área de Engenharia que vai desde a área de 
Mecânica e Ondas até Termodinâmica, Eletromagnetismo e Eletrônica. 

A área de equações diferenciais como modelagem vem assumindo um pa- 
pel bastante importante nos últimos tempos, sendo uma das áreas mais pro- 
missoras da matemática aplicada. 

Agora que você está mais familiarizado com o termo equação diferencial, 
deve estar se questionando o seguinte: 

Como é possível resolver uma equação que envolve uma derivada? 

Bem, para responder a essa pergunta, faremos uso de seus conhecimento 
de Cálculo Diferencial e Integral, para mostrar como, de uma forma genérica, 
podemos resolver uma equação diferencial. 

1 .2 Soluções de equações diferenciais. 

Existem muitos caminhos para resolver uma equação diferencial. Para enten- 
der isso, basta que você retroceda nos seus estudos básicos de funções e lem- 
brará que para cada tipo de função existe um protocolo para resolução. 

O mesmo ocorre com as equações diferenciais e tais protocolos, modelos de 
resolução serão detalhados oportunamente neste livro. 

Por hora, vamos apenas ilustrar de forma geral o que ocorre quando parti- 
mos para a solução de uma equação diferencial. 

Para iniciar nosso debate sobre a resolução, temos que, obviamente, dispor 
de uma equação diferencial para trabalhar. Já discutimos que uma equação di- 
ferencial, em sua forma mais simplificada, agrupa uma variável dependente y, 
uma variável independente x e a derivada de y em relação a x que pode ser repre- 
sentada por — , pory' ou pory'. 
dx 

No exemplo 1 temos a representação de um tipo de equação diferencial. 

Todas as equações diferenciais resolvidas neste capítulo utilizam o método da solução 
por integração direta, conforme exposto por EDWARDS e PENNEY (1995). 
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C7 EXEMPLO 

Resolver a equação diferencial — = 7 x 2 + 2x 

dx 

De uma forma geral, o objetivo na resolução da equação diferencial é reduzi-la para a 
forma y = f(x), ou seja, encontrar a função primitiva (original) cuja derivada é exatamente a 
equação diferencial, como pontuado em ZILL e CULLEN (2001). 

Para isolar a variável dependente y na expressão, temos que inicialmente isolar a dife- 
rencial de y, da forma: 

d ^ = 7x 2 +2x 
dx 

dy = (7x 2 +2x)dx 


Como o termo à esquerda é uma diferencial, para resolver a questão temos que utilizar 
uma integral indefinida, de ambos os lados da equação, como explicitado em EDWARDS e 
PENNEY (1 995). Assim, temos: 

dy = (7x 2 +2x)dx 
jdy = \(7x 2 +2x)dx 


Resolvendo a integral em ambos os lados, temos: 

jdy = j(7x 2 +2x)dx 

jdy = j(7x 2 +2x)dx 

y + C, =— x 3 +x 2 +C„ 

1 3 2 

y = lx 3 + x 2 + (C 2 -C 1 ) 

ó 

y = l X 3 +x ‘2 +C 

3 

Assim, a equação diferencial está resolvida e essa solução, como apontado em BOYCE 
e DiPRIMA (2014) é chamada de Solução Geral. 
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Duas considerações sobre a solução de uma equação diferencial devem ser feitas neste 
momento. 

Em primeiro lugar, para verificar a validade dessa solução, se está ou não correta, você 

pode simplesmente derivar ambos os lados da expressão e verificar se reproduz a função 

d ^ = 7x 2 +2x. 
dx 


Considerando o resultado anterior, temos: 


= -x 3 + x 2 +C 
3 


df-x 3 +x 2 +C 

± = A3 

dx dx 


dy 

dx 

dy 

dx 

dy 

dx 


d-x 3 

3 

dx 



dC 

dx 


7 

3 


3x 2 + 2x + 0 


= 7x 2 + 2x 


Veja que o resultado é exatamente a equação diferencial original, o que confirma que a 
resposta está correta. 

A segunda observação a ser feita é que a solução é geral, ou seja, dependendo do valor 
da constante C, pode assumir infinitas soluções. 

Uma forma de ilustrar a questão é através da figurai .1, considerando 3 valores diferen- 
tes para C. 



Figura 1.1 - Possíveis soluções para o exemplo 1 , com C = 1 , C = 4 e C = 10. 


CAPÍTULO 1 ■ 1 I 


Note que, para possibilitar a construção dos gráficos ilustrados na figura 1.1, fixamos 
valores escolhidos de forma aleatória para a constante C, gerando assim uma solução 
particular da equação diferencial para cada C. Contudo, essa não é única forma de obter 
valores particulares para a equação diferencial. 

Uma forma de obter uma Solução Particular para a equação diferencial, de acordo com 
EDWARDS e PENNEY (1 995) é associar um resultado numérico para a mesma. 

No caso do exemplo 1 , podemos estar interessados em uma solução que torne y= 6 para 
quando x = 1, ou sejay(1) = 6. Nesses casos, o valor de C não é escolhido aleatoriamente, 
mas definido com base na condição. Assim, temos: 


Assim, uma solução particular para a equaçao diferencial do exemplol, considerando a 
condição imposta, é: 



18-3-7 

3 


3 


y = —x 3 +x 2 + 


7 s . o . 8 


3 3 


Vejamos outro exemplo. 
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XJ EXEMPLO 

Resolver a equação diferencial < ^- = 1 0 - — 

dx 3 

Procedendo da mesma forma do exemplo 1 , isolamos dy na equação, de forma que: 

^ = 10-' 
dx 3 

dy _ 30- y 

dx 3 

3dy = {30- y)dx 

dy dx 

(30-y) ~ ~3~ 

A expressão anterior pode ser escrita convenientemente da forma: 
dy dx 

(y-30) _ 3~ 

Integrando em ambos os lados, temos: 
r 1 dy c 1 dx 
J (y-30) ” ~3 

ln|y-30| + C 1= -| + C 2 
ln|y-30|=-| + (C 2 -C 1 ) 
ln|y-30| = -| + C 

Colocando ambos os lados na base exponencial, temos: 

e /n|y-30| _ e (“3 +C ] 

g/n|y-30| _ g 3 gC 

Como e também é uma constante, podemos escrever e c = C. 

X 

y-30 = Ce 3 

E fínalmente chegamos a solução geral da equação diferencial que é: 

X 

y — Og 3 + 30 
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Variando o valor da constante C, podemos encontrar diferentes soluçoes para a equaçao 
do exemplo 2, como ilustrado na figura 1.2. 



c=i — c=2 — •— C=4 


Figura 1 .2 - Possíveis soluçoes para o exemplo 2, com C= 1 , C=2 e C=4. 


[7 EXEMPLO 


Resolver a equação diferencial 4x - y'= 1 

Isolando a derivada e escrevendo y’em termos de dy e dx, temos: 


= 1-4* 
dx 


dx 


= 4 * -1 


Calculando a integral, temos: 


\dy=j(4x-f)dx 
n 4x2 n 

y +c i - -x+c 2 

y = 2x 2 -x + (C 2 ~C ] ) 
y = 2x 2 -x + C 

Esse resultado é a solução geral da equação diferencial, com C uma constante qualquer. 
É fácil verificar que a solução está correta, basta derivar e verificar que o resultado é igual 
a equação diferencial original. 

Da mesma maneira já feita nos exemplos 1 e 2, variando o valor de C temos soluções 
particulares para a equação diferencial, como exibido na figura 1 .3. 
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Figura 1 .3 - Possíveis soluções para o exemplo 3, com C=0, C=3 e C=5 

Ainda, é possível determinar uma solução particular para a equação diferencial colocan- 
do uma condição. Por exemplo, colocando que o valor de y(2) = 5. Colocando essa condição 
na solução geral, temos: 
y = 2x 2 - x + C 
y(2) = 2.2 2 - 2 + C = 5 
8-2+C=5 
6 + C = 5=>C = -1 

Assim, uma solução particular da equação diferencial, a partir da condição inicial pro- 
posta é: 

y = 2x 2 -x-1 


[y EXEMPLO 

_ _ . , dy x 

Resolver a equaçao diferencial — = — 

dx y 


Reescrevendo a equação diferencial e calculando a integral, temos: 
Jy dy = | - x dx 


K 2 x 2 

— + G = + C 9 

2 i 2 2 

I / 2 X 2 

Z_ + _ = (c 2 _ Cl ) 

^+^ = c 
2 2 

y 2 + x 2 = 2C 
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Agora, como C é uma constante, podemos considerar, sem incorrer em nenhum erro 
nas condições das equações diferenciais, que a constante 2C pode ser substituída por c 2 , 
obtendo a seguinte solução geral: 


y 2 + x 2 = C 2 

Note que com a mudança na constante, a solução geral corresponde a equação de uma 
circunferência de C. Assim, variando o valor de C, temos diferentes circunferências que re- 
presentam soluções para a equação diferencial. 

Como a solução geral de uma equação diferencial possuirá uma constante de integra- 
ção desconhecida, muitas vezes, para melhor representação da equação será necessário 
escrever a constante de maneiras especiais. Veja o exemplo 4, escrevemos a constante C 
como sendo C2. Isso não causa nenhum problema, pois a verificação da solução é feita 
com base da derivada da solução geral. Assim, se C é uma constante, C2 também é uma 
constante, ou eC também, e sendo assim a derivada de qualquer uma delas será igual a 
zero. Desta forma, não importa como a constante C será representada. 

A partir do capítulo 2, veremos diferentes métodos para a resolução de uma equação di- 
ferencial, baseados em certas propriedades presentes nas equações diferenciais. Mas, antes 
disso, vamos verificar inicialmente quais as características de uma equação diferencial que 
lhe conferem diferentes classificações. 


1.3 Classificações das equações diferenciais. 

Você deve se lembrar de que uma função matemática pode ter diferentes classi- 
ficações dependendo de alguns parâmetros que as identifiquem. Por exemplo, 
as funções polinomiais são classificadas em função do grau atribuído a elas, e 
dependendo do grau possuem características diferenciadas. Se o grau for igual 
a um, temos uma função linear que sempre apresenta um zero nos reais, de 
grau 2, uma função quadrática, que pode ou não apresentar dois zeros nos re- 
ais, e assim por diante. 

Com as equações diferenciais ocorre o mesmo. Dependendo do parâmetro 
escolhido para caracterizar a equação, diferentes classificações são obtidas e 
essas classificações podem influenciar decisivamente a existência de solução 
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e, caso exista, na forma de resolução das mesmas. 

Temos, de uma forma geral, cinco diferentes formas de classificar uma 
equação diferencial, a saber: 

• Classificação quanto ao número de variáveis e da derivada presente na 
equação; 

• Classificação quanto ao número de equações diferenciais; 

• Classificação quanto a ordem das derivadas; 

• Classificação quanto a existência ou não de solução; 

• Classificação quanto a linearidade da equação. 

Vejamos detalhadamente cada uma das classificações destacando a forma 
geral em cada uma das possibilidades. 

1 .3.1 Classificação quanto ao número de variáveis e o tipo da 
derivada. 

Fazendo novamente uma analogia aos estudos realizados anteriormente, você 
deve se lembrar de que quando uma função matemática apresenta apenas de 
uma variável dependente (y) e uma variável independente (x), y=f(x), sua deri- 
vada é dita total e depende apenas de x, da forma: 

dy 

y ' = Tx 

Contudo, também, existem funções matemáticas que relacionam uma vari- 
ável dependente (z) com duas variáveis independentes ( xey ), z = f(x,y), ou mais, 
e com isso temos as chamadas derivadas parciais que podem depender de uma 
ou outra variável independente, da forma: 

õz õz 

z -= — z-= 

õx õy 

Desta forma, como as equações diferenciais dependem das derivadas, se a 
equação envolver apenas uma variável independente e derivadas totais, con- 
forme colocado por ZILL e CULLEN (2001), temos uma Equação Diferencial 
Ordinária , que chamaremos muitas vezes de EDO. Nos exemplos 4, 5 e 6 ilus- 
tramos algumas equações diferenciais ordinárias: 
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C7 EXEMPLO 

Seja a seguinte equação diferencial: 


Esta equação diferencial é dita ser ordinária pelo fato de conter apenas derivadas totais e 
depender apenas de x. 


E y EXEMPLO 

Seja a seguinte equação diferencial: 

dy_ 

dx 


Este exemplo também apresenta uma equaçao diferencial ordinária, pois não temos a 
presença de derivadas parciais. 


E7 EXEMPLO 

Seja a seguinte equação diferencial 


y-+3y = e 5t 

Como neste caso também temos apenas derivadas totais, trata-se de uma EDO. 

Agora, se a equação diferencial envolver mais de uma variável independente e derivadas 
parciais, temos uma Equação Diferencial Parcial, que podem ser representadas por EDP. Nos 
exemplos 7 e 8 temos algumas equações diferenciais parciais: 
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XJ EXEMPLO 

Seja a seguinte equação diferencial: 

õz õz 
— + 4- — 
õx õy 


Esta equaçao diferencial é dita ser parcial, pois depende de duas variáveis independen- 
tes e contém derivadas parciais. 


XJ EXEMPLO 

Seja a seguinte equação diferencial: 

õ 2 y _ õ 2 y 
õt 2 õx 2 

É uma equação diferencial parcial pelo fato de contar derivadas parciais e depender de 
mais de uma variável independente, xe t. 

O objetivo da disciplina é estudar todas as propriedades e métodos de resolução das 
equações diferenciais ordinárias e não serão abordadas as equações diferenciais parciais 
na disciplina. 


J CONEXÃO 


As equações diferenciais parciais ocorrem quando a equação diferencial contém mais de 
uma variável independente e derivadas parciais dessas variáveis. Recomendamos que leia as 
referências abaixo para ampliar seus conhecimentos sobre as equações diferenciais parciais: 
http://www.mat.ufmg.br/~lima/apostilas/EDP1 -aulas-22-8-1 3.pdf 
http://www.im.ufrj.br/~medeiros/LinkedDocuments/livrometclassicos.pdf 
Você pode também assistir aos vídeos produzidos pela IMPA sobre esse assunto: 
http://video.impa.br/index.php?page=doutorado-201 3-equacoes-diferenciais-parciais- 
-e-aplicacoes 
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1.3.2 Classificação quanto ao número de equações diferenciais. 

Da mesma forma como ocorre com as equações, existem cenários onde são 
necessárias mais de uma equação ou mais de uma função para caracterizar o 
problema, o que provoca uma mudança de uma equação para um sistema de 
equações. 

Com as equações diferenciais ocorre da mesma forma. Podemos trabalhar 
com apenas uma equação diferencial, como no caso dos exemplos de 1 a 5, ou 
então, devido à necessidade de trabalhar com duas ou mais funções, como de- 
finido por BOYCE e DiPRIMA (2014), termos um sistema de equações diferen- 
ciais, como o exibido no exemplo 9. 


ty EXEMPLO 

Em uma situação prática qualquer foram definidas duas equações diferencias, gerando um 
sistema de equações diferenciais, da forma: 

^- = 5x + 4 xy 


Na disciplina de Cálculo III limitaremos nossa abordagem aos casos de uma única equa- 
ção diferencial, mas você pode consultar a referência recomendada para saber mais sobre 
os sistemas de equações diferenciais. 


Ey CONEXÃO 

Os sistemas de equações diferenciais ocorrem quando temos duas ou mais equações di- 
ferenciais sendo trabalhadas de forma simulfânea. Recomendamos que leia as referências 
abaixo para ampliar seus conhecimentos sobre os sistemas de equações diferenciais: 
http://www.dm.ufscar.br/profs/waldeck/sourceforge/pngtest.php 
http://www.ime.uerj.br/~calculo/LivrolV/sistemas.pdf 
http://www.mtm.ufsc.br/~daniel/sem1_05/edo/farlow/sec7.pdf 
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1 .3.3 Classificação quanto a ordem das derivadas. 

Essa classificação é feita de forma similar a classificação dos polinómios, ou 
seja, o grau mais elevado de uma variável em um polinómio caracteriza a or- 
dem do mesmo. Por exemplo, o polinómio P(x)=4x 5 +7x í -2x 2 é de ordem 5, ou de 
quinto grau, já que o maior grau da variável xé 5. 

Para uma equação diferencial, basta examinarmos qual a derivada de maior 
ordem envolvida na equação, e teremos o grau da equação diferencial, como 
apresentado em EDWARDS e PENNEY (1995). 

Vejamos alguns exemplos. 


XJ EXEMPLO 


Qual a ordem da equaçao diferencial abaixo? 


dt 3 dfi 


Como a derivada de maior grau é 
ceira ordem. 


d 3 y 

dt 3 


a equação diferencial é de ordem 3, ou de fer- 


XJ EXEMPLO 


Qual a ordem da equaçao diferencial abaixo? 

y' y + 2y - 6xy = 10 

Como a derivada de maior grau é y\ a equação diferencial é de ordem 1 , ou de primeira 
ordem. 


EXEMPLO 


Qual a ordem da equaçao diferencial abaixo? 


dy at 

■ — + ay ■ e at 
dt 


e at g(t) 


CAPÍTULO 1 - 21 





Como a derivada de maior grau é Oi , a equaçao diferencial é de ordem 1 , ou de primeira 

dt 

ordem. 

Resumidamente, na tabela 1 temos a seguinte as seguintes relações de ordem das 
equações diferenciais e o respectivo grau, que está associado ao expoente da derivada de 
maior grau tomando a equação diferencial na forma de um polinómio. 


EQUAÇÃO DIFERENCIAL 


GRAU 

y’ = 2x 3 y 

1 

i 

y" = 2y + sen(x) 

2 

i 

y”' = y(x+1) 

3 

i 

(y’) 2 + x.cos(y) = 0 

1 

2 

(y”) 2 + 5xy = 0 

2 

2 

(y”’) 2 + 3y’+2y =e x 

3 

2 

(y’) 3 + x.e x = 0 

1 

3 

(y”) 3 + 6xy' = 0 

2 

3 

(y”’) 3 + y"+2y' =e x 

3 

3 


Tabela 1 - Ordem e grau de equações diferenciais. 


1 .3.4 Classificação quanto a existência de solução. 

Você pode observar nos exemplos 1, 2 e 3 que trabalhar com equações diferen- 
ciais envolve a resolução de integrais, e no caso dos exemplos, de forma direta. 
Notou ainda que temos dois tipos de soluções, uma solução geral e uma solu- 
ção parcial, obtido mediante uma condição prévia. 

• Assim, como para as equações reais, para as equações diferenciais pode- 
mos ter duas situações, equações diferenciais com solução e equações diferen- 
ciais que não possuem soluções. 
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O grande problema em trabalhar as soluções de equações diferenciais é o 
grau de complexidade da mesma que pode exigir um esforço significativo na 
resolução. Contudo, existe um elemento importante neste contexto, caso você 
obtenha uma solução, não é difícil verificar se a mesma é correta ou não, como 
mostrado nos exemplos anteriores. 

É particularmente saudável que você se acostume a fazer a verificação para 
se certificar que a solução é de fato verdadeira. Lembre-se, determinada a so- 
lução, basta derivar o resultado e verificar se corresponde a equação original. 

Ainda, existe outra maneira de trabalhar com equações diferenciais que 
pode auxiliar nos casos mais simples. Não é uma técnica muito eficiente, mas 
é simples, como explicitado em BOYCE e DiPRIMA (2014). Consiste em suge- 
rir uma solução geral através da observação da própria equação diferencial. 
Vejamos dois exemplos. 


& EXEMPLO 

Seja a equação diferencial abaixo: 


d 2 y 

~dx 2 


-y = o 


Será que a função y = e x é uma solução para essa equação diferencial? 

Para avaliar a questão, basta substituir na equação diferencial e realizar os cálculos ne- 
cessários: 


0,2 ( e *) 

dx 2 

d ( e *) c » 
dx 

e x - e x = 0 


Assim, a função y=e*é uma solução, assim como a função y = e _x , também, é uma solu- 
ção. Desta forma, as funções y = Ce* e y = Ce~* também são soluções gerais, ou seja, temos 
uma família de soluções para a equação diferencial. 
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C7 EXEMPLO 

Seja a equação diferencial abaixo: 

Será que a função y = sen x é uma solução para essa equação diferencial? 
Substituindo, temos: 


d 2 y 


+ sen x = 0 


d 2 (sen x) 

d7 2 

d cos x 
dx 

-sen x + sen x = 0 => 0 = 0 


+ sen x = 0 


Assim, a função y = senx é uma solução, assim como a função y = -cosx , também, é 
uma solução. Desta forma, as funções y = C ■ (sen x) e y=C ■ (- cos x) são soluções gerais 
da equação diferencial. 


1 .3.5 Classificação quanto a linearidade da equação. 


A linearidade de uma equação diferencial, fator extremamente relevante para 
a seleção de técnicas de resolução, está relacionada à forma como as variáveis 
e as derivadas se associam na equação. Simplificadamente, sempre que /es- 
tiver multiplicando dy, ou / se apresentar não linear, a equação diferencial é 
não linear, mas se/e dpestiverem em forma somativa com/linear, a equação 
diferencial será linear. Vejamos um exemplo de cada tipo. 


C7 EXEMPLO 


Seja a equaçao diferencial abaixo: 

4— ^- + 5y = 10 
dm d 


É uma equação diferencial linear, pois não temos o produto da variável dependente y pela 
sua derivada e y aparece na forma linear com expoente igual a 1 . 
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XJ EXEMPLO 

Seja a equação diferencial abaixo: 
y.y" - y ' + sen(yt) = 1 5 

É uma equação diferencial não linear por dois motivos, a ocorrência do produfoy.y”e pelo 
fato de y estar associada a uma função não linear (seno). 

Bem, neste primeiro capítulo nos concentramos em fazer uma apresentação dos conceitos 
de uma equação diferencial, bem como dar os primeiros passos para a resolução das mesmas. 

Ainda, e igualmenfe importante, discutimos os elementos que utilizamos para classificar 
essas equações, especialmenfe a ordem e a linearidade das mesmas. 

Para finalizar a apresentação das equações diferenciais, é importante mencionar que o 
estudo das soluções de uma equação diferencial de primeira ordem pode ser feito através 
de campos de direção. Conforme colocado por BOYCE e DiPRIMA (2014) os campos de 
direção são construídos calculando-se o valor da função para uma grande variação das 
variáveis e dispondo os resultados em um plano cartesiano, criando assim uma malha que 
possibilita avaliar o comportamento da função e da equação diferencial. Como depende de 
muitos cálculos, a construção é feita computacionalmente. 


XJ EXEMPLO 

Supondo a equação diferencial y’= 6xy + 3x. O campo de direções associado a equação 
diferencial fica como ilustrado na figura 
y 

10 
5 
0 
-5 
-10 

-10 -5 0 5 10 

x (10.7424,-10.2523) 

Figura 1.4 - Saída do aplicativo Máxima com o Campo de direções da equação diferencial 
do exemplol 7. 
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Os campos de direção podem ser construídos a partir do aplicativo livre Máxima, que 
pode ser obtido a partir do seguinte endereço: 

http://br.ccm.net/download/baixaki-2 1 68-maxima 

Para a utilização do mesmo recomendamos a leitura dos seguintes materiais: 
www.ime.unicamp.br/~marcio/ss2006/grupo1 0.pdf 
www.sobralmatematica.org/monografias/milsangela.pdf 


No capítulo 2 iniciaremos as técnicas, as principais pelo menos, existentes para a reso- 
lução das equações diferenciais ordinárias de primeira ordem, EDO de primeira ordem, mas 
antes, vejamos alguns exemplos ilustrativos que originam as equações diferenciais. 


1.4 Situações práticas modeladas via 
equações diferenciais. 


Antes de iniciarmos a discussão sobre os métodos de resolução, vamos, a título 
de ilustração, discutir três situações práticas que podem ser representadas por 
meio de uma equação diferencial. 


EXEMPLO 

[Adaptado de EDWARDS JR e PENNEY (1995)]. Suponhamos que você acenda o forno 
do fogão que existe na cozinha de sua casa e coloque, em seu interior, uma assadeira de 
alumínio e a deixe por lá tempo suficiente para atingir uma temperatura elevada. Ao retirar 
a assadeira do forno e colocá-la sobre uma superfície qualquer dentro da cozinha, o que 
ocorrerá com a assadeira? Se você está pensando que a assadeira deixará de ficar quente 
e passará para um estado mais frio, está com a razão. Na verdade, a assadeira está em um 
processo de resfriamento, troca de temperatura com o ambiente e com a superfície. 

Se considerarmos que o ambiente possui uma temperatura 7, igual em qualquer ponto 
do ambiente, e que o objeto, a assadeira, também possui uma temperatura, T ob , constante 
ao longo de todo o objeto, e levando em consideração que 7 < T ob , temos, então, que com 
a variação do tempo f, existe uma perda de calor do objeto para o ambiente que ocorre em 
função do tipo do material do objeto. 
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Assim, podemos considerar que a taxa de variação da temperatura da assadeira, T ob , em 
função do tempo t é uma função da diferença entre a temperatura da assadeira e a tem- 
peratura do ambiente (T ob - F), e essa diferença, o decaimento da temperatura, a troca de 
calor, ocorre de forma ponderada pelas propriedades físicas do material. Chamando de a a 
constante que relacionada o comportamento do tipo de material da assadeira, alumínio, por 
exemplo, podemos representar a situação através de uma equação diferencial (considerando 
o sinal negativo pelo fato de ser uma perda de calor por parte da assadeira): 




EXEMPLO 

[Adaptado de ZILL e CULLEN (2001)] Estudos relacionados a planejamento urbano, 
impactos ambientais, entre outros dependem fundamentalmente da projeção de crescimento 
da população que ocupa a área geográfica em questão. Situações que envolvem cresci- 
mento populacional, não necessariamente humana, possuem uma característica peculiar, a 
presença de fatores que tendem a desacelerar o ritmo de crescimento, como predadores da 
população. Em relação a população humana, não temos predadores, mas temos epidemias, 
violência urbana e acidentes de larga escala que funcionam como tal. 

Dessa forma, suponhamos que em uma dada população, com uma quantidade inicial de 
P pessoas, a taxa de variação da população em função do tempo é proporcional a população 
atual, ou seja, quanto maior a população, a taxa de crescimento tende a ser maior também, 
de forma proporcional. 

Contudo, é importante considerar que a população também apresenta uma taxa de de- 
caimento provocada pelos fatores de desaceleração. Suponhamos que nessa população 
ocorra uma epidemia de dengue e que W indivíduos contraiam a doença e que (P-W) não. 
Assim, a taxa de crescimento da população doentes será proporcional a uma relação ponde- 
rada entre os doentes e os não doentes, da forma: 


dW 

dt 


k-W-(P-W) 


Onde k é uma constante de proporcionalidade. 


CAPÍTULO 1 ■ 27 



C7 EXEMPLO 

[Adaptado de BOYCE e DiPRIMA (2014)] Do alto de um prédio de altura H 0 , um objeto 
é lançado em direção ao solo. Assim, como podemos descrever a queda desse objeto em 
função do tempo? Bem, para representar o problema de uma forma simplificada, podemos 
supor inicialmente que o objeto possua uma massa m que permanece constante ao longo 
de toda a queda; em segundo lugar, podemos supor que a queda do objeto não seja afetada 
pela resistência do ar, mas apenas pela ação da gravidade g. Ainda, é importante lembrar que 
a ação da gravidade age acelerando a queda do objeto, mas a aceleração é uma derivada da 
velocidade que por sua vez é uma derivada do deslocamento. Assim, temos duas derivadas 
envolvidas na representação da queda do objeto, que pode ser representada por uma equa- 
ção diferencial da forma: 


d 2 H 



Neste momento, antes de prosseguir para o capítulo 2 e iniciar os estudos das técnicas 
de resolução de uma EDO de ordem 1 , convidamos você a resolver os exercícios abaixo que 
têm o intuito de reforçar alguns conceitos trabalhados neste capítulo. A resolução detalhada 
dos mesmos você encontrará no final dos exercícios propostos. 


E§T ATIVIDADES 

01. Resolva as equações diferenciais dadas abaixo, encontrando a sua solução geral, ve- 
rificando a veracidade da mesma e obtendo a solução particular em função da condição 
informada. 

a) — + 3 = senlx , com y(0) = 3. 
dx 

b) Y ~ Xe , com y(0) = 1 . 

c) cos^xy' = Se ” — , com y(0) = 2. 


28 ■ CAPÍTULO 1 




02. Classifique cada uma das equações diferenciais abaixo quanto ao seu grau 


x 2 y" + xy' + 7xy = Ln(x + 1) 

b) ^. (2 , + 1) + ^(e*)=1 

dv V ’ dv 3y ’ 

c) y' + 3y = 5 

03. Classifique as equações diferenciais em lineares ou não lineares. 

a) y” +( x -i)y" + ^y( x - = 4 

d 3 t 

b) — - + e 3 *=cosA' 
dt 3 

c ) xy” + xy + 3y = e x 

04. Resolva a equação diferencial 2y' + y = 0 e represente graficamente a solução geral 
para C igual a 5 e 10. 


ET REFLEXÃO 

Você deve ter notado que trabalhamos com equações diferenciais de primeira ordem e todos 
os exercícios foram resolvidos da mesma forma, utilizando integrais indefinidas. Contudo, 
podemos refletir sobre o seguinte: se existem diferentes classificações para as equações 
diferenciais e diferentes formas da equação ser representada, será que podemos utilizar 
essas particularidades para obter formas diferenciadas e talvez mais simples de resolver uma 
equação diferencial? 


^ CONCEITO 

Alguns termos citados neste capítulo são de especial importância e aparecerão nos capítulos 
futuros, sendo assim reforçamos seus conceitos neste glossário. 

1 . Equações diferenciais: são equações especiais que unem, na mesma sentença, uma 
função matemática, suas variáveis e suas derivadas; 
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2. Equações diferenciais de primeira ordem: são equações diferenciais que contém 
apenas derivadas de primeira ordem na sua formulação; 

3. Equações diferenciais de segunda ordem: são equações diferenciais onde a deri- 
vada de ordem mais alfa é a derivada de segunda ordem; 

4. Equações diferenciais ordinárias: são equações cuja função desconhecida depen- 
de apenas de uma variável. 

5. Equações diferenciais parciais: são aquelas equações cuja função desconhecida 
depende de mais de uma variável. 

6. Solução geral de uma equação diferencial: é a solução obtida com a resolução 
inicial da equação diferencial e que contém a constante de integração; 

7. Solução particular de uma equação diferencial: é a solução obtida atribuindo um 
valor para a constante de integração, geralmente calculado através de uma condição inicial. 


u LEITURA 

Para se aprofundar nos conceitos de equações diferenciais, suas classificações e aplicações, 
recomendamos que você leia as obras listadas nas referências bibliográficas do capítulo, em 
especial: 

BASSANEZI, R.C. Equações Diferenciais Ordinárias - Um Curso Introdutório. Coleção 
BC&T - Textos Didáticos - Volume 1 - Universidade Federal do ABC. Capítulos 1 e 2. Dis- 
ponível em: http://gradmat.ufabc.edu.br/disciplinas/listas/iedo/notasdeaulas/equacoes-di- 
ferenciais-ordinrias-rodney.pdf [Acessado em 26/05/2015] 

BOYCE, W.E e DiPRIMA, R.C. Equações Diferenciais Elementares e Problemas de Valo- 
res de Contorno. 9 a Ed. Rio de Janeiro: LTC, 201 4. [Ler o capítulo 1] 

DIACU, Florin. Introdução a Equações Diferenciais - Teoria e Aplicações. LTC: Rio de 
Janeiro. 2004. [Ler o capítulo 1 , páginas 01 - 12] 

EDWARDS JR, C.PI. e PENNEY, D.E. Equações Diferenciais Elementares com Problemas 
de Contorno. 3 a Ed. Rio de Janeiro: Prentice-Hall do Brasil, 1 995. [Ler o capítulo 1 , páginas 
02-15] 

NAGLE, R. Kent; SAFF, Edward B.; SNIDER, Arthur David. Equações Diferenciais. 8 ed. 
Pearson: São Paulo. 201 2. [Ler o capítulo 1 , páginas 01-19] 

ROSA, M.S.R. Equações Diferenciais. Instituto de Matemática. Universidade Federal do 
Rio de Janeiro, 2009. Capítulo 1. Disponível em: http://www.dma.im.ufrj.br/~rrosa/dvifiles/ 
aposfila-ed.pdf [Acessado em 26/05/201 5] 

ZILL, D.G. e CULLEN, M.R. Equações Diferenciais. Volume 1. 3 a Ed. São Paulo: Makron 
Books, 2001 . [Ler o capítulo 1], 
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Métodos de 
Resolução 
de Equações 
Diferenciais de 
Primeira Ordem 



Neste segundo capítulo vamos discutir as principais técnicas de resolução das 
equações diferenciais ordinárias, EDO, de primeira ordem, detalhar cada caso 
particular e como proceder em cada um. 


^7 OBJETIVOS 


■ Entender o processo de resolução de uma EDO de ordem 1 via variáveis separáveis. 

■ Entender o processo de resolução de uma EDO de ordem 1 via equações homogêneas. 

■ Entender o processo de resolução de uma EDO de ordem 1 via equações exatas. 

■ Entender o processo de resolução de uma EDO de ordem 1 via equações lineares. 
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2.1 Considerações iniciais 


Antes de partir para técnicas particulares de resolução de equações diferen- 
ciais de primeira ordem com condições especiais, é importante discutir alguns 
pontos. 

Primeiro, será que podemos garantir a existência de solução de uma equação 
diferencial? Melhor ainda, será que podemos garantir de essa solução é única? 

Sim, podemos, isso graças a um teorema enunciado e provado por Charles 
Émile Picard (1856-1941), de acordo com ZILL e CULLEN (2001), que associa a 
existência de uma única solução a continuidade das funções e de suas deriva- 
das a existência de uma região R no plano. Não vamos verificar tais condições 
nos casos discutidos neste livro, mas não é difícil avaliar a continuidade das 
funções e de suas derivadas, nem tampouco determinar uma região R válida de 
acordo com o comportamento das variáveis. Apenas a título de melhor entendi- 
mento do teorema, vejamos um exemplo de como obter a região R. 

Teorema da existência de uma única solução: Charles Émile Picard (1 856-1941) 

Seja R uma região retangular no plano xy definida por a<x<bec<y<d que con- 
tém um ponto (x 0 ,y 0 ) em seu interior. Se f(x,y) e ^ são contínuas em R, então existe um 
intervalo / centrado em x 0 e uma única solução y(x) definida em / que satisfaz o problema 
de valor inicial. 


C7 EXEMPLO 

Seja a equação diferencial fíZ. = 
dx 

Neste caso temos que f(x,y) = ^/y , que será contínua para qualquer número real x e para 
qualquer valor de y > 0. Calculando a derivada parcial de f(x, y) em relação a y, temos: 
df(*,y) _êy/y _3 rj 
Õy 8y 2 ' y 

Da mesma forma, a derivada parcial será contínua para qualquer x pertencente aos reais 
e para valores não negativos de y. Assim R = {xe R e ye R + } 

Bem, é oportuno enfatizar que o teorema citado anteriormente é importante no sentido de 
tranquilizar os usuários de equações diferenciais no sentido de garantir que existirá uma solu- 
ção única dada uma condição inicial para o problema, mas infelízmente o teorema não garante 
que a descoberta de tal solução será simples. Ainda, não nos orienta como chegar até ela 
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É por isso que na evolução histórica do estudo das equações diferenciais foram cons- 
truídos muitos modelos de resolução, conforme discutido em BOYCE e DiPRIMA (2014), 
baseados em características particulares das equações, ou seja, dependendo da forma como 
a equação é apresentada, temos um, ou em alguns casos mais de um, procedimento padrão 
para sua resolução. 

Preparado? Vamos lá, então! Começaremos com uma família de equações diferenciais 
chamada de EDO de Variáveis Separáveis. 


2.2 Método das Variáveis Separáveis 

Para melhor compreensão do método das variáveis separáveis, vamos iniciar 
nossa discussão com uma equação diferencial de primeira ordem bem simples 
e básica. Acompanhe os exemplos a seguir e procure associar com o discutido 
no capítulo 1. 


XJ EXEMPLO 

Seja a equação diferencial = ^ 
dx 

Note que nesta equação diferencial, temos uma função que depende apenas da variável 
x, e com isso é simples isolar de um lado da igualdade apenas termos associados a y e do 
outro lado termos associados apenas a variável x, da forma: 

-^ = Vx=>l-dy = Vx-dx 

dx 

Como cada lado depende de apenas uma variável, é fácil resolver a equação diferencial, 
encontrar a sua solução geral, utilizando díretamente integrais, da forma: 
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Jl-dy = JVx-dx 
y+C t = Jx 1/2 dx 

y + C, = — x 3/2 +C, 
J 1 3 2 

y+Ci=^+c 2 

+c z~ c i 

y^sí^ + C 


XJ EXEMPLO 

dv x 2 

Seja a equaçao diferencial ~dL = — Resolvendo de maneira similar ao exemplol , temos: 
dx 2y 

^ = — =>2yvdy = x 2 -dx 

dx 2y 

Í2ydy = jx 2 -dx 

y 2 +C. =-x 3 +C, 

1 3 2 

y 2 =ix 3 +C 

3 

y=Jf^ 

É possível notar que nos dois exemplos foi possível separar uma função que depende 
apenas de ajunto com a diferencial dx da variável y que ficou junto com a diferencial dy, fa- 
vorecendo o uso direto das integrais. 

Essa família de EQO de primeira ordem são solucionadas, portanto, pela técnica de Se- 
paração de Variáveis. 

Segundo ZILL e CULLEN (2001), sempre que for possível separar as variáveis x e y em 
funções h(y) e g(x) de forma que: 

dy _ g(x) 

dx h(y) 

a equação diferencial é separável ou tem suas variáveis separáveis, e a solução geral 
será igual a: 
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dy _ g(x) 
dx h(y) 
h (y)-dy=g(x)-dx 
|h(y)'dy + C 1 =|g(x)dx + C 2 

Assim, obtemos a solução geral para o caso de variáveis separáveis: 

jh(y)dy = jg(x)dx + C 
onde C = C,+C 2 . 

Assim, se a equação diferencial for de variáveis separáveis, basta identificar as funções 
h(y) e g(x) e substituir na solução geral. 


XJ EXEMPLO 

Resolva a equação diferencial e y y'=cos x 

Podemos representar a equação diferencial da seguinte forma: 
dy 

e y — = cosx=>e y dy = cosxdx 
dx 

É imediato notar que h(y)=e / e que g(x)=cos x, assim a solução geral será: 

je y dy = Jcosxdx + C 
e y =senx + C 


XJ EXEMPLO 

Resolva a equação diferenciai x 2 y . -2xy 2 =0 ■ 

dx 

De início, vamos determinar as funções h(y) e g(x). Passando um dos termos para o lado 
direito da igualdade e passando dx multiplicando, temos: 

x 2 y ■ — -2xy 2 =0 

dx 

x 2 y -dy = 2xy 2 - dx 

Note que para gerar a função h(y) temos que dividir o termo à esquerda por x 2 e para 
gerar a função g(x) temos que dividir o termo à direita por y 2 . Desta forma, basta dividir a 
equação diferencial por^y 2 : 
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x 2 y 

x 2 y 2 


dy = 


2xy 2 

x 2 y 2 


■dx 



Assim, chegamos a: h(y) = — e g(x) = — 

y x 

Então, nota-se que a equação diferencial é de variáveis separáveis e podemos substituir 
as duas funções encontradas na solução geral, obtendo: 

J— dy = f-dx + C 

J y J X 

Lny =2Lnx+C 

gLny _g2 Lnx+C 
y — — g2.Ln x t gC 

Na expressão, como e é uma constante (e = 2,72), podemos escrever da forma: 
y = Ce 2Lnx 


EXEMPLO 

dy y 

Seja uma EDO de primeira ordem ^=— com uma condição inicial y(2) = 1. Resolva via 
separação de variáveis. 

Temos que h(y)=1/y e g(x)=1/x. Assim: 

f-dy= f-dx+C 

j y j x 

Lny=Lnx+C 

pLny _pLnx+C 
y = C-x 

Com a condição y(2)= 1 , temos: 

y = Cx=>l = C-2=>C = — 

2 


Assim, a solução particular da EDO é: y 


x 

2 


CAPÍTULO 2 ■ 39 



Vejamos um exemplo aplicado envolvendo equações com variáveis 
separáveis. 


C7 EXEMPLO 

[adaptado de BOYCE e DIPRIMA(201 4)\. Um exemplo clássico de aplicação das equações 
diferenciais é o estudo decaimento radioativo em função do tempo. 0 decaimento radioativo 
ocorre quando os núcleos de isótopos instáveis são rompidos por conta da instabilidade atô- 
mica. Como ilustração, temos o Urânio, cujo isótopo U-238 é desintegrado até a produção 
do Tório, cujo isótopo Th-234 é desintegrado até a produção do Protactínio, cujo isótopo 
Pa-234 é desintegrado até a produção do Chumbo, que é estável e com isso o decaimento 
encerra. Assim, dependendo do tempo de decaimento, novos isótopos podem ser formados. 

Para representar o decaimento radioativo de uma quantidade inicial Q em função do 
tempo t, temos a seguinte representação diferencial: 

?--*<*> 

Assim, suponhamos um determinado isótopo que no início do processo de decaimento 
(t=0) tinha massa igual a Q=500g e que após 10 minutos ( t=10 ) sofreu uma redução de 
1 5% em sua massa ( Q=425g ). Considerando tais condições iniciais, qual a equação que 
determina a massa do isótopo em um tempo qualquer? 

Vamos inicialmente encontrar a solução geral da equação diferencial, utilizando a técnica 
de separação das variáveis. 


ÜQ=-k. 

dt 
dQ 
Q 

í- = í- 

J Q J 


= -kdt 


kdt 


Ln Q = -kt + C 
Q = Ce -kt 


Para determinar a solução particular, vamos considerar a condição inicial de que a massa 
Q era igual a 500g com t=0. 


Q = Ce kt => 500 = Ce~ k0 => 500 = Ce° 500 = C 
Portanto a solução da equação diferencial de decaimento é: 

Q = 500e kt 
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Agora, para o isótopo em questão, qual a constante de decaimento k Para encontrar o 
valor de k temos que utilizar a outra condição do problema, uma massa Q = 425g para um 
tempo t = 10. Temos: 

425 

Q = 500e -kt => 425 = 500e- lok => = e- 10k => 

500 

=> 0,85 = e- 10k => ln(0,85) = -10k => 

=> -0,162518929 = -10k=>k = 0,0163 

Assim, temos a solução final igual a: Q = 500e~° ,0163t . 


2.3 Método das Equações Homogêneas 

Você pode ter notado que se uma equação diferencial ordinária de primeira or- 
dem for separável quanto às variáveis, fica bem facilitada a sua resolução quan- 
to a solução geral, correto? 

Sim, mas será que toda EDO de primeira ordem é separável? 

Infelizmente não. Mas isso pode não ser exatamente um problema. Existe 
uma família de EDO de primeira ordem que podem ser reescritas de maneira a 
serem separáveis. São as equações diferenciais homogêneas. 

A homogeneidade da equação diferencial está associada a homogeneidade 
das funções presentes na equação. 

Como ilustrado por ZILL e CULLEN (2001), uma função é homogênea quan- 
do satisfaz a seguinte relação: Se uma função f(x,y) satisfaz a relação f( tx,ty)=t n 
■ f(x,y)é dita ser homogênea de grau n. 


U EXEMPLO 

Sejaafunção f(x,y) = x 3 +2xy 2 -2y 3 . 

Substituindo x por tx e y por ty, temos: 

f(tx,ty) = (tx) 3 + 2tx(ty) 2 -2(ty) 3 
f(tx,ty) = t 3 x 3 + 2txt 2 y 2 - 2t 3 y 3 
f(tx,ty) = t 3 x 3 + 2t 3 xy 2 - 2t 3 y 3 
f(tx,ty) = t 3 (x 3 + 2xy 2 — 2y 3 ) 

f(tx,ty) = t 3 f(x,y) e, portanto, é homogênea de grau 3. 
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C7 EXEMPLO 


Seja a função f(x,y) = x 3 +y 3 +5. . 

Substituindo x por tx e y por ty, temos: 
f(tx,ty) = (tx) 3 +(ty) 3 +5 
f(tx,ty) = t 3 x 3 +t 3 y 3 +5 
f(tx,ty) = t 3 (x 3 + y 3 +5/t 3 ) 

f(tx,ty)?tt 3 f(x,y) e portanto não é homogênea. 

Sobre as funções homogêneas existe uma forma rápida de avaliar a satisfação da 
relação. Se a soma dos graus das variáveis de cada termo for sempre o mesmo, a função 
será homogênea. Note que no exemplo 7, nos três termos a soma dos graus das variáveis é 
igual a 3 (x 3 y 0 / 1 y 2 e x° y 3 ) e com isso a equação é homogênea de grau 3. Já no exemplo 8, 
devido a presença da constante 5, nem todos os termos possuem o mesmo grau, portanto 
não homogênea. 

Supondo, então, uma função homogênea, conforme apontado por EDWARDS JR e PEN- 
NEY (1 995), podemos reescrever a função original em termos de duas novas funções homo- 
gêneas, mas de grau zero, da forma: 


f(x,y) = y n f^,i] 

onde n representa o grau de homogeneidade da função. 

Com os resultados anteriores e considerando ainda o exposto em ZILL e CULLEN 
(2001), uma equação diferencial do tipo: 
g(x, y) dx + h(x, y) dy =0 

será homogênea se os coeficientes g e h forem funções homogêneas de mesmo grau e 
se isso ocorrer significa que a EDO pode ser separada pelas variáveis. 

Para encontrar a solução geral de uma equação diferencial de primeira ordem homogê- 
nea pelo método de separação de variável, é necessário fazer uma mudança de variável, que 
pode ser executada de duas formas diferentes: 

Primeira forma: substituir/ por ux,y=ux, o que gera como consequência que a diferen- 
cial dy será obtida pela derivada do produto ux, ou seja dy=u.dx+x.du. 

Segunda forma: substituir x por uy, x=uy, o que gera como consequência que a diferen- 
cial d*" será obtida pela derivada do produto uy , ou seja dx=u.dy+y.du. 
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XJ EXEMPLO 

Para verificar a validade das formas anteriores, vamos considerar a forma geral de uma equa- 
ção diferencial homogênea: 
g(x, y) dx + h(x, y) dy = 0 

Fazendo uso da segunda forma, x = uy e dx = u ■ dy + y ■ du, temos: 
g(uy, y)(u ■ dy + y ■ du) + h(uy, y)dy = 0 

Aplicando os resultados sobre funções homogêneas, podemos transformá-las em novas 
funções homogêneas de grau zero, da forma: 

g(uy,y) = y n g(u, 1) 
h(uy, y) = y n h(u, 1 ) 

Substituindo, temos: 
y n g(u,l)(u.dy + y.du) + y n h(u, l)dy = 0 
y n [g(u,l)(u.dy + y.du) + h(u, l)dy] = 0 
g(u,l)(u.dy + y.du) + h(u, l)dy = 0 
u.g(u,l)dy + y.g(u,l)du + h(u,l)dy = 0 
[u.g(u,l) + h(u,l)]dy + y.g(u,l)du = 0 
[u.g(u,l) + h(u,l)]dy = -y.g(u,l)du 

dy = M du 

y u-g(u,l) + h(u,l) 


Assim a equação diferencial está devidamente separada quanto às variáveis e os resul- 
tados discutidos no item 2.1 podem ser aplicados na resolução da EDO. O mesmo resultado 
pode ser obtido pela primeira forma, como ilustrado por ZILL e CULLEN (2001 , p 55). 


XY EXEMPLO 

Seja a equação diferencial dy_ x 2 +y 2 u ti|j Z a n do os conceitos de equação diferencial ho- 
dx x 2 -xy 

mogênea. 

Vamos verificar a homogeneidade da equação. Temos: 
g(x,y) = x 2 +y 2 =x 2 y° +x°y 2 
h(x,y) = x 2 -xy = x 2 y° +x 1 y 1 
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Assim, tanto a função g quanto a função h sao homogêneas de grau 2. 

Fazendo a substituição y = ux e escrevendo dy como dy = u ■ dx + x ■ du, temos: 

g(x,y)dx + h(x,y)dy = 0 

(x 2 + u 2 x 2 )-dx+(x 2 -x 2 u)(udx + xdu) = 0 

(x 2 + u 2 x 2 )-dx+(x 2 -x 2 u)udx+(x 2 -x 2 u)xdu = 0 

(x 2 + u 2 x 2 )-dx+(ux 2 -x 2 u 2 )-dx+(x 3 -x 3 u)du = 0 

(x 2 +u 2 x 2 +ux 2 -x 2 u 2 )-dx+(x 3 -x 3 u)du = 0 

(x 2 + ux 2 ) ■ dx + (x 3 - x 3 u) ■ du = 0 

x 2 (l+u)-dx+x 3 (l-u)du = 0 


Dividindo a expressão por x 2 , temos: 
(1 + u 2 ) ■ dx + 2ux ■ du = 0 


Dividindo a expressão por 2u, temos: 
(1 + u) 

(1-u) 


■dx+xdu =0 


Multiplicando a expressão por 

i dx-^.du.O 
x (1+u) 


2u 


-(i +u q 


temos: 


Assim, transformamos a equação diferencial original em uma equação diferencial com 
variáveis separadas, que pode ser resolvida pela técnica discutida no item 2.1. 

Utilizando a técnica anterior, temos que integrar ambos os lados, mas note que a integral 
referente a u trará alguma dificuldade. Assim, antes de integrar, vamos reescrever a fração 
original através de uma alteração por frações parciais, da forma: 

(l-u) A B 

— 1 

(l + u) (l + u) (l + u) 


Para que a igualdade seja verdadeira, é necessário que: 
A + B=1-u=>A=1-u-B 
Supondo B = 2, temos: 

A= 1 - u - B = 1 - u - 2 = - 1 - u 


Substituindo na equação original, temos: 
du = 0 


— -dx + 
x 


— dx + 


-l-u 2 
l + u l + u 
2 


-1 + - 


1 + u 


■du = 0 


Agora fica mais fácil o cálculo das integrais. Prosseguindo com o cálculo, temos: 
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f— dx- fldu + 2f— !— du =0 
x ■* ■* 1 + u 

Lnx - u + 2 ■ Ln(l + u) = C 
-u + Lnx + Ln(l + u) 2 = C 
-u + Ln[x-(l + u) 2 ] = C 


Neste momento temos que lembrar a mudança de variável feita no início do problema, ou 

seja, y = yx, o que significa dizer que u = — 

x 

Assim, temos: 


= C 




f 

2 

+Ln 

x- 

l+y 


X 


V x; 


+ Ln 

X- 

^x + y' 

f 

X 


V x 


-- + Ln 

X 


(x+y) 

2 “ 


X 2 


+ Ln 

"(x+y) 2 ' 

= < 

X 


X 



=c 


Supondo ainda que a constante C possa ser escrita como Ln C, temos: 


+ Ln 


( x+ y ) 2 


Ln 


Ln 


(x+y ) 2 

X 

(x+y ) 2 

X 


= LnC 


-LnC = — 


Ln 


Ln 


(x+y ) 2 i 

x C 

(x+y ) 2 

Cx 
\2 


(x+y ) 2 z 

— — çx 

Cx 


Temos fínalmente uma solução geral para a equação diferencial dada por: 

2 l 

(x + y) =Cxex 
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2.4 Método das Equações Exatas 


O termo equações exatas está relacionado com o comportamento das derivadas 
de segunda ordem que podemos obter supondo que as funções associadas à 
equação diferencial sejam do tipo f(x,y). 

Para melhor compreensão, vamos retomar a representação de uma equação 
diferencial conforme exposto no tópico 2.2, ou seja: 
g(x,y)dx + h(x,y)dy = 0 

A equação diferencial ordinária de primeira ordem será exata se e somente 
se ocorrer a seguinte igualdade, como explicitado por ZILL e CULLEN (2001): 

dg( x >y) _ d h ( x .y) 
õy õx 


Independentemente de ser uma equação diferencial separável ou homogê- 
nea, o fato de ser exata acaba por indicar um caminho para o cálculo da solução 
geral da EDO. De acordo com EDWARDS e PENNEY (1995), a solução ocorre em 
função de um teorema que garante o seguinte resultado: 

Se uma equação diferencial de primeira ordem é exata, deve existir uma fun- 
ção f(x,y) tal que 2^lZ} = g( X; y) e t!tElZl = h(x,y) , tendo como solução geral 
õx õy 

f(x,y)=C, com C uma constante real qualquer. 

Vejamos alguns exemplos: 


XJ EXEMPLO 

Resolver a seguinte equação diferencial: 

(5y - 2x) dx + (5x + 3y 2 ) dy = 0 

Vejamos inicialmente se a equação diferencial é ou não exata. Temos: 

, . 5g(x,y) 

g(x,y) = 5y-2x=> M ’ =5 

õy 

, , , , 5h(x,y) 

h(x,y) = 5x + 3y 2 => ^ ' = 5 

ox 
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Como as derivadas parciais são iguais, a equação diferencial é exata. Sendo assim, deve 
existir uma função f(x,y) tal que ínicíalmente: 


df(x,y) 

dx 


=g( x -y)= 


df(x,y) 

ax 


= (5y-2x)=>af(x,y)=(5y-2x)dx 


Calculando a integral em ambos os lados, em relação ax, temos: 

JSf(x,y)dx = J(5y-2x)dx 
f(x, y) = 5xy - x 2 + u (y) 

O aparecimento de uma função u(y) é por conta da derivada parcial do lado esquerdo. 

df(x,y) 

Continuando, não podemos esquecer da outra condição, = h(x,y) Contudo, 

õy 

neste caso já temos um resultado parcial para f(x,y) e utilizaremos agora, fazendo: 

õf(x,y) . õÍ5xy-x 2 +u(y)) 

-^ = hx,y^ J t — = 5x + 3y 2 

9y õy 

5x-0 + u'(y) = 5x + 3y 2 =>u(y) = 5x + 3y 2 -5x + 0 = 3y 2 
Integrando u’(y) em função de y, temos: 


u (y) = 3y 2 ^Ju(y)dy = j3y 2 dy^.u(y) = y 3 +C 1 

Bem, pense um pouco, o que fazer com o valor de u(y) Se você pensou em substituir na 
expressão anterior de f(x,y), está com razão, é exatamente isso que devemos fazer: 


f ( x -y) = 5xy-x 2 +u(y)^>f(x,y) = 5xy-x 2 +y 3 +C 1 

Pergunta: essa é a solução geral da equação? Não, não é. Retornando ao conceito de 
equações exatas, deve existir uma função f(x,y) que satisfaz as condições anteriores e dada 
essa função, a solução geral é da forma f(x,y)=C. Assim, escrevemos a solução como colo- 
cado abaixo: 


C = 5xy-x 2 +y 3 -i-Cj 
5xy-x 2 +y 3 -C = 0 


Se desejar verificar se a solução está correta, basta derivar a f(x,y) em relação a cada 
variável e verificar se retorna para as funções g e h, da forma: 
õf (x,y) _ õ(5xy -x 2 + y 3 + Cj ) 


õx õx 

õf (x,y) _ s(5xy-x 2 + y 3 + Cj ) 
õx õy 


= 5y-2x + 0 = 5y-2x = g(x,y) 

= 5x-0 + 3y 2 =5x + 3y 2 =h(x,y) 


Portanto a solução está correta. 
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C7 EXEMPLO 

, -2xy 3 -4 

Resolver a seguinte equação diferencial: y = 3 X 2y2 _ 2 


Reescrevendo temos: 

— = 2xy — -=>Í3x 2 y 2 ~2)dy = (-2xy 3 -4)dx 
dx 3x 2 y 2 -2 V y ’ 

Colocando na forma inicial, ficamos com: 

(3x 2 y 2 -2)dy+(2xy 3 + 4)dx = 0 

Vejamos inicialmente se a equação diferencial é ou não exata. Temos: 

g( x »y)= 2x y 3 + 4 ^ =6xy 2 

õy 

h(x,y) = 3x 2 y 2 -2^> C ' ih ^ X ’ y) =6xy 2 


Como as derivadas parciais sao iguais, a equaçao diferencial é exata. Sendo assim, deve 
existir uma função f(x,y) tal que inicialmente: 


5f(x,y) 


dx 


=g( x -y)= 


5f(x,y) 


0X 


= 2xy 3 +4=>5f(x,y) = (2xy 3 +4)dx 


Calculando a integral em ambos os lados, em relaçao a x, temos: 

f( x .y)=J( 2 x y 3 + 4 )d x 

f(x,y) = x 2 y 3 +4x + u(y) 

Fazendo agora: 

^=h(„)» í í í ’^l=3,V-! 

õy õy 

3x 2 y 2 +0 + u'(y) = 3x z y 2 -2=>u (y) = 3x z y 2 -2 + 3x z y 2 +0 = -2 
Integrando u’(y) em função dey, temos: 
u (y) = -2=> Ju (y)dy = J-2dy=>u(y) = -2y + C 1 

Bem, pense um pouco, o que fazer com o valor de u(y)7 Se você pensou em substituir na 
expressão anterior de f(x,y), está com razão, é exatamente isso que devemos fazer: 
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f(x,y) = x 2 y 3 +4x + u(y)^>f(x,y) = x 2 y 3 +4x-2y + C 1 


Assim, escrevemos a solução como colocado abaixo: 

C = x 2 y 3 + 4x-2y + C 1 
x 2 y 3 +4x-2y-C = 0 


2.5 Método das Equações Lineares. 


Como definido no capítulo 1, classificamos uma equação diferencial como li- 
near se todos os coeficientes das derivadas são funções apenas da variável inde- 
pendente x, e ao mesmo tempo a variável dependente ye suas derivadas sejam 
de primeira ordem. 

Como apresentado por ZILL e CULLEN (2001), uma equação diferencial li- 
near pode ser representada de maneira geral pela seguinte expressão: 


(x)£ + a 0 (x)y = g(y) 


Dividindo a equação diferencial em sua forma geral por a 1 (x) temos uma 
maneira mais simplificada de representar uma EDO linear: 

dy t a o( x ) .. g( x ) 
dx a 1 (x) a 1 (x) 


^ + P(x)y = f(x) 
dx 


Nesta representação, as funções P(x) e f(x) são contínuas em uma dada re- 
gião, lineares e dependentes apenas de x. Ainda, podemos reescrever a equação 
diferencial da seguinte forma: 
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^ + P(x)y = f(x) 
dx 

^ = -P(x)y+f(x) 
dx 

dy = [-P(x)y + f (x)] dx 
dy + [p(x)y-f(x)]dx = 0 


Como afirmado em ZILL e CULLEN (2001), se a equação diferencial é line- 
ar, como a colocada anteriormente, é possível determinar uma função p/y/que 
seja um coeficiente dos dois termos da equação diferencial de maneira a não 
modificar a igualdade, ou seja: 

p(x)-dy + p(x)-[p(x)y-f(x)]dx = 0 


Note que, representada desta forma, temos uma equação diferencial exata 
e vale o seguinte resultado: 

8p(x) _ g{p(x)[P(x)y-f(x)]} 
õx õy 

8p(x) _ 8[|i(x)P(x)y] 8[p(x)f(x)] 

õx õy õy 

íd^ = n( x )P(x)-0 
MíU(x)dx 

p(x) 

Note que saímos de uma equação diferencial linear para uma equação di- 
ferencial exata e no final em uma equação diferencial separável. Aplicando a 
integral dos dois lados, temos: 

/i^) dt ‘ = í p(x)dx 

Ln|i(x) + C 1 = Jp(x)dx 

Lnp(x) = Jp(x)dx-C 1 


Como não conhecemos a função P, a integral do lado direito não tem como 
ser resolvida, mas isso não impede de determinarmos a forma final para a fun- 
ção u (x): 
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e Lnn(x) =e jP(x)dx-C 1 

n(x) = C-J P(x)dx 


Assim, a função u (x) representa o Fator de Integração da equação diferen- 
cial e, como bem explicado em ZILL e CULLEN (2001), a constante Cnão afeta 
o resultado final. Assim, podemos escrever o fator integrante u (x) da forma: 

pi(x) = e/ p ( x > dx 


Você deve estar se perguntando, qual a importância do fator de integração? 
A importância é que se multiplicarmos a equação diferencial linear por u 
(x), teremos como resultado uma equação separável e integrável diretamente, 
por isso o nome fator de integração. 

Para utilizar o fator de integração, basta identificar a função P(x). 


EXEMPLO 

Resolver a equação diferencial xy' — 3y = x 5 e x ■ 

Reescrevendo a equação diferencial, temos: 

x— -3y = x 5 e x 


Dividindo a equação por x: 

dy 3y x 5 e x 
dx x x 

— - — y = x 4 e x 
dx x 

Note que a equação diferencial já está na forma padrão de uma equação line- 

q 

ar e é possível verificar que p(x) = — o que resulta no seguinte fator integrador: 


p(x) = e 


r 3 j-dx 

= £»•' x — p x — p _ 3Ln(x) _ , 
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De posse do fator integrador, multiplicamos na equação diferencial: 


X 3^--yx 3 =X 3 X 4 e x 
dx x 

x~ 3 — -3yx _J =xe x 
dx 


Note que o termo do lado esquerdo pode ser escrito como a derivada do pro- 
duto entre y e u(x) em relação a x, ou seja: 

dx dx 


Assim, ficamos com: 


d(yx 3 ) 

— - = xe x 


Integrando dos dois lados, 
yx~ 3 = jxe x dx 

Resolvendo a integral do lado direito por partes, temos: 

u = x => du = dx 
dv = e x dx => v = e x 
Ju-dv = u-v-Jv-du 

Jxe x dx = xe x -Je x dxjxe x dx = xe x -(e x +C) 

Assim, a solução geral da equação diferencial será: 
yx ~ 3 =x e x -e x -C 

Dividindo a expressão por x~ 3 : 
y = x 4 e x -x 3 e x -x 3 C 

y = x 4 -e x -x 3 (e x -C) 

Que é a solução geral da equação diferencial 
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u EXEMPLO 

[Adaptado de BOYCE e DiPRIMA (2014)] Suponha um corpo de massa m, posicionado no 
alto de uma torre, em repouso. Se esse corpo cair em uma trajetória ortogonal ao solo, sen- 
do afetado de forma positiva pela ação da gravidade gf e de forma negativa por uma força 
resístíva proporcional a sua velocidade v de queda, determine uma expressão que aponte a 
velocidade de queda em função do tempo t. 

Bem, temos duas forças atuando sobre o corpo. Uma força positiva F 1 = m.g que contri- 
bui para a aceleração do corpo no sentido do solo e uma outra força F 2 = k.v que desacelera 
a queda em função da resistência encontrada na trajetória, onde k é uma constante. Pela se- 
gunda lei de Newton, a somatória das forças F e F 2 é igual ao produto da massa m do corpo 
pela sua aceleração, e a aceleração por sua vez é a derivada da velocidade 1 / pelo tempo t. 
Assim, temos a seguinte equação diferencial: 

, dv 

msr-k-v=m — 
s dt 

Reescrevendo a equação diferencial, temos: 
dv 

m — +k-v = m-sr 
dt S 

Dividindo por m, ficamos com: 
dv k 

+ — v = g 


Observe que a forma finai da equação diferencial está linearizada e a função P(t) — — 
e, portanto o fator integrador será: 

r k . k r. k 

/ \ [— 1 dt — [dt — t 

p(t) = e J m =em J =em 


Multiplicando na equaçao diferencial, temos: 

k j ki k 

— t dv — t k — t 

em he m — v = ge m 

dt m 
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Note que o lado esquerdo pode ser escrito como a derivada 


dv 


/ 


A equação diferencial passa a ser escrita como: 

d 


( 

vem 


dv 


= gem dt 


Integrando dos dois lados, temos: 


ve m 




em + 


C 


em 

Dividindo por 

k 

m _ — t 
v = s ; — hCe m 
5 k 


Agora, você deve lembrar que no instante t=0 o corpo estava em repouso (v=0). Assim: 

k 

~ ^ — 0 m ^ ^ m 

0 = g — + Ce m =>0 = g — + C=>C = -g — 
k k k 


Substituindo o valor da constante na equação diferencial, temos a solução geral procura- 
da para a velocidade da queda do corpo em função do tempo, dada por: 

k 

m m — t 

v = g g — e m 

k k 


1-e m 


\ 


y 
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Ey CONEXÃO 

O estudo das técnicas de resolução de equações diferenciais ordinárias de primeira ordem 
realizado neste capítulo não leva em consideração a ocorrência de equações especiais ou 
particulares. As técnicas procuram identificar características que levem a um procedimento 
de resolução, que se aplique a qualquer equações com as mesmas características (separá- 
veis, homogêneas, exatas ou lineares). Contudo, existem algumas equações diferenciais bem 
específicas que possuem métodos próprios de resolução, como as equações diferenciais de 
Jacques Bernoulli (1 654-1 705), as equações de Jacob Ricattí (1 676-1 754) e também de 
Alex Clairaut (1713-1765). Recomendamos que leia a referência abaixo para ampliar seus 
conhecimentos sobre essas equações diferenciais: 
http://www.feg.unesp.br/~ernesto/guiaedo/Tcc.pdf 
Ver capítulo 2, páginas 26 a 32. 


^ CONEXÃO 

Da discussão realizada no capítulo, você deve ter percebido que em alguns casos é possível 
fazer uma transformação na equação diferencial original a fim de coloca-las em um padrão 
que atenda a uma determinada técnica de resolução. Para saber mais sobre as transforma- 
ções e ou substituições para simplificar a solução de uma EDO, recomendamos que leia a 
referência abaixo e assista ao vídeo indicado para ampliar seus conhecimentos sobre o uso 
de substituições na resolução de EDO de primeira ordem: 

http://www.feg.unesp.br/~ernesto/guiaedo/Tcc.pdf [Ver capítulo 2, páginas 1 0 a 1 4], 

Vídeo parte 1 : 

https://www.youtube.com/watch?list=PL444C47ACE456066B&v=el8Pwa6ztAQ 
Vídeo parte 2: 

https://www.youtube.com/watch?v=VGslyKQInyk&index=2&list=PL444C47ACE456066B 


EP ATIVIDADES 

05. Resolva as equações diferenciais abaixo através de separação de variáveis. Verifique em 
cada caso se de fato a equação é separável, encontre a solução geral e, quando for o caso, a 
solução particular em função da condição inicial imposta. 

a) y 1 (x 2 + 9) = xy, com y(4) = 1 0 

b) y' = y~ 3 sen x 
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06. Resolva as equações diferenciais abaixo através da técnica para equações homogêne- 
as. Verifique em cada caso se de fato a equação é homogênea, encontre a solução geral. 

a) X y 3 ^ = 2y 4 + x 4 

dx 

b ) dy _ y 2 - x 2 
dx 2xy 

07. Resolva as equações diferenciais abaixo através da técnica para equações exatas. Veri- 
fique em cada caso se de fato a equação é exata e encontre a solução geral. 


a) 


/ = 


2x-3y "l 
2y-3x J 


b) 


ydy -xdx 
x 2 + y 2 x 2 + y 2 


08. Resolva as equações diferenciais abaixo através da técnica para equações lineares. Ve- 
rifique em cada caso se de fato a equação é linear, encontre a solução geral e, quando for o 
caso, a solução particular em função da condição inicial imposta. 

a) y’ = e 3x + 5y 

b) dy = x 2 (3y + 1) dx 


REFLEXÃO 

Todas as técnicas trabalhadas neste capítulo impõem uma condição importante sobre a 
equação diferencial, que sejam de primeira ordem, ou seja, a derivada de maior grau é a 
derivada primeira. Contudo, podemos refletir sobre o seguinte: se existem situações práticas, 
como apresentadas no capítulo 1, que exigem derivadas de segunda ordem ou de ordem 
superior, será que é possível determinar métodos especiais para a resolução desse tipo de 
equação diferencial? 


^ CONCEITO 

Alguns termos citados neste capítulo são de especial importância e aparecerão nos capítulos 
futuros, sendo assim reforçamos seus conceitos neste glossário. 

1. Equações diferenciais separáveis: são equações diferenciais onde podemos sepa- 
rar de um lado da igualdade funções que dependam apenas de x e dx e do outro lado da 
igualdade funções que dependam apenas deye dy ; 


50 ■ CAPÍTULO 2 




2. Equações diferenciais homogêneas: são equações diferenciais onde a soma do 
grau das variáveis de cada funções sempre resultam no mesmo grau total; 

3. Equações diferenciais exatas: são equações diferenciais que possuem em sua for- 
mulação função do tipo f(x,y) de maneira que a derivada parcial de f em relação ay é igual a 
derivada parcial em relação ax. 

4. Equações diferenciais lineares: todos os coeficientes das derivadas são funções 
apenas da variável independente x, e ao mesmo tempo a variável dependente y e suas deri- 
vadas sejam de primeira ordem; 

5. Fator integrante: é uma função que multiplicada por uma equação diferencial a torna 
solucionável com a aplicação direta de integrais. 


ET LEITURA 

Para se aprofundar nos conceitos e técnicas de resolução das equações diferenciais de 
primeira ordem, recomendamos que você leia as obras listadas nas referências bibliográficas 
do capítulo, em especial: 

BOYCE, W.E e DiPRIMA, R.C. Equações Diferenciais Elementares e Problemas de Valo- 
res de Contorno. 9 a Ed. Rio de Janeiro: LTC, 2014. [Ler o capitulo 2] 

EDWARDS JR, C.H. e PENNEY, D.E. Equações Diferenciais Elementares com Problemas 
de Contorno. 3 a Ed. Rio de Janeiro: Prentice-Hall do Brasil, 1 995. [Ler o capítulo 1, páginas 
15-82] 

ZILL, D.G. e CULLEN, M.R. Equações Diferenciais. Volume 1. 3 a Ed. São Paulo: Makron 
Books, 2001 . [Ler os capítulos 2 e 3] 

BASSANEZI, R.C. Equações Diferenciais Ordinárias - Um Curso Introdutório. Coleção 
BC&T - Textos Didáticos - Volume 1 - Universidade Federal do ABC. Capítulo 3. Disponível 
em: 

http://gradmat.ufabc.edu.br/disciplinas/listas/iedo/notasdeaulas/equacoes-diferen- 
ciais-ordinrias-rodney.pdf [Acessado em 26/05/201 5] 

ROSA, M.S.R. Equações Diferenciais. Instituto de Matemática. Universidade Federal do 
Rio de Janeiro, 2009. Capítulo 4. Disponível em: 

http://www.dma.im.ufrj.br/~rrosa/dvifiles/apostila-ed.pdf [Acessado em 26/05/201 5] 
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Equações 
Diferenciais de 
Segunda Ordem e 
de Ordem Superior 



Até agora, nos capítulos 1 e 2, vimos o que são as equações diferenciais, suas 
várias classificações e estudamos de forma mais detalhada as equações dife- 
renciais de primeira ordem. 

Neste capítulo faremos uma discussão sobre as equações diferenciais de se- 
gunda ordem, tratando com maiores detalhes suas soluções, assim como as 
equações de ordem acima de dois, ditas equações diferenciais de ordem su- 
perior. E como as aplicações são sempre o foco do cálculo diferencial e in- 
tegral, apresentaremos também alguns exemplos clássicos como o caso do 
oscilador harmônico, com e sem amortecimento, na área de mecânica, e da 
aplicação em eletrônica para circuitos RLC (resistor - indutor - capacitor). 

^ OBJETIVOS 


■ Encontrar as soluções das equações diferenciais lineares de segunda ordem. 

■ Reconhecer e desenvolver as soluções das equações diferenciais lineares de ordem su- 
perior. 

■ Verificar algumas das principais aplicações das equações diferenciais lineares de segunda 
ordem. 


60 


CAPÍTULO 3 



3.1 Equações diferenciais de segunda 
ordem. 


Como já dissemos antes, uma equação diferencial de segunda ordem é repre- 
sentada por uma equação que possui a segunda derivada de uma função como 
o termo que apresenta a derivada de maior ordem, ou seja, a derivada de segun- 
da ordem de uma função ( y" = C L ^ ), como, por exemplo: 


y"-3y' + 4y=0 


ou 

+ x 2 — = õsenx 
dx 2 dx 


Você pode conferir com razoável facilidade que uma função (y 0 ) é solução 
de uma equação diferencial de segunda ordem substituindo y 0 na equação di- 
ferencial, como pode observar no exemplo abaixo [Adaptado de SIMMONS e 
KRANTS (2008)]: 


EXEMPLO 

Observe a equação diferencial de segunda ordem 
y" + y = - cosx 

1 

e sua solução sugerida y 0 =-— xsenx . Verifique se a solução y0 é realmente solução 
da equação diferencial. 

Resolução: Podemos derivar a função solução por duas vezes: 

, 1 1 
y = — senx — xcosx 

y 2 2 

„ 1 1 1 

y = — cosx — cosx + — xsenx 
2 2 2 

E substituir o resultado na equação diferencial, juntamente com a própria solução, ou 
seja, com 

y" + y = - cosx 


CAPÍTULO 3 ■ 61 



então 


1111 

— cosx — cosx + — xsenx — xsenx = -cosx 
2 2 2 2 


Afinal, — cosx — cosx = -cosx e —xsenx — xsenx = 0- 
2 2 2 2 

Mas você deve estar se perguntando: 

E como podemos encontrar a solução de uma equação diferencial de segunda ordem? 


Neste capítulo vamos abordar algumas das técnicas que determinam a solução dessas 
equações. Mas antes, vamos definir a forma geral de uma equação diferencial de ordem n 
como sendo: 


h 0 ( x )y( n) +hi(x)y( n +h 2 (x)y( n 2 ) +... + h(x)y = g(x) 

Onde h 0 (x)y"+h 1 (x)y’+h 2 (x)y = 0 são funções de x contínuas no intervalo I e h 0 (x)^ 0. 
Se g(x) = 0 e n = 2, então essa equação é chamada de equação diferencial lineares de se- 
gunda ordem homogênea, bastante utilizada para resolver as equações diferenciais lineares 
de segunda ordem não-homogêneas, como veremos mais à frente, e cujas soluções são 
obtidas com maior facilidade, ficando na forma: 
h 0 (x) y" + h 1 (x) y' + h 2 (x) y = 0 

Em seguida, abordaremos algumas definições importantes para se chegar às várias for- 
mas de resolução das equações diferenciais de segunda ordem, assim como as de ordem 
superior. 


3.2 Teorema de Existência da unicidade. 

Segundo ZILL e CULLEN, 2009, o Teorema que trata sobre a existência de uma 
solução única para uma equação diferencial sujeita às condições iniciais im- 
postas pelo problema, diz que: 
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“As funções h Q (x), h 1 (x), h, (x),...,h n (x) e g(x) sendo contínuas no intervalo I e 
sendo h 0 (x) y 0 para todo x neste intervalo, se x = xO como qualquer ponto neste 
intervalo, então uma solução y(x) do problema de valor inicial existe no interva- 
lo e é única”, de forma que as condições iniciais y(xO) = yO, y’(xO) = yl, y”(xO) = y2, 
...,y(n-l) =yn-l, assim como a equação diferencial de ordem n, dada por: 
h 0 (x)y (n) +h 1 (x) y 01 ’ 11 + h, (x) y (n ^ 2) + ... +h(x)y = g(x) 

devem ser satisfeitas pela função solução yCxJ, com x g dentro do intervalo I. 


XJ EXEMPLO 

Se a equação diferencial de segunda ordem dada por 
y”-9y = -18x 

tiver como solução a função y = e 3x - 2e -3x + 2x, sujeita as condições de valor inicial y(0) 
= -1 e y'(0) =11, podemos admitir que a solução y é única, pois satisfaz o teorema da exis- 
tência da unicidade da solução, uma vez que as funções h 0 (x) e h, (x) são contínuas, assim 
como g(x) e ainda h Q (x) é diferente de 0 em qualquer intervalo que contenha x = 0. 

Resolução: Veja que não é difícil verificar que a função y(x) é realmente a solução da 
equação diferencial. Ou seja: 

Dado y = e 3x - 2e~ 3x + 2x, derivando esta função por duas vezes, teremos: 
y' = 3e 3x + 6e _3x + 2 
e 

y" = 9e 3x - 1 8e~ 3x 

Substituindo o resultado da segunda derivada da função e a própria função solução na 
equação diferencial de segunda ordem y” - 9y = -1 8x, teremos: 

9e 3x - 1 8e~ 3x -9(e 3x - 2e~ 3x + 2x) = 

=9e 3x - 1 8e~ 3x - 9e 3x + 1 8e~ 3x - 1 8x = - 1 8x 

Conforme esperávamos após as substituições necessárias para a comprovação do resul- 
tado, tanto a solução como a equação diferencial são satisfeitas com as condições iniciais. 

Uma outra situação que pode ocorrer é quando for dado não apenas um ponto ( x 0 , y g) 
por onde a função solução passa, mas dois pontos diferentes ou valores de suas derivadas, 
ou seja, além do ponto (x 0 , y^, mais um outro ponto {x v y ( ) ou suas derivadas nesses pontos. 
Nesse caso teremos os chamados valores de contorno ou condições de contorno ou de 
fronteira. 
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E agora precisaremos determinar as soluções gerais das equações diferenciais lineares 
de segunda ordem ou, a partir dos valores iniciais ou das condições de contorno, determinar 
as soluções particulares dessas equações, apresentadas na próxima unidade. 

Inicialmente, apresentaremos as soluções para os casos em que a equação diferencial 
linear de segunda ordem é homogênea, uma vez que as equações não-homogêneas podem 
ser resolvidas a partir das homogêneas associadas, ou seja: 

Equação diferencial linear de segunda ordem não-homogênea: 
h 0 (x) y" + h, (x) y 1 + h 2 (x)y = g(x) 

Equação diferencial linear de segunda ordem homogênea associada: 
h 0 (x) y" + h, (x)y' + h 2 (x)y = 0 

Consideraremos também que as funções coeficientes h 0 (x), h, (x), h 2 (x) são contínuas, 
g (x) = 0 e que h 0 (x) 0. 


3.3 Soluções gerais, equações de 
coeficientes constantes e não constantes e 
redução de ordem. 

Para iniciar a discussão sobre as soluções gerais das equações diferenciais li- 
neares homogêneas devemos analisar o teorema da superposição de soluções 
que diz: 

“Considerando as k soluções da equação diferencial de segunda ordem y 1 
(x), y 2 (x), y 3 (x), ...y k (x), individualmente, a superposição, ou combinação linear 
dessas soluções, também será uma solução da equação linear”. Ou seja, 

y = c 2 y t (x) + C 2 y 2 (x) + C 3 y 3 (x) + ... + c n y n (x) 


com c ± ,c 2 , c 3 ,..., c n as constantes arbitrárias. 

Esta função é a solução geral da equação diferencial linear de segunda or- 
dem homogênea, quantas forem as soluções. Mas estamos interessados apenas 
nas soluções que são linearmente independentes, uma vez que a solução trivial 
y= 0 é sempre uma solução de uma equação diferencial linear homogênea. 
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E você se lembra o que sao funções linearmente independentes? 


É importante que você reflita sobre os conceitos de funções linearmente dependentes 
e independentes. As funções y 1 (x),y 2 (x),y 3 (x),. . .y (x) são linearmente dependentes dentro 
de um intervalo estipulado I se existir constantes c.,, c 2 , c 3 , c n , nem todas iguais a zero, 
de forma que, a combinação linear dessas funções seja igual a zero, ou seja, c, y, (x) + c 2 
y 2 (x) + c 3 y 3 (x) + ... + c n y n (x) = 0. 

Se as funções não forem linearmente dependentes, elas serão linearmente indepen- 
dentes, ou seja, as funções y, (x), y 2 (x), y 3 (x), . . ,y n (x) serão linearmente independentes se 
nenhuma delas for múltiplo de qualquer outra. 


Existe uma forma bastante prática para determinar se as funções y, (x), y 2 
(x), y 3 (x), ...y n (x) são ou não linearmente independentes. Esta forma consiste 
na montagem e resolução de um determinante de ordem n contendo na pri- 
meira linha as n funções originais e nas demais linhas suas derivadas indo 
até a derivada de ordem n, sendo que, após a primeira linha do determinante 
contendo as funções originais, a cada linha do determinante são adicionadas 
as suas derivadas sucessivas de primeira ordem para a linha dois, segunda or- 
dem para a linha três, e assim por diante. Este determinante é conhecido como 
WRONSKIANO (W), em homenagem ao matemático Hoenè Wronski (1778 - 
1853) e é dado por: 



Yi 

y 2 

y 3 --y n 

w = 

y’i 

y'a 

yá-/,. 


(n) 

(n) 

, ,( n ) ,/n) 


yl 

y 2 

Ta ■■■y [ u > 


Se o cálculo do determinante Wronskiano for igual a zero (W = 0) para to- 
dos os valores de x dentro do intervalo I, então as funções y, (x), y 2 (x), y 3 (x), ...y n 
(x) serão linearmente dependentes. Mas seW^O dentro desse intervalo, então 
as funções serão linearmente independentes, resultado que nos interessa para 
determinar a solução geral 

y=c 1 y 1 M + c 2 y 2 M + c 3 y 3 (x) + ... + c n y n (x) 


Vamos verificar com um exemplo as considerações a respeito da solução ge- 
ral construída a partir de várias funções solução linearmente independentes. 
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C7 EXEMPLO 


A equação diferencial linear de segunda ordem 


y" - 1 6y = 0 


admite cada uma das seguintes soluções individualmente: y = e 4x e y 2 = e _4x . Verifique se 
a combinação das soluções também é uma solução geral da equação diferencial, ou seja, se 
y = c, y, + c 2 y 2 , então, y = c, e 4x + c 2 e~ 4x . 

Resolução: Podemos testar individualmente as soluções y, = e 4x ey,= e~ 4x na equação 
diferencial y" - 1 6y = 0 e verificamos que ambas são, realmente soluções da equação. Mas 
será que a combinação linear das duas soluções também será uma solução da equação 
diferencial? 

Para responder isso construímos o determinante Wronskiano, dado por: 


W = -4e 4x “ 4x - 4e 4x “ 4x = -4e° - 4e° = -4 - 4 = -8 

Como W = -8 ^ 0, então as soluções são linearmente independentes e, portanto, as solu- 
ções y t = e 4x e y 2 = e~ 4x podem ser escritas como uma combinação linear, gerando a solução 
geral da equação diferencial de segunda ordem, ou seja, 


Sendo a e a 2 constantes arbitrárias. Se substituirmos a solução geral determinada na 
equação diferencial, você poderá ver que ela também satisfaz a equação. Então, não se es- 
queça de fazer a verificação, tanto das funções individuais como da solução geral. 

No caso em que temos a equação diferencial linear homogênea de segunda ordem 
h 0 (x) y" + h, (x) y' + h 2 (x) y = 0 

com coeficientes h 0 (x) = 1 , h, (x) = b e h 2 (x) = c, teremos uma equação diferencial mais 
simplificada, ou seja: 
y" + by' + cy = 0 

Onde b e c são coeficientes constantes, ou seja, independentes de x. 

Dos resultados anteriores, vemos que a solução para esta equação pode ser escrita na 
forma y = C, y, + C 2 y 2 , onde C, e C 2 (sendo C maiúsculo) são constantes arbitrárias. 

Para verificarmos uma solução padrão para a equação diferencial, utilizaremos a solução 
parcial y = e mx como uma tentativa para a solução. 


é^4x 0 -4x 

w= e , 

4e 4x 



y = c 1 e 4x + c 2 e~ 4x 
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Encontrando a primeira e a segunda derivada da solução parcial, temos: 
y' = me™ e y" = m 2 e™ 

E se substituirmos estes resultados na equação diferencial y" + by' + cy = 0, teremos: 
m 2 e™ + bme™ + ce™ = 0 

Colocarmos em evidência o termo e™, teremos a forma: 
e™ (m 2 + bm + c) = 0 

A solução dessa equação é determinada a partir da resolução da equação do segundo 
grau m 2 + bm + c = 0, uma vez que e™^ 0, sempre. Esta equação quadrática é chamada de 
equação auxiliar da equação diferencial. E talvez você já tenha notado que a equação auxiliar 
pode ser encontrada simplesmente se fizermos a associação ou substituição de y" por m 2 , 
de y 1 por m e y por 1, mantendo os coeficientes a = 1, b e c em suas posições na equação. 
Assim, devemos resolver a equação 
m 2 + bm + c = 0 

para que m, o coeficiente da solução parcial y=e™ no expoente da solução, seja deter- 
minado. E como a equação auxiliar é uma equação do segundo grau, sabemos que pode 
ocorrer três possibilidades quanto à solução da equação. Podem ocorrer duas raízes reais 
iguais e duas raízes diferentes podendo ser reais ou complexas conjugadas. Vamos, então, 
analisar cada uma das possibilidades: 

1°. caso - duas raízes reais distintas: 

Neste caso, a solução geral será da forma 

y = C 1 e m i x + C 2 e m 2 x 

Onde C, e C 2 são as constantes arbitrárias, determinadas a partir das condições iniciais 
ou de contorno e m 1 e m 2 são as constantes determinadas a partir da solução da equação 
auxiliar. 
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C7 EXEMPLO 

Seja a equação diferencial de segunda ordem 

y" - 2y' - 24y = O 

Encontre a solução da equação a partir da equação auxiliar. 

Resolução: Sabemos que a associação da equação diferencial y" - 2y' - 24y = 0 com 
a equação do segundo grau auxiliar m 2 - 2m - 24 = 0 é praticamente direta, podendo ser 
escrita dessa forma, onde y" = m 2 , - 2y' = -2 m e -24 y = - 24. 

Assim, a solução (m e m 2 ) da equação do segundo grau m 2 - 2m - 24 = 0 pode ser 
encontrada de forma bastante simples através da técnica da soma e produto das raízes da 
equação, obtendo: m = - 4 e m 2 = 6, uma vez que A > 0. 

Como a solução geral para a equação diferencial de segunda ordem reduzida à equação 
auxiliar pode ser dada por: 

y = C^ -4 * + C 2 e 6x 

Então, podemos substituir as soluções m | = -4em 5 = 6 nos expoentes da forma geral, 
obtendo: 


y = Cje 4x + C 2 e 6x 

Partindo dos valores iniciais y(0) = 3 e y'(0) = 8, como exemplo, podemos encontrar a so- 
lução particular para a equação diferencial, ou seja, C, e C 2 , ou seja, se 

y' = -4C 1 e _4x + 6C 2 e 6x 

E sua derivada como sendo 
y ’ = -4C 1 e _4x + 6C 2 e 6x 

Substituindo em y e em y' os valores iniciais, teremos: 

3 = C 1 e ( - 4( 0)) + C 2 e (6( 0)) ou 3 = C 1 +C (2) 
e 

8 = -4C 1 e(- 4 (0)) + ?6C 2 e( £3 (0)) ou 8 = -4C 1 +6C (2) 
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Resolvendo o sistema linear por uma forma bastante simples, por substituição ou soma, 


C, + C 2 — 3 

teremos \ 

[-4C 1 + 6C 2 =8 

A solução encontrada para o sistema será C t = 1 e C 2 = 2, como você pode facilmente 
verificar. 

Obtemos, assim, a solução particular para a equação diferencial y" - 2y' - 24y = 0 da 
forma 

y = e~ 4x + 2e 6x 

Caso você queira festar a solução na equação diferencial, seria um excelente exercício, 
assim como as solução individuais y 1 = e~ 4x e y 2 = 2e 6x . 

2°. caso - duas raízes reais iguais: 

Assim como no caso anterior, quando temos duas raízes iguais, teremos uma forma pa- 
drão para a solução geral, sendo dada por: 

y = C 1 e mx + C 2 xe mx 

Onde m é o valor da raiz dupla da equação auxiliar e C, e C 2 são as constantes arbitrárias. 

Então, você deve estar se perguntando: 

Como saber que a solução geral possui essa forma e não oufra qualquer? 


Isso será explicado quando falarmos sobre a técnica de redução de ordem logo à frente. 
No exemplo 5 a seguir vamos verificar a aplicação para este caso. 
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C7 EXEMPLO 

Seja a equação diferencial de segunda ordem 
y" - 8y' + 1 6y = O 

Encontre a solução da equação a partir da equação auxiliar, determinando ainda a solu- 
ção particular no caso em que y(0) = 3 e y' (0) = 17. 

Resolução: A equação auxiliar associada à equação linear de segunda ordem é dada por: 
m 2 - 8m +16 = 0 


Cujas soluções são duas raízes iguais a 4. Portanto, m = 4. 

E como a solução geral para a equação diferencial com duas raízes iguais é dada por: 

y = C, e mx + C 2 xe mx 

Então, a solução será escrita como: 

y = C, e 5x + C 2 xe 5x 

Se substituirmos os valores y = 3, quando x = 0 na solução geral ey'= 17, quando x = 
0 na primeira derivada dessa solução, então encontraremos duas equações. Resolvendo o 
sistema linear, teremos: 

r c,=3 

[5C, + C 2 = 1 7 

Encontramos, então C, = 3 e C 2 = 2, sendo a solução particular da equação diferencial 
y" - 8y' + 1 6y = 0, da forma: 

y = 3e 4x + 2xe 4x 

Esta solução também pode ser testada integralmente na equação diferencial ou ainda 
suas solução parciais, no caso y 1 = 3e 4x ou ainda y 2 = 2xe 4x , individualmente. 

Podemos ainda ter o caso em que o coeficiente h 0 (x) da equação diferencial seja dife- 
rente de 1 . Nesse caso, podemos dividir a equação toda por h 0 (x), recaindo em um dos casos 
anteriores, com raízes iguais ou distintas ou, ainda, complexas. 

Vamos desenvolver um exemplo que teremos raízes não exatas e distintas com coefi- 
ciente h 0 (x) + 1 para ver que o procedimento utilizado será o mesmo. 
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XJ EXEMPLO 

A equação diferencial dada por 
4y" + 9y' + 5y = 0 

terá uma solução particular se y = 5 quando x = 0 e y' = 13 quando x = 0. Encontre a 
solução particular dessa equação diferencial. 

Resolução: A equação auxiliar associada à equação diferencial é dada por 
4m 2 + 9m + 5 = 0 

Cuja solução é encontrada facilmente resolvendo a equação do segundo grau, obtendo 

5 

A = 1 e as raízes m = - 1 e m„ = — . 

1 2 4 


Dessa forma, a solução geral da equaçao diferencial será dada por: 

_5 

y = C 1 e _x + C 2 e 4* 

A solução particular poderá ser encontrada substituindo na solução geral as condições 
iniciais y = 5 e x = 0. Ou seja, 


5 = C,e 0 + C 2 e 4 


-#(°) 


E na sua primeira derivada da solução geral y' = 1 3 quando x = 0. Portanto, 

5 --x 

y' = “C 1 e _x - — C 2 e 4 


E com a substituição 

13 = -C 1 e-°--C 2 e _ 4 (0) 
4 


Obtendo o seguinte sistema linear: 
í C, + C 2 = 5 

-C.--C, =13 
l 1 4 2 


cuja solução será C, = 77 e C 2 = -72 . Assim a solução particular da equação particular 
será 


y = 77e x -72e ^ 


3°. caso - duas raízes complexas conjugadas: 

No caso em que a solução da equação auxiliar fornece raízes complexas, de forma que 
A < 0, teremos duas raízes distintas dadas por z, e z 2 , ou seja: 
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z, = s + ti 
Zj = s - ti 

Onde set são números reais que compõem os números complexos e i = V-T ou ainda 
i 2 = - 1 . A solução geral será da forma 
y = e sx (C, cos (tx) + C 2 sen (tx)) 


Para entender melhor a formação dessa solução geral devemos nos lembrar um pouco 
sobre séries infinitas, números complexos e a representação de e iz em termos de seno e 
cosseno. 

A série infinita que define as funções e x , cos x e sen x são da forma: 

00 y n y i y 3 

» > A . A A 

e x = > — = ] + x + — + — ... 

^ 2! 3! 


n = 0 


v 2n+1 


Z / lV 7 ^ X 

(- 1 ) — = * + ... 

v ’ (2r? + 1)l 3! 5! 


n = 0 

CO 

senx = Y J {~ 1 ) 

n-0 


n X 2n _ X 1 ^ X A 

(2n)l _ 2T + ~4d 


E se tivermos x = iz, então a série e x = e ix será escrita na forma: 


. » (iz) n (izf (iz? (iz? (iz? 

e' z = ^ — — = 1 + iz + v v v ’ + x—!— 


2! 


3! 


4! 


5! 


Obtemos: 


1 + /z / — h h / — ... 

2! 3! 4! 5! 


Se agruparmos os termos contendo i e os termos que não 0 contém, encontramos: 


z 2 z 4 | / z 3 z 5 

e' x = | 1 + ... +/ z + .. 

2! 4! j 3! 5! 


De forma que 0 primeiro termo da expressão reproduz a série da função sen z e o 
segundo termo a função cos z, podendo dessa forma escrever a função e ix em termos do 
seno e do cosseno de z, ou seja: 
e ix = cos z + isen z 


Esta expressão é conhecida como Fórmula de Euler, em homenagem ao matemático 
Leonhard Euler (1707 - 1783). 
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Se tomarmos a equação diferencial associada à equaçao auxiliar 
m 2 + bm + c = 0 

cuja solução são os números complexos z t e z 2 , então a solução geral na forma: 
y = C 1 e z < x +C 2 e z = x 

Teremos 

y = C, e (s+ti)x + C 2 e (s ' ti)x 

Cujo desenvolvimento algébrico resulta na expressão: 
y = e sx (C, e itx + C 2 e ibí ) 

Se a Fórmula de Euler pode ser escrita como 
e itx = costx + isentx 

Então, 

e~ ltx = costx - isentx 

Que associadas e substituídas na solução geral temos: 
y = e sx (C 1 cos (tx) + C 2 sen (tx)) 

Conforme descrito anteriormente. 

E agora podemos resolver uma equação diferencial que se enquadra nos termos descri- 
tos anteriormente. 


EXEMPLO 

Resolver a equação diferencial 
y" - y' + 4,25y = 0 


Sujeita às condições iniciais de y = 2 quando x = 0 e y' 


7 

3 


quando x = 0. 
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Resolução: Fazendo a associação da equação diferencial com a equação auxiliar tere- 
mos: m 2 - m + 4,25 = 0 

Onde determinamos A = -1 6, ou seja, serão duas raízes complexas e distintas, dadas por: 


m, = — + 2i 
1 2 



Em associaçao com a solução complexa, temos S = — e t = 2, resultando na solução 
geral da equação diferencial: 

2 

y = e2 (C, cos(2x) + C 2 sen(2x)) 

Submetendo a solução geral às condições iniciais teremos: 

2 = e2^ (c, cos(2(0)) + C 2 sen(2(0))) 

Ou seja, 2 = e° (C, cos (0) + C 2 sen(0)) 

Ou ainda 2= 1 (C, (1) + C 2 (0)) 

Então, C, = 2 

7 

E substituindo y’ = — quando x = 0 na derivada da solução geral, dada por: 

1 4 x 3 2 X 

y' = — e2 (C,cos(2x) + C 2 sen(2x)) + e2 (2C 2 cos(2x)-2C,sen(2x)) 

Teremos: 

— = — e 2 * * (c, cos(2(0)) + C 2 sen(2(0))) + e^ * (2C 2 cos(2(0)) - 20,560(2(0))) 

3 2 

Ou seja, 

^ = ^e°(C, cosO + C 2 senO) + e°(2C 2 0050-20,5600) 

Resultando em: 

I = ^(1)(C, (1) + C 2 (0)) + (1)(2C 2 (1) - 20, (0)) 
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Sendo 



■2 


Substituindo C, = 2 na expressão acima teremos 


7 1 

— (2) = 2C, 
3 2 V ; 2 




E a solução particular da equação diferencial será: 


y = e 2 * 2cos(2x) 


2 

3 


sen(2x) 


Cuja solução, assim como as dos outros exemplos podem ser testadas junto às equa- 
ções diferenciais. 

Para este exemplo, construímos o gráfico da solução particular exemplificando o caso de 
uma função que simula o efeito de ressonância pura, uma vez que a amplitude dos picos au- 
menta com o aumento de x, com y tendendo a infinito quando x tende a infinito e y tendendo 
a zero quando x tende menos infinito. 


Figura3.1 -Comportamentodafunçãosoluçaoparticulardaequaçãodiferencialy"-y'+4,25y=0 
simulando um oscilador harmônico forçado. 


-0,5 0,5 1 C5 2 2,5 3 3,5 1 0,5 
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-100 
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C7 EXEMPLO 

Resolver a equação diferencial 
2y" + y' + 8,1 25y = 0 

7 

Encontrando a solução particular sabendo que y = 1 e y’ = — quando x = 0 e construin- 
do o gráfico. 

Resolução: Fazendo a associação da equação diferencial com a equação auxiliar tere- 
mos: 2m 2 + m + 8,125 = 0 

Onde determinamos A = -64, ou seja, serão duas raízes complexas e distintas, dadas por: 


m. = — + 2i 
1 4 



Em associaçao com a solução complexa, temos s = — e t = 2, resultando na solução 
geral da equação diferencial: 

-1 

y = e 4 * (C, cos(2x) + C 2 sen(2x)) 


Submetendo a solução geral às condições iniciais teremos: 

1 = e 4 ^ (C, cos(2(0)) + C 2 sen(2(0))) 

Ou seja, 1 = e° (C t cos (0) + C 2 sen(0)) 

Ou ainda 1 = 1(C, (1) + C 2 (0)) 

Então, C, = 1 


E substituindo y' = — quando x = 0 na derivada da solução geral, dada por: 


-1 —X —X 

y' = — e 4 (C 1 cos(2x) + C 2 sen(2x)) + e 4 (2C 2 cos(2x)-2C 1 sen(2x)) 


Teremos: 

\ = ^ eT<0) ( C i cos(2(0))) + e^ (0) (2C 2 cos(2(0))) 


76 ■ CAPÍTULO 3 



Ou seja, 

^- = -^e° (C, cos O) + e° (2C 2 cosO) 


Sendo 


7 

4 



+ 2C 2 


Substituindo C, = 1 na expressão acima teremos 
7 1 

— + — (1) = 2C 2 
4 4 W 2 


— = 2C 2 
4 2 


C 2 = 


E a solução particular da equaçao diferencial será: 

_2 x 

y = e 4 (cos(2x) + sen(2x)) 


O gráfico que representa a solução da função tem o comportamento similar a um osci- 
lador harmônico amortecido, onde verificamos que a função tende a zero quando x tende a 
infinito. 



Figura 3.2 - Comportamento da função solução particular da equaçao diferencial 2y"+y’+8, 
1 25y = 0 simulando um oscilador harmônico amortecido. 
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Outra forma de resolver uma equaçao diferencial da forma 
y" + u(x) y' + v(x) y = 0 


é construir uma segunda solução a partir de uma solução não trivial reduzindo a equação 
diferencial de segunda ordem em uma equação diferencial de primeira ordem. Esta técnica 
é conhecida como Método de Redução de Ordem e pode ser usado para encontrar a forma 
geral da solução de uma equação diferencial de segunda ordem homogênea com coeficien- 
tes constantes ou não. 

Utilizaremos também este método para o caso em que as soluções de uma equação 
auxiliar tem raízes iguais e reais, de forma que a demonstração da forma geral será uma das 
ilustrações do método. 

O Método de Redução de Ordem parte de uma primeira solução y, conhecida e, a partir 
dessa, encontramos uma segunda solução y 2 , de forma que ambas devem ser linearmente 
independente, ou seja, a razão entre as soluções não pode ser constante dentro de um in- 
tervalo definido. Então, 


Ou ainda, como y 2 = g(x) y,, sendo y t e y 2 funções de x. 

Calculando a primeira e segunda derivadas de y 2 , e após algumas manipulações algébri- 
cas chegamos à forma geral da solução complementar y 2 como: 


De forma que a solução geral pode ser escrita a partir das duas soluções y, e y 2 na forma: 
Y = C, y, + C 2 y 2 . 

Para o caso da demonstração da solução geral no caso em que as raízes da equação 
auxiliar são iguais e reais, escrevemos a equação diferencial padrão na forma: 
y" + by' + cy = 0 

E tomamos como primeira solução a função y t = e mx , uma vez que rrq = m 2 = m. Segundo 
a solução da equação auxiliar 
m 2 + bm + c = 0 
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Uma vez que A = 0, tem como solução 


-b 



Como a = 1, podemos escrever b = -2m. 

Podemos assim, associar os termos da equação com a expressão para o cálculo da se- 
gunda solução como sendo u(x) = b = -2m. Daí, a solução poderá ser escrita como: 

e -|u(x)dx 

-J-2mdx 

y ( X ) = e mx .f— — -dx 

J ’ J ( e mx)2 

Resolvendo a integral presente no expoente do numerador temos 

e -j-2mdx = e 2mx 0 ( e mx)2 = e 2mx 

Portanto, 

- e 2mx r 

y 2 (x) = e mx • J 2mx = emx ' j^x = x '6 mx 

E a solução geral para este caso, conforme havíamos apresentado é da forma: 
y = C, e mx + C 2 xe mx 

É importante verificarmos a solução para o caso em que os coeficientes são constantes, 
sabendo previamente uma das soluções que satisfaz a equação diferencial, como no exem- 
plo 8 abaixo: 


XJ EXEMPLO 

[retirado de Zill e Cullen] Encontre a solução geral da equação diferencial homogênea: 
y" + 2y' + y = 0 


Tendo como uma das soluções a função y 1 = xeA 

Resolução: Fazendo a associação com a forma geral da equação diferencial 
y" + u(x) y 1 + v(x) y = 0 
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Tendo como u(x) = 2, teremos fazendo as substituições na expressão geral para a se- 
gunda solução, devendo ambas serem linearmente independente, então; 


y 2 ( x ) = xe- x -J ( 


-|2dx 


-dx 


(xe x ) 2 

E com as substituições, teremos: 


y 2 (x) = xe- x -J. 


3 -J 2dx 

xe~ x ) 2 


dx 


Como a integral presente no expoente do numerador será: 


J2dx = 2x 
Então, 

r e _2x 

y 2 (x) = xe - x .J- 5 -^dx 

Simplificando o integrando, teremos: 
y 2 ( x ) = xe _x • j-^-dx 


Cujo valor da integral é dado por: 

y 2 (x) = xe- x -^J 



-1 

— , então: 
x 


Ou seja, 

y 2 (x) = -e x 

Será a outra solução para a equação diferencial. A combinação das duas solução nos 
resulta em: 

y = C, xe x - C 2 e x 

Sendo essa a solução geral da equação diferencial y" + 2y' + y = 0. 

Seria um excelente exercício verificar se as soluções individuais, assim como a composta 
com as soluções y e y , são realmente soluções da equação diferencial. 
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Mas também é importante verificarmos a solução para o caso em que os coeficientes não 
são constantes. 0 método também é válido, desde que se saiba uma das soluções da equa- 
ção diferencial e, a partir dela, encontramos uma segunda solução, formando uma solução 
que comporta as duas individuais, sendo essas linearmente independentes. Então, vamos 
desenvolver um exemplo em que os coeficientes não são constantes, utilizando assim o 
Método de Redução de Ordem para encontrar a solução geral de uma equação diferencial 
de segunda ordem homogênea. 


Ache a solução geral da equaçao diferencial, sabendo que uma das soluçoes da equação é 
é Y, = x 4 : 


Para podermos aplicar o Método de Redução de Ordem devemos ter um formato padrão 
da equação diferencial, assim como saber uma das soluções que satisfaça a equação. Para 
transformar na forma utilizada no método é necessário que a equação diferencial seja divi- 
dida por x 2 , obtendo: 


EXEMPLO 


x 2 y" - 7xy' + 1 6y = 0 



Ou seja, 



Portanto, u(x) 


-7 


ey ) = x 4 . Substituindo os dois resultados na forma geral: 


x 



-ju(x)dx 


Teremos: 
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Se lembrarmos da propriedade de logaritmos que diz a log » = x 7 , então: 

y 9 (x) = x 4 • f í-dx = x 4 ■ f-dx = x 4 -Inx 
J x a J x 


Portanto, a segunda solução será: y 2 (x) = x 4 ■ Inx. 

A solução geral que é uma combinação linear das soluções y, e y 2 será: 
y = C t x 4 + C 2 x 4 ■ Inx 

Até agora resolvemos equações diferenciais de segunda ordem homogêneas. Na pró- 
xima seção discutiremos a solução para as equações diferenciais de segunda ordem não 
homogêneas. 

A seguir vamos ver algumas aplicações para equações diferenciais homogêneas e não 
homogêneas, como o caso do oscilador harmônico simples não amortecido ou de um circuito 
elétrico RLC. 

Todo sistema físico em equilíbrio, quando submetido a forças externas que tendem a 
desequilibrá-lo, tende a voltar ao estado original graças a forças restauradoras existentes no 
sistema por natureza. 

Este é o caso de um sistema massa-mola em que a massa m é retirada de seu estado de 
repouso e equilíbrio a partir de uma força externa. No estado de equilíbrio ou repouso, a mola 
não exerce qualquer força sobre a massa, mas se a massa for submetida a uma força externa 
que a retira de seu estado de equilíbrio, então a mola aplicará uma força sobre a massa a fim 
de tentar restaurar a situação original. A equação de movimento para esse sistema físico é 
representado por uma equação diferencial de segunda ordem, dada por: 
m ■ x" = -kx 

Podendo ainda ser escrita na forma: 
m ■ x" + kx = 0 


Onde m é a massa do corpo que foi retirada de seu estado de equilíbrio, x(t) é o deslo- 
camento da massa a partir do repouso e k é a constante elástica da mola ou constante de 
Hooke, sendo k>0 

Se dividirmos a equação diferencial por m e substituirmos a 2 =— , então teremos uma 

m 

equação diferencial bastante conhecida da forma: 

„ k 

x h — x = 0 
m 
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Ou ainda 
x" + a 2 x = 0 

Cuja solução é dada por: 
x(t) = C, sen(at) + C 2 cos(at) 

Se submetermos o sistema às condições iniciais de x(0) = x 0 , que representa a massa 
ser deslocada para a direita e solta a partir do repouso, ou seja, sua velocidade inicial é zero, 
ou seja, v(0) = 0. Como a velocidade é a primeira derivada da função posição em relação ao 
tempo x'(t) = v(t) então, podemos encontrar a solução particular da equação diferencial que 
representa o deslocamento da massa no tempo, que permanece nesse estado de movimento 
desde que não haja forças externar para cessar o movimento, como a força de atrito, uma 
situação ideal. 

Então, substituindo na solução geral as condições iniciais, temos: 
x(t) = C, sen(at) + C 2 cos(at) 

Oue substituída, nos dá: 

x 0 = C, sen(a(0)) + C 2 cos(a(0)) 

Ou ainda, 

x 0 =C, (0) + C 2 (1) 

Portanto, C 2 = x Q 

Para a primeira derivada da solução, representada pela velocidade, teremos: 
x' (0) = 0 = aC, cos(a(0)) - aC 2 sen(a(0)) 

Dando, 

0 = aC, cos(a(0)) 

Ou 

0 = aC, (1) 


Portanto, C, = 0. 

Então, a solução particular será da forma: 
x(t) = x Q cos(at) 


CAPÍTULO 3 ■ 83 


Caso existam forças externar, o sistema é chamado de amortecido, o que torna a situa- 
ção mais próxima do que realmente ocorre. 

No caso de termos forças externas presentes no sistema, então podemos modificar a 
equação diferencial fazendo a soma de todas as forças atuanfes, como por exemplo: 

mx" + cx' + kx = F 0 cos t 

Também tomada como exemplo de uma equação diferencial de segunda ordem, agora 
não homogênea. 

Uma outra situação também bastante conhecida é o caso de um circuito elétrico com- 
posto por um capacitor (C), um indutor (L) e um Resistor, também chamado de circuito RLC 
em que atua uma força eletromotriz E = E 0 cos t, cuja carga Q no capacitor é calculada pela 
equação diferencial: 

LQ" + RQ '+ — Q = E n cos cot 
C 0 

Onde L é a indutância o indutor, R é a resistência do resistor e C é a capacitância do ca- 
pacitor. Essa equação pode ser associada facilmente à equação do oscilador harmônico for- 
çado dada anteriormente, pois são sistemas similares devido à atuação das forças externas. 


Ey CONEXÃO 

As aplicações para equações diferenciais de segunda ordem são bastanfe utilizadas em 
sistemas físicos, computacionais e de engenharia. 

E importante que você leia um pouco mais sobre as várias aplicações. Abaixo alguns links em 
que você encontra algumas das aplicações mais clássicas. 
http://www.ime.uerj.br/~calculo/LivrolV/edoseg.pdf 

http://efisica.if.usp.br/mecanica/universitario/movimento/ocilador_harm_amortecido/ 
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3.4 Equações diferenciais lineares de 
segunda ordem não homogêneas. Soluções 
particulares. 


A equação diferencial de segunda ordem não homogênea tem a forma: 
h 0 (x)y" + tq (x)y' + h 2 (x)y = g(x) 

Onde g(x) ? 0. A solução geral é dada pelo conjunto das soluções que satisfa- 
zem a equação diferencial homogênea, chamada de solução geral da equação 
complementar y , dada por: 
h 0 (x)y" + hj (x)y' + h 2 (x)y = 0 

Além da combinação linear com a solução particular y . Ou seja, 

y=y c +y P 

Como sabemos resolver as equações diferenciais de segunda ordem homo- 
gênea, de acordo com nosso estudo feito até agora, então fica fácil determinar 
a solução para uma equação diferencial de segunda ordem não homogênea, 
principalmente se for conhecida a solução particular dessa equação. Vamos en- 
tender melhor o método resolvendo mais um exemplo. 


EXEMPLO 


Seja a equação diferencial e sua solução particular, dadas por: 
y" - 6y' + 5y = -9e 2 * com y p = 3e 2x . 


Primeiramente, precisamos resolver a equação diferencial homogênea dada por: 
y" - 6/ + 5y = 0 


Encontrando a equaçao auxiliar associada a ela teremos a seguinte equaçao do segundo 
grau: 

m 2 - 6m + 5 = 0 
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Cuja solução desenvolvida através de soma e produto das raízes, será m = 1 
e m 2 = 5. E como as raízes são reais e distintas, sabemos que a solução geral será dada 

por: 

y c = C 1 e m,x + C 2 e m 2 x 

Substituindo as raízes m_1 e m_2, teremos: 
y c = C, e 1x + C 2 e 5x 
Então, a solução geral será da forma: 
y = C, e ,x + C 2 e 5x + 3e 2x 

Existe ainda a possibilidade de termos uma sequência semelhante de equações diferen- 
ciais de segunda ordem, cujo termo g(x) se diferenciam entre elas, assim como suas soluções 
particulares y . Neste caso, a solução da equação diferencial composta por todas as funções 
g(x) é uma sobreposição de todas as soluções também, além da solução complementar. 


^ CONEXÃO 

E no link abaixo você poderá observar uma séries de imagens reproduzindo os osciladores 
harmônicos, como nos pêndulos de pintura: 
https://www.youtube.com/watch?v=3rigdgiHmnQ 
ou uma aula sobre osciladores harmônicos 
https://www.youtube.com/watch?v=y4sB6oGZel_s 


3.5 Equações diferenciais de ordem superior. 

Em outras situações que a equação diferencial é de ordem superior a dois 
podemos escrever a forma geral da equação diferencial como: 

h 0 (x)y tn, + h 1 (x) y (n ' 1) + h 2 (x) y 1 ”' 25 + ... + h(x)y = g(x) 

Que escrita na forma homogênea ou na forma da equação complementar: 

h Q (x)y (n) + fq (x) y (n ' 1) + h 2 (x) y 1 "' 2 ^ ... + h(x) y = 0 
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Podendo a ela associar uma equação polinomial auxiliar de grau n, dada 
por: 

amM+a ,m( n_1 )+a 0 m* n_2 ) + ... + a.m + a n =0 

n n— 1 n -2 1 U 

De forma que os coeficientes a_n são constantes dentro dos números reais. 

A forma da solução geral da equação diferencial homogênea de ordem n de- 
penderá do número de soluções de mesmo valor da equação auxiliar tiver. Ou 
seja: 

• se todas as raízes da equação auxiliar forem diferentes entre si, o formato 
da solução geral será: 

y = C ie m,x + C 2 e m 5 x + C 3 e m 3 x + . . . + C n e m " x 


EXEMPLO 

[retirado de Zill e Cullen] Dada a equação diferencial de terceira ordem homogênea, encontre 
sua solução geral 
y'" - 4y" - 5/ = 0 

Resolução: Encontrando a equação auxiliar associada à equação diferencial de ordem 
3, temos: 

m 3 - 4m 2 - 5m = 0 

Resolvendo a equação colocando inicialmente m em evidência e resolvendo a equação 
do segundo grau, teremos: 
m(m 2 - 4m - 5) = 0 

Cujas soluções serão: m = 0, m 2 = -1 e m 3 = 5, as duas calculadas por soma e produto 
das raízes da equação do segundo grau. Dessa forma, a solução geral será: 
y = C, e 0x + C 2 e ( ' lx) + C 3 e 5x 

ou ainda 

y = C, + C 2 e w + C 3 e 5x 

Para o caso em que duas das raízes da equação auxiliar são iguais ou complexas, então 
a combinação de possibilidades é bastante grande, dependendo da ordem da equação dife- 
rencial. Dessa forma a equação diferencial terá a forma: 


CAPÍTULO 3 ■ 87 



y = C 1 e m < x + C 2 xe m * x + C 3 x 2 e m 3 x + . . . + C n x n “ 1 e m n x 

Uma vez que as soluções individuais devem ser línearmente independentes, podendo se 
repetir o termo apara raízes diferentes. 


XJ EXEMPLO 

[retirado de Zíll e Cullen] Resolva a equação diferencial: 
y"" + 2y" + y = 0 

Encontrando a solução geral da equação: 

A equação auxiliar associada será: 
m 4 + 2m 2 +1=0 

Podemos escrever a equação na forma de um produto notável para que a solução seja 
mais facilmente encontrada, ou seja: 
m 4 + 2m 2 + 1 = (m 2 + 1 ) 2 = 0 


As raízes da equação auxiliar serão: rrq = m 2 = i e m 3 = m 4 = -i. Portanto, a forma geral 
da solução será: 

y = C, e“ + C 2 xe“ + C 3 e ( ix) + C 4 xe (4x) 

E se aplicarmos a Fórmula de Euler, podemos reescrever a solução em termos de senos 
e cossenos, ou seja: 

y = C t cosx + C 2 x ■ cosx + C 3 senx + C 4 x ■ senx 


ET ATIVIDADES 

09. Resolva a equação diferencial particular abaixo: 

y" - 4y' + 1 3y = 0, com y = 6 quando x = 0 e y' = 10 quando x = 0. 

1 0. Encontre a solução geral da seguinte equação diferencial: 
y" - 2y' + 2y = 0 
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1 1 . Dada a equação diferencia xy" + y' = 0 e sabendo que y t = Inx é uma solução, encontre 
uma segunda solução da equação diferencial escrevendo uma solução geral. 

1 2. Seja a equação diferencial de segunda ordem 

y"+y'+^y = ° 

Encontre a solução da equação a partir da equação auxiliar, determinando ainda a solu- 
ção particular no caso em que y(0) = 2 e y'(0) = 1 . 


E7 REFLEXÃO 

Todos os métodos desenvolvidos até agora são métodos para apenas alguns tipos de equa- 
ções diferenciais. E será que existe alguma outra forma de resolver essas equações diferen- 
ciais? Ou será que os métodos aqui estudados seriam capazes de resolver todos os tipos de 
equações diferenciais? Nos próximos capítulos estaremos discutindo um pouco mais sobre 
isso, mostrando as técnicas de Transformada de Laplace e Fourier. 


^ CONCEITO 

6. Equação diferencial de segunda ordem: é uma equação que apresenta uma segun- 
da derivida de uma função como a derivada de maior ordem na equação. 

7. Equação auxiliar: é uma equação algébrica associada a uma equação diferencial de 
segunda ordem ou de ordem superior em que uma variável equivale à derivada na equação, 
como m 2 = y". 

8. Oscilador harmônico: é um sistema físico cujo corpo em movimento alcança a mesma 
posição repetidas vezes, de forma harmônica, podendo ser Simples, Amortecido ou Forçado. 

9. Circuito RLC: é um sistema composto por três componentes eletrônicos, um resistor 
(R), um indutor (L) e um capacitor (C). Esses componentes, quando ligados entre si formam 
um circuito elétrico conhecido como RLC. 

1 0. Wronskiano: é uma função desenvolvida a partir do cálculo de um determinante com- 
posto por funções originais e suas derivadas para definir se um conjunto de funções são 
linearmente dependentes ou independentes. 
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EM LEITURA 

É importante que você leia um pouco mais sobre os Osciladores harmônicos. O link em des- 
taque abaixo mostra algumas definições para que você possa estar mais informado sobre o 
assunto. Além dessa aplicação outras aplicações e outros métodos de resolução de equa- 
ções diferenciais de segunda ordem que não foram abordados aqui podem ser encontrados 
nos textos indicados abaixo. Além também de existirem muitas aplicações das equações 
diferenciais, principalmente na área de Física e Engenharia. Os modelos que envolvendo 
equações diferenciais foram capazes de desenvolver diversos avanços nessas áreas, assim 
como explicar diversos efeitos como o famoso caso da Ponte Tacoma Narrows, falando um 
pouco mais sobre o efeito de Ressonância. Abaixo apresentamos também alguns links e 
textos para que você possa se ínterar um pouco mais sobre essas aplicações, príncipalmente 
na área de tecnologia. 

(ZILL, D.G. e CULLEN, M.R. Matemática Avançada para Engenharia. Equações diferen- 
ciais elementares e transformada de Laplace. Volume Único. 3 a Ed. Porto Alegre: Bookman, 
2009. Capítulo 3, unidade 3.4 (Coeficientes indeterminados), pag 139 a 146.) 

(ZILL, D.G. e CULLEN, M.R. Equações Diferenciais. Volume 1. 3 a Ed. São Paulo: Makron 
Books, 2001. Capítulo 5, unidade 5.3 (Movimento forçado), pag. 248 a 256.) 
http://www.seara.ufc.br/tintim/fisica/ressonancia/ressonancia6.htm 
http://www.sbmac.org.br/cmac-se201 1/trabalhos/PDF/332.pdf 
http://www.peb.ufrj.br/cursos/eel420/EEL420_Modulo7.pdf 
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A Transformada de 
Laplace 



Neste quarto capítulo vamos discutir os conceitos e propriedades da 
transformada de Laplace. Como todo procedimento de transformação, a 
transformada de Laplace tem especial importância na simplificação de algu- 
mas equações diferenciais. Para a aplicação do processo de transformação, 
faremos uso de vários conceitos, como, por exemplo, o de integral imprópria 
e de função que sejam de ordem exponencial. 

H7 OBJETIVOS 


■ Entender os conceitos de integrais impróprias. 

■ Desenvolver a transformada de Laplace. 

■ Entender as propriedades e as condições de existência da transformada de Laplace. 

■ Conhecer a transformada das derivadas e a inversa da transformada. 

■ Entender a aplicação das transformadas e o uso da tabela das transformadas. 
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4.1 Conceitos e noções básicas. 


O que é uma transformada? 

Para que serve uma transformada? 

Bem, para dar início a nossa discussão da Transformada de Laplace, é im- 
portante que antes tenhamos respostas para as duas perguntas iniciais. 

Uma transformada nada mais é do que um processo de transformação de 
um problema inicial, que pode ser mais complexo do que gostaríamos, em um 
problema equivalente mas com menor grau de complexidade. Assim, ao re- 
solvermos o problema transformado, supostamente mais fácil, obtemos mais 
rapidamente e com menor esforço a resposta do problema, mas não só do pro- 
blema transformado, mas também do problema original. Essa é a validade e a 
utilidade de uma transformada, facilitar o cálculo de um dado problema. 

Temos muitos exemplos na Matemática de problemas que podem ser trans- 
formados em outros mais fáceis de serem resolvidos. 

Por exemplo, no cálculo integral você aprendeu pelo menos duas técnicas 
de integração que podem ser interpretadas como rotinas de transformações. 
Você se lembra de alguma delas? 

Sim, a técnica de substituição de variáveis é uma forma de transformar uma 
integral não muito simples de ser resolvida em outra mais simples, como no 
caso do exemplo 1. 


EXEMPLO 


Uso de mudança de variável na resolução de integrais. 
Seja a integral definida I dada por: 


z. 

1= Jsen(x)-cos(x)dx 


Fazendo a transformação u=sen(x), temos que du = cos(x)dx e os limites passam para: 
= sení — 1 = — - eu = sení —1 = 1. Assim a inteqral transformada é: 

UJ 2 sup 1,2 J a 


1 

I = j udu 

•Jã 

2 
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Ainda no campo da resolução de integrais, você também aprendeu a técnica da integra- 
ção por partes, que possui a mesma finalidade, facilitar a resolução da integral, como no caso 
do exemplo 2. 


ty EXEMPLO 

Uso da técnica da integração por partes na resolução de integrais. 
Seja a integral definida I dada por: 


1 = 


4 


Jx-e 3x dx 

i 


Fazendo a transformaçao u=x e dv=e3x, du=dx e v = ~ e3x ■ Assim a integral transfor- 
mada é: 


1 = 


x-e' 


,3x 



Outro caso que podemos listar que representa uma transformação é quando tratamos de 
integrais numéricas. Além das técnicas de integração já conhecidas por você, temos também 
as técnicas de resolução de integrais definidas por meio de aproximações numéricas, como 
por exemplo, a regra dos trapézios ou a regra dos 3/8 de Simpson, que são regras fechadas 
de Newton-Cotes. Apesar de ser um método de aproximação, também podemos entender 
como um tipo de transformação, já que transformamos a função original a ser integrada em 
outra função, no caso um polinómio, para simplificar o cálculo. Veja o exemplo 3. 


E y EXEMPLO 

Uso da regra dos 3/8 de Simpson para aproximar uma integral. Seja a integral I igual a dada 
no exemplo 2: 

4 

l = Jx-e 3x dx 

Utilizando uma aproximação da função através de polinómios de grau 3, considerando 
7 pontos igualmente espaçados de h=0,5, temos uma transformação da integral da forma: 
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1 5 

1 = y • ( V 0 + 3 Yi + 3y 2 + 2 Y 3 + 3 Y4 + 3 Y 5 + Ve ) 


Com: 

yJ=x_í-e A (3x_í ) 

x_0= 1 

x_í=x_(i-1)+0,5, para i=1,2, .... 6 

Assim, a integral original é transformada em uma expansão baseada na interpolação por 
um polinómio de grau 3. 

Outro caso que podemos citar e que serve muito bem para entender o que ocorrerá com 
os problemas após o uso da Transformada de Laplace, são as transformações geradas pelo 
uso do logaritmo em uma função exponencial. Lembre-se que o logaritmo tem a capacidade 
de transformar uma operação de exponenciação em um produto, uma multiplicação em uma 
soma e uma divisão em uma diferença, ou seja, é capaz de simplificar os cálculos associados 
a um problema. Observe o exemplo 4. 


XJ EXEMPLO 

Uso de logaritmos para simplificar uma equação exponencial. Seja a equação 

20 = 1 5 ■ (4,5) x 1 

Aplicando o logaritmo natural na equação, temos: 

In20 = ln[1 5 ■ (4,5) x '] In20 = In 15 + In (4,5) x 1 In20 = In 15 + (x - 1) In4,5 


E assim a equação inicial é transformada em uma equação linear que pode ser facilmen- 
te resolvida. 

É claro que quando falamos em transformadas ou transformação não nos limitamos ape- 
nas a situações de cálculo. 

Temos também transformações associadas a imagens e sons, como, por exemplo, as 
modificações feitas pelo aplicativo photoshop em uma foto, ou então, as alterações na grava- 
ção de uma música feita por um mixer. 

Esses casos envolvem a utilização de transformadas, seja nas configurações da imagem 
ou na configuração de sinais e ondas. 
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Por exemplo, observe a transformação da onda representada pela linha pontilhada em 
outra onda mais suave representada pela linha contínua (figura 5.1), ou a mudança em uma 
imagem através de uma transformação por um filtro gaussiano (figura 5.2). 



Figura 5.1 - Onda original com linha pontilhada transformada em outra onda em linha con- 
tínua. 

É possível notar pelos exemplos discutidos anteriormente que é natural a utilização de 
transformações para se obter resultados melhores ou facilitações em termos de cálculos. A 
transformada de Laplace é um desses casos, uma transformada que possui como objetivo 
representar de maneira particular um determinado problema a fim de facilifar seu cálculo. 



Imagem normal Imagem transformada 

Figura 5.2 - Imagem original com transformação por um filtro gaussiano. 


Contudo, temos mais uma discussão a fazer anfes de começar a trabalhar com a trans- 
formada de Laplace. Você se lembra das integrais impróprias? Façamos uma recordação 
rápida. 
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Toda vez que calculamos uma integral definida, mesmo que não seja do seu interesse 
avaliar tal resultado, mas sempre podemos associar essa integral a uma medida de área. 
Assim, sempre é possível analisar uma integral definida do ponto de vista geométrico a uma 
área sob a curva, levando-se em conta uma região delimitada pelos limites da integral. Veja 
uma representação desse fenômeno no exemplo 5. 


EXEMPLO 


B 

l = jf(x)dx 

A 

Geometricamente, a integral pode ser representada pela figura 3. 


A B 

Figura 5.3 - Área sob a curva f(x) no intervalo [A,B], que equivale ao valor de /. 

Note que essa representação é de fácil interpretação pois a integral é definida em um 
intervalo fechado. Agora, sempre teremos intervalos fechados nos limites de uma integral? O 
que pode ocorrer com esses limites de forma a termos um intervalo aberto? 

Basta que o limite inferior seja igual a infinito negativo ou o superior seja igual a infinito 
positivo. 

Considerando essas possibilidades, temos, na prática, três possibilidades de escrever 
uma integral definida em termos de limites abertos, mesmo que seja semiaberto, vejamos: 

Caso 1 : integral partindo do limite inferior A e indo até infinito, [A,+¥) ouxwA. 

Assim, temos: 

00 

1 = jf(x)dx 

A 
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Caso 2: integral partindo do limite superior B, partindo de infinito negativo, (- oo, B] ou 
x £ B. 

Assim, temos: 


B 


J f(x)dx 


Caso 3: integral partindo de infinito negativo e indo até infinito, (- ¥, ¥) ou x ÍR. 
Assim, temos: 

CO 

I = j f(x)dx 

—co 


É natural admitir que nesses três casos não seja razoável pensarmos em uma medida 
de área, contudo, isso não impede que as integrais impróprias possam ser calculadas e até 
mesmo associadas a uma medida aproximada do que poderia ser admitida como a área. 
Vejamos um exemplo de como trabalhar com as integrais impróprias. 


XJ EXEMPLO 

Suponhamos uma integral definida do tipo exposto no caso 1, mas ao invés de denotarmos 
o limite superior por infinito, vamos denotar por M. Assim, seja a seguinte integral definida: 


'-I 


i 


1 


dx = 




1 

M 


+ 1 


Como não sabemos o valor de M, não temos como determinar a área sob a função, mas 
podemos intuir que quanto maior for o valor de M, mais próximo estaremos da verdadeira 
área. Assim, se tomarmos o limite do valor da I com M tendendo para infinito, podemos en- 
contrar uma boa aproximação para o resultado. Assim: 
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Área Aproximada = lim 



— - + 1 = 0 + 1 = 1 


00 


Note que esse raciocínio nos leva a concluir que para trabalhar com integrais impróprias 
não basta dominar as técnicas de integração, mas também temos que dominar as operações 
com limites no infinito. 

Do exposto, temos três importantes propriedades derivadas dos três casos apresentados 
anteríormente e que serão úteis na operação das transformadas de Laplace. 

Propriedade 1 : Seja f(x) uma função integrável em um intervalo [A,+¥). Então: 


+00 


B 



A 


A 


Propriedade 2: Seja f(x) uma função integrável em um intervalo (- ¥, B]. Então: 


B 


B 


1= f f(x)dx = lim íf(x)dx 

J A— >-co J 


—00 


A 


Propriedade 3: Seja f(x) uma função integrável em um intervalo (- ¥, +¥). Então: 



o 


o 


B 


EXEMPLO 


Calcular a integral imprópria I dada por: 



—00 
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Aplicando a propriedade 3, temos: 


/= fe X dx = íe X dx+ fe X dx = lim fe X dx+ lim fe X dx = 

J J J A—±— qqJ R— >oo ^ 


A->-a 


lim — - 


0 

+ 

lim f 

A 

B — >co y 


A->-a 


1 1 


lim — -H + lim — — h — — = 


B — >oo 


1 1 


/A— >-ool e' 


lim — -1 + lim 1 — -\ = 


-1 + 1 


= (e” -l) + (1-0) = (oo-1) + 1 = oo 

Convencidos de que o uso de transformadas é bem comum na matemática e com essa 
breve revisão sobre as integrais impróprias, podemos partir para a transformada de Laplace. 
Vamos lá. 


4.2 Propriedades da Transformada de 
Laplace. 

Em 1812, Pierri Simon de Laplace (1749-1827) publicou uma obra intitulada 
Teoria Analítica, e nesta apresentou a Transformada de Laplace de uma função 
f(t), que passará a ser denotada por £{f(t)}. 

Da mesma forma que os cenários discutidos no tópico 4.1, a transformada 
de Laplace aparece como uma maneira de simplificar alguns tipos de proble- 
mas matemáticos, em especial algumas equações diferenciais. 

Neste contexto, o que se espera com a aplicação da transformada é gerar, a 
partir de uma equação diferencial ordinária, uma equação algébrica mais fácil 
de ser trabalhada e resolvida. Com a solução dessa equação algébrica, fazemos 
o caminho de volta, através da inversa da transformada de Laplace para obter a 
solução da equação original. 

Esquematicamente podemos entender o processo da seguinte forma: 
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Figura 5.4 - Esquema de utilização da transformada de Laplace. 

Bem, como então é a transformada de Laplace? Vejamos. 

Definição: Dada uma função integrável f(t), tal que f:[ 0 ,°°) — > R, a transfor- 
mada de Laplace, , é dada por: 

co 

£{/(?)} = F(s) = Je“ s, f(í)dí 
o 

Para todo s > 0 de maneira que a integral tenha convergência e com s = o +ico 
uma variável do plano complexo 

É possível notar que a transformada de Laplace envolve uma das formas dis- 
cutidas anteriormente de integral imprópria. Assim, fica evidente que as pro- 
priedades associadas aos limites desse tipo de integral serão essenciais para se 
trabalhar com as transformadas. 

Por motivos de conveniência e concordância com a literatura específica das transfor- 
madas de Laplace, utilizaremos a variável t como argumento da função original e a variável 
s para a função transformada. Ainda, a função original será sempre representada por letras 
minúsculas, f(t), g(t), h(t), enquanto que a função transformada será representada por mai- 
úsculas, F(s), G(s), H(s). 
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Ey CONEXÃO 

Você pode notar que a variável s associada a transformada de Laplace é uma variável comple- 
xa. Não vamos abordar detalhes sobre os complexos neste capítulo, mas você pode aprofun- 
dar seus conhecimentos sobre as variáveis complexas lendo as referências listadas abaixo: 
www.icmc.usp.br/~szani/complexa.pdf 

www.impa.br/opencms/pt/biblioteca/cbm/06CBM/6_CBM_67_09.pdf 


Vejamos um exemplo inicial de utilização da transformada de Laplace. 


C7 EXEMPLO 

Encontrar a transformada de Laplace para a função f(t)=1. 
Aplicando a definição, temos: 


£{f(f)} = F(s) = je st f(t)dt = je sí -1 dt 


Aplicando agora as propriedades de integrais impróprias, temos: 

B 


CÍf(t) } = F(s)= lim fe S M dt = lim Í-- — 1 

1 ’ W B-> oo^ s J 


0 


Í_£ÜL 

L— ] 

-1 

f g-sco 

L— ] 

-I 

r_2 + l1 

1 s 

s J 

i 

L s 

s J 

rl 

1 S sj 


s 


Assim F(s) = - , para s > 0. 
s 

Da mesma forma anterior, podemos encontrar a transformada de Laplace para uma 
gama enorme de funções, sejam polinomiais, trigonométricas ou exponenciais, ou até mes- 
mo uma mistura das mesmas. 

Faremos isso em um tópico posterior. Neste momento, a fim de facilitar a obtenção das 
transformadas de outras funções, é importante estabelecer algumas propriedades para a 
transformada de Laplace. 
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Propriedade da Linearidade: 

Se a transformada de Laplace de f(t) é F(s) , para s > ai e a transformada de Laplace de 
g(t) é G(s) , para s > a2, então considerando duas constantes a e b, temos: 

£{af(t) + Pg(t)} = a£{f(t)} + p£{g(t)}=ccF(s) + |3G(s) 
para s > maxja^aj. 

Propriedade do Deslocamento: 

Se a transformada de Laplace de f(t) é F(s), para s > a e considerando uma constante a, 
então a transformada da função 

g(t) = e at f(t) 
será igual a 

G(s) = F(s-a) 
para s > a + a. 

Propriedade da Convolução: 

Se a transformada de Laplace de f(t) é F(s) , para s > a 1 e a transformada de Laplace de 
g(t) é G(s) , para s > a 2 , então considerando duas constantes a e b, temos: 

£{f(t).g(t)} = £{f(t)}.£{g(t)} = F(s).G(s) 


para s > maxja^aj. 

Condições de Existência: 

Para que exista a transformada de Laplace F(s) de f(t) é necessário que: 

I) a função f(t) seja contínua em cada intervalo entre dois pontos quaisquer de desconti- 
nuidade, caso existam; 

II) a função f(t) seja de ordem exponencial, ou seja, deve existir uma constante a, com a 
pertencente aos reais de modo que exista lim |/(f)|e' aí 

t-> 00 1 1 

III) o domínio de F(s) de f(t) será s > a. 

Função de ordem exponencial é toda função tal que existam constantes c, M > 0 e 
T > 0 de maneira que |f(t)| < M.e 0 *, para todo t > T. 
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Para melhor entendimento do critério (i), observe a figura 5.5. Note que a função possui 
diversos pontos de continuidade, mas entre cada dois pontos a função é contínua. Isso basta 
para possibilitar o cálculo da transformada de Laplace. 



Figura 5.5 - Função f(t) descontínua nos pontos t e t 2 . 

A verificação da validade das condições de existência da transformada de Laplace parte 
do conceito de função de ordem exponencial e da propriedade da soma de integrais defini- 
das. Vejamos: 

CO T co 

£{/(?)} = Je“ s, f(f)c/f = je -s 7(f) c/f + ^e~ st f(t)dt = /, +/ 2 

0 0 T 


Nota-se claramente que a primeira integral está sendo calculada em um intervalo con- 
tínuo então o seu valor existirá sem problemas. Contudo, na segunda integral não necessa- 
riamente, mas supondo que a função que está sendo integrada seja de ordem exponencial, 
valendo a relação |f(t)| < M.e ct , para todo t > T, podemos afirmar que: 


|/ 2 | < J|e _s '/(f)|c/f < M - je~ sf -e c, dt = M ■ Je c '>‘dt = M- 


-(s-c)r 


Assim, paras > c, como existe a convergência de / 2 , a transformada também apresentada 
convergência, existindo, portanto, a transformada de Laplace para a função. Mas é importan- 
te atentar para o fato de que o resultado só é válido se a função for de ordem exponencial 
e se s > c. 
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Assim, para funções cuja velocidade de crescimento seja superior a função exponencial, 
não sendo, portanto de ordem exponencial, não existirá a transformada de Laplace. 

Vejamos algumas representações gráficas para melhor compreensão. Nas figuras de 6 
a 1 0 estão representadas funções exponenciais, através da linha contínua, e outros tipos de 
funções nas linhas pontilhadas. 



Figura 5.6 - Função f(t)=e°' 4l +2 é de ordem exponencial. 



Figura 5.7 - Função f(t)=2,5t é de ordem exponencial. 



Figura 5.8 - Função f(t)=e ‘+1 é de ordem exponencial. 
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Figura 5.9 - Função f(t)=1 ,5cos(t) é de ordem exponencial. 



Figura 5.1 0 - Função f(t) = e J ' 2 não é de ordem exponencial. 

Aplicar a definição e as propriedades da transformada de Laplace é possível obter a 
transformadas para um conjunto bem amplo de funções, o que será útil nos casos de aplica- 
ções, visto que executando as transformadas a priori, podemos gerar uma tabela de resulta- 
dos que podem ser consultados a posteriori. 

Preparado? Vamos lá. 

Vamos começar com a função constante f(t)=k. Anteriormente já vimos o desenvolvi- 
mento da transformada de Laplace para a função f(t)= 1 , mas como será a transformada para 
uma função genérica f(t)=k? 

A definição será da forma: 

00 B 

C{k\ = [e~ st kdt = k lim [e~ st dt 

1 S->ooJ 

0 0 

Uma vez que a função f(t)=k é constante em t, então podemos resolver a integral por 

substituição, fazendo a substituição u=-st e, com isso, dt = —du ■ Então, 

s 
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Aplicando os limites de integração, lembrando que e~ st — > 0 quando t — » a> e que e~ st = 1 
quando t = 0, quando s > 0, temos: 


C{k) = — lim 

S B— »co 


f 

B > 


e~ st 


= A Hm ( e -^_ e -s(0)\ 



s B-> oo' ' 

V 

OJ 



'T< 


Ou seja, a transformada de Laplace para uma função do tipo f(t) = k será: 

Válida para qualquer k real e com s>0. 

Vejamos agora o caso para a função linear f(t) = t. 

Pela definição, temos: 

00 B 

£{k} = Je“ sf fc/í = lim Je -s, fc/f 
o s_> °°o 


Podemos dispensar a notação de limite para agilizar um pouco mais os cálculos, mas 
sempre lembrando da ideia de limite que deve ser levada em consideração para o cálculo da 
integral. 

Para a resolução da integral devemos usar a técnica de integração por partes, onde: 

^udv = uv - Jvduj 


Fazendo u = t, du = dt, dv = e~ st dt e 


i/ = 


2 e — sf i temos: 

s 


C{t} = je st tdt 
o 
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Ou ainda, 


00 

— - e st — 1- e st 

f e~ st tdt = 

J 

0 

s s 2 j 


00 

0 


Substituindo os limites de integração e lembrando que: 
lim Be~ sB = 0 

B — >oo 

e 

lim e~ sB = 0 

B — >oo 

Temos, 


£{t} = je s, tdt 
o 


-e~ s (°) + -U- s ( 0 ) = 0 + -1 = 4 

s s 2 s 2 s 2 


Portanto, a transformada de Laplace, F(s) (em maiúsculo e em função de s) para a fun- 
ção f(t) = t (em minúsculo e em função de t), será: 

F(s) = £{?}= ^ 

s 


Se desenvolvermos a transformada de Laplace para uma função quadrática, cúbica e 
assim por diante, ou seja, generalizando a função polinomial, com n = 1 , 2, 3, 4,. . teremos 
a forma geral: 


F(s) = C{t n ) 


n\ 

s n+ 1 


Então, para uma função do tipo: f(t) = at 2 , a transformada de Laplace será: 
C{af}=a.C{f} = a.^ = | 


3) Conforme visto na seção anterior, para o caso em que temos somas de funções, 
quantas forem elas, aplicamos a transformada de Laplace em cada função (f(t) + g(t)), mul- 
tiplicadas pelos fatores a e b. A transformada de Laplace, que neste caso é chamada de 
transformada linear, pode ser escrita na forma: 
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C{af(t) + bg(t)} = je sl (af(t) + bg(t))dt 
0 


Ou ainda, como: 

OO 00 

F(s) + G(s) = aC[f(t)} + b£{c/(f)} = a.je~ s, f(t)dt + ò.jV sí gf(f )df 

o o 


Vamos resolver um exemplo para verificar como é mais fácil o desenvolvimento no caso 
em que algumas funções já têm sua transformada conhecida. 


U EXEMPLO 

Seja a função quadrática dada por: f(t) = 3x 2 + 2x - 4, encontre a transformada de 
Laplace para este caso. 

Resolução: A forma geral da transformada de Laplace para a função dada será: 

00 

F(s) = C{f(t)} = \e- st (3x 2 +2x-4)dt 
0 

Ou ainda na forma: 

00 00 00 

F(s) = £{/(?)} = je' sí (3A" 2 )c/f + Je~ s, (2x)dí- Je“ s, (4)c/f 
ooo 

Mas, como já sabemos a forma geral para a função polinomial, então podemos desenvol- 
ver a transformada de Laplace com os resultados já obtidos, ou seja: 

F(s) = £{3* 2 + 2* - 4} = £{3* 2 } + £{2*} + £{-4} 
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então teremos: 


Lembrando que F(s) = £{f n } = -^-j- , 

£{3^} + £{M + £!-4}=3| + 2.^4.1 


Resultando em: 

F(s) = £{ 3 x 2 +2x-4} = 4 + 4-- 
' ’ S d s 2 s 


como sendo a transformada de Laplace para a função f(t) = 3x 2 + 2x - 4. 

Vejamos agora o caso em que a função original é da forma exponencial, ou f(t) = e a1 
Podemos desenvolver a transformada de Laplace fazendo: 

OO 00 00 

F(s) = £{f(f)j = £{e af j = Je -sí (e a, )c/f = Je _sí+a 'c/f = Je"( s-a ) f c/f 


Resolvendo a integral por substituição, fazendo u = -(s - a) t, então: 

00 00 .. U 00 

F(s) = £{e a '} = Íe^ s_a ^c/í = í e u du = f e u du 

o o ( s ~ a ) -( s - a ) J o 


F( s) . £ { e »J._L_( e -l.-.>) ._L_(o-e-(-l») 

Então, a transformada de Laplace da função exponencial f(t) = e at será: 

F(s) = £{e a, }= — 

' ’ s-a 


Considerando agora a função original como uma função trigonométrica do tipo 
f(t) = sen(at). Podemos desenvolver a transformada de Laplace fazendo: 


F(s) = £|f(f)| = £{sen(af)} = je S, (sen(af))df 
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Resolvendo a integral por partes, fazendo u = sen at, du = acos(at)dt e ~ e s,Q,f , então 


\e~ st dt = — e~ st , e substituindo em J udv = uv - J vdu, temos: 


je Sí (sen(af))c/f = 


-e S, senat 


+ — [e sí cos (at)dt 

Q J 


Sabendo que: 


lim e sB sen(aB) = 0 
6— >00 


Então, 

-e~ st senat 

s 


00 

= 0 

0 


e 


je st (sen(at))dt = — ji 


cos(af)df 


Por partes, mais uma vez, teremos u = cos (at), du = -a sen(at) dt e dv = e" st dt, então 

y — f Q~c^dt — e _sí * então. 

s 


Je st (sen(atj)dt = — 
o s 



00 

00 

— e s 'cos(aí) 

o 

+ J — e _sí (sen(af))c/f 

O 

0 

O 

O 


Ou ainda na forma: 

00 

Je' sf (sen(at))dt = -^-e~ st cos(at) 


--Q-fe sl (sen(at))dt 


Como 

lim e~ sB cos(aB) = 0 

B-> oo v ’ 
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Então, 


je Sí (sen(af))df = je Sí (sen(af))c/f 

o s s o 

Agrupando as integrais, teremos: 


Je Sí (sen(af))df + ^- je Sí (sen(af))df = 


Ou na forma: 


^2+1 Je st (sen(at))dt = -^ 


a 2 a 2 + s 2 
72 + = ^2 


Isolando a integral e fazendo s s teremos: 

00 

^ st inseri [atfjdt - 


as * 


? (s 2 + a 2 ) 


Ou ainda, 


F(s) = £{sen(af)| = je S, (sen(af))c/f = - 


ty EXEMPLO 

Qual a transformada de Laplace para a função f(t) = sen (3t)? 

Resolução: Poderíamos resolver a integral e seguir todo o processo, mas é mais prático 
utilizarmos o resultado genérico encontrado para a função f(t) = sen(at), então, 

F{s) = C{sen{3t)}=^-^ 

0 desenvolvimento para a função f(t) = cos(at) é similar à f(t) = sen(at) e sendo assim, 
o resultado encontrado será: 
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s 


F(s) = C{cos(at)} = Je st (cos(at)^dt 
o 


U EXEMPLO 

Qual a transformada de Laplace para a função f(t) = cos(4t)? 

Resolução: Utilizando o resultado já encontrado, temos que a transformada de Laplace 
para a função f(t) = cos(4t), será: 

F(s) = C{cos{4t)} = -^ 

Dessa forma vemos que, sabendo os resultados das transformadas de Laplace de fun- 
ções bases, podemos encontrar de forma bastante fácil o resultado para outras funções 
similares desenvolvidos genericamente. 

Podemos ainda verificar a forma da transformada para funções contínuas por partes. 
Neste caso, devemos analisar separadamente parte da função, somando os resultados indi- 
viduais. Então, vejamos. 


U EXEMPLO 

Para a função com duas sentenças dada por: 



paraO < t < 5 
para t > 5 


calcule a sua transformada de Laplace. 

Resolução: A partir da definição: 

CO T 00 

£{f(t)} = Je“ s 'f(f)c/í = Je -s, f(f)c/f + Je“ sí f(í)c/í = /, +/ 2 

0 0 T 

Podemos construir a transformada em termos das integrais da forma: 
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£{/(?)} =Je s, f(f)df = Je s '/(f)df + Je sl f(t)dt = /, + / 2 


o 


o 


T 


Como a primeira das duas integrais é igual a zero, então será necessário resolver apenas 
a segunda derivada, resultando em: 


É interessante observar que para cada tipo de função existe uma regra definida para 
a forma da transformada de Laplace. Dando continuidade à discussão das propriedades e 
regras associadas às transformadas, com os resultados que serão apresentados nos dois 
tópicos a seguir, é possível gerar uma diversidade bem significativa de formas para as trans- 
formadas, que serão apresentadas posteriormente em forma de tabela, de maneira a facilitar 
a aplicação, da mesma maneira que ocorre com as tabelas de derivadas e integrais. 

Por falar em derivadas e integrais, vejamos o que ocorre quando operamos as transfor- 
madas nesses aspectos. 


Vamos iniciar com as derivadas. De início, o que ocorre quando aplicamos o 
conceito da transformada de Laplace na derivada de uma função? 

Você pode estar se perguntando: qual a importância disso? Bem, como já 
mencionado anteriormente, um dos objetivos da transformada de Laplace é 
simplificar o manuseio de equações diferenciais transformando-as em equa- 
ções algébricas. Para isso, as transformadas devem ser aplicadas nas equações 
diferenciais e você vai se lembrar que além das variáveis e das funções, as equa- 
ções diferenciais também possuem derivadas. Assim, ao calcular a transforma- 
da de uma EDO, estamos transformando também as derivadas. 



00 


Portanto, 



3e 5s 


4.3 Derivada da transformada e 
transformada da derivada. 
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Vejamos o que ocorre com a transformada de uma derivada de primeira or- 
dem. Suponhamos uma função f(t) integrável, tal que f:[0,«q ->R,ea transfor- 
mada de Laplace L{f(t)}. Considerando a derivada de primeira ordem de f(t) 
como sendo f’(t) então, pela definição da transformada de Laplace, temos: 

oo 

C{f’(t)} = je- St f’(t)dt 
o 


Resolvendo a integral por parte, temos u=e~(-st) e dv=T ' (t). 

Assim, du = - s • e' st e v = f(t). Substituindo na regra da integral por partes, 
temos: 


00 


00 


£{f’(t)}^e st f(i) 


+ sje st f(t)dt 
0 0 


O primeiro termo pode ser resolvido utilizando limite, da forma: 
B 

lim [e _s, f(í)] = lim [e~ sS f(B)-e“ st V(0)] = [0-/(0)] = -f(0) 

o e ^°° 


O segundo termos, se você observar bem, é exatamente a transformada da 
função f(t). Assim, o resultado fica reduzido a: 

£{f’(t)} = -f(0) + s-C{f(t)} 

ou da forma: 

C{f’(t)}=s-F(s)-f(0) 


Fazendo agora para a derivada de segunda ordem, temos: 

oo 

jC.{f"(t)} = \e- st f"(t)dt 
o 

Resolvendo a integral por parte, temos u = e' st e dv = f ' (t). 

Assim, du = -s • e st e v = f (t). Substituindo na regra da integral por partes, 
temos: 
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00 


+ sje s ‘f'(t)dt 
0 0 


O primeiro termo pode ser resolvido utilizando limite, da forma: 
B 


lim í" e sí /’(f)l 
e->» L J 


0 


= Jim [e- sB f’{B) - e- s0 f (0)] = [O - f’( 0)] = -f’( 0) 


Novamente no segundo termo temos uma transformada, mas neste caso da 
função f (t). Assim, o resultado fica reduzido a: 


C{f"(t)}=-f’(0) + s-C{f’(t)} 


Substituindo a transformada da derivada de primeira ordem, temos final- 
mente que: 

£{/"(f)} = -r(0) +s .[ s .F(s)-/(0)] 

E portanto: 

£{/"(f)} = -/'(0) + s 2 ■ F(s) - s ■ /(0) 


Ou então, da forma: 

£{/"(f)}=s 2 -F(s)- S -/(0)-f(0) 

Fazendo agora para a derivada de terceira ordem, temos: 


£{F"(f)} = Je- sí F”(?)df 

o 

Resolvendo a integral por parte, temos u = e st e dv = f " (t). 

Assim, du = - s • e st e v = f ’ (t). Substituindo na regra da integral por partes, 
temos: 
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£{f”(0} = e- s, r(f) 


00 


00 


+ sje~ st f"(t)dt 
0 0 


O primeiro termo pode ser resolvido utilizando limite, da forma: 
B 


lim [~e sí F’(f)l 


= jim [V sS r(s) - e - s0 r(o)] = [o - r(o)] = -r ( o) 


Novamente no segundo termo temos uma transformada, mas neste caso da 
função f”(t). Assim, o resultado fica reduzido a: 

£{r(t)} = -f"(o) +s-c{r(t)} 


Substituindo a transformada da derivada de segunda ordem, temos final- 
mente que: 

£{f ”(í )} = -r(0) + s • [s 2 ■ F(s) - s ■ f (0) - f( 0)] 


E portanto: 

£{f”’(f)} = + s 3 • F(s) - s 2 • f (0) - s • /'(O) 


Ou então, da forma: 

£{r(t)}=s 3 -F(s)-s 2 -f{0)-s-f'^-f"(0) 


Comparando os três casos, podemos generalizar o que ocorre com a transfor- 
mada de Laplace de uma derivada de ordem n, gerando a seguinte propriedade: 

Transformadas de Derivadas: 

Sejam as funções f(t), f‘(t), f“(t), f“‘(t), ..., f(n-l)(t) contínuas em [0, ¥) e todas 
de ordem exponencial e a função f(n)(t) contínua por partes em [0, ¥), então a 
transformada de qualquer derivada será igual a: 

£ {f( n) (f )} = s n F(s) - s^V (0) - si"” 2 V'(°í fí"- 1 ) (0) 
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Essa propriedade será muito útil em situações de transformada de equa- 
ções diferenciais, como veremos mais adiante. 

Da mesma maneira que temos uma transformada de uma derivada, pode- 
mos ter também a derivada de uma transformada. Conforme colocado por ZILL 
e CULLEN (2001), tomando a transformada de uma função f(t) e calculando a 
sua derivada em relação a s, temos: 

- S - = — Je" sí /(f)dí = j— [e -sí f(f )]df = ~^e~ st tí(t)dt = 
ds ds Q 0 ds 0 

Logo, temos que: 

Apesar do objetivo do resultado ser mostrar o que ocorre com a derivada da 
transformada, o mesmo acaba sendo mais uma forma de generalizar os pro- 
cessos de transformada, no caso quando temos a transformada do produto da 
função f(t) por um monómio tn. 

De maneira generalizada, temos outra importante propriedade: 

Transformada de funções do tipo t n f(t): 

d'"'(0} £;F(‘) 


E CONEXÃO 


Além das derivadas, podemos trabalhar também com as integrais das transformadas, ou 
as transformadas de integrais. Para maiores detalhes sobre essas propriedades, estudo o 
artigo de VIANA, R.L Transformadas Integrais. DF, UFPA, Curitiba (PA), 201 3, disponível em 
fisica.ufpr.br/viana/metodos/transformadas.pdf 
Acessado em 29/06/2015. 
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4.4 A transformada Inversa. 


Vimos anteriormente, na figura 4 em especial, que a transformada de Laplace 
pode ser utilizada para simplificar a operação matemática de algumas equa- 
ções diferenciais mais simples e os resultados anteriores em especial quanto às 
derivadas serão fundamentais para isso. 

Contudo, temos ainda outro conceito a ser discutido. Na figura 4, é possí- 
vel notar que a transformada será capaz de reduzir uma EDO em uma equação 
algébrica, facilitando a sua resolução, mas a partir dessa resolução, temos que 
transformar de novo, fazer a volta da transformação, a inversa da transforma- 
ção para associar a solução da equação algébrica à solução da EDO. 

Como já fizemos anteriormente o cálculo da transformada para algumas 
funções típicas é claro que a inversa da transformada, L' 1 , deve retornar para a 
função original. Vamos verificar cinco casos principais. 

Resultado 1. 


Para f(t)=k temos que } = C{k\ = - 
1 ’ s 


Então, a inversa da transformada será igual a: = k 



Resultado 2. 


Para f(t)=ktn temos que 

Então, a inversa da transformada será igual a: £{^(0} = £{kt n ) = 


Resultado 3. 


Para f(t)=eat temos que £{f(0} = £{e aí } 


s-a 


Então, a inversa da transformada será igual a: Cr ] = kt n 



Resultado 4. 
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Então, a inversa da transformada será igual a: \ = sen(kt) 

\s 2 +k 2 


Resultado 5. 


Para f(t)=cos(kt) temos que £{f(f)} = £{cos(/cf)} = ■ 


; 2 +k 2 


Então, a inversa da transformada será igual a: C 1 


s 2 + k 2 


m 


Vamos exemplificar o uso das transformadas inversas para melhor 
compreensão. 


C7 EXEMPLO 

1 2 

Calcular a transformada inversa da função F(s) = 

çO 


Podemos reescrever a função da seguinte forma: 

F(s) = — ~ 1 2- — - 12-—L - 1 2 ■ — -_L - 1E . _Ü!_ 

vv s e s 6 S 5+1 51 S 5+1 5! S 5+1 


A função está exatamente como mostrado no resultado 1 . Assim, a transformada inversa 


será igual a: 




1 9 

— -f 5 =0,1f 5 
120 


ty EXEMPLO 


Calcular a transformada inversa da função 


3 

s 2 +1 6 


Podemos reescrever a função da seguinte forma: 


F(s) 


_J_ = 3 .^ = 3 .^ = 3 .4.^_ 

s 2 + 1 6 s 2 + 1 6 s 2 + 4 2 4 s 2 + 4 2 


3 4 

4 s 2 + 4 2 
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A função está exatamente como mostrado no resultado 4. Assim, a transformada inversa 
será igual a: 


Cr 1 



3 

— ■sen 

4 


(4t) = 0,75sen(4t) 


XJ EXEMPLO 

Calcular a transformada inversa da função F(s) = — 0 + 

w s 2 +9 

Podemos reescrever a função da seguinte forma: 

\_-4 + 5s_ -4 5s _ -4 3 _ s 

lS, “^ 2 79““?73 2+ ?73 2 “T'^ 2 73 2+ '^T 2 


A função está exatamente como mostrado nos resultados 4 e 5. Assim, a transformada 
inversa será igual a: 


-4 + 5s 
s 2 +9 


• sen(3f) + 5 • cos(3f) = 5 • cos(3í) - 0,75 ■ sen(3f) 


XJ EXEMPLO 

Encontre a função original f(t), dada a transformada de Laplace na forma: 
rí \ 7s ~ 16 

^ s 2 - 5s + 6 

Resolução: Devemos inicialmente fatorar o denominar, dado por: 
s 2 -5s + 6 =(s-2)(s-3) 


e em seguida desenvolver as frações parciais, ou seja: 
7s-16 _ 7s-1 6 _ A B 

s 2 -5s + 6 (s-2)(s-3) s-2 s-3 
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Desenvolvendo a fração parcial para encontrar os valores de A e B, temos: 


7s-16 _ A(s-3) B(s-2) _As-3A + Bs-2B 

(s-2)(s-3) ~~ s-2 + s-3 ” (s-2)(s-3) 


E de forma comparativa, agrupando os termos em s e constantes, temos: 

7s-16 _{A + B)s + (-3A-2B) 

(s-2)(s-3) _ (s-2)(s-3) 

Ou seja, (A + B) = 7 e (-3A - 2B) = -16. Então, resolvendo o sistema linear para deter- 
minar os valores de A e B, temos: 


í A+B= 7 
[-3A-2B = -16 


Isolando o termo A em uma das equações e substituindo na outra equação, encontramos 
o valor de B = 5. A partir do valor de B podemos encontrar o valor de A = 2. 

Então, a expressão na forma de fração parcial será escrita como: 

7s-1 6 _ 2 5 

s 2 - 5s + 6 s-2 s-3 


Lembrando que 



s-a 


Podemos escrever a função f(t) fazendo a analogia entre as expressões, considerando 
ainda as constantes multiplicativas, ou seja: 


o 

— = 2e 2í 
s-2 


5 


s-3 


= 5e 3f 


Então, 


f(t) = 2e lt +5e 3í 
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4.5 Tabela da transformada de Laplace. 


Considerando todos os resultados obtidos anteriormente, podemos organizá- 
-los de uma forma adequada para consultas posteriores em forma de tabela. 

Além dos resultados já obtidos e levando-se em conta resultados colocados 
em ZILL e CULLEN (2001), apresentada uma lista de funções e suas respectivas 
funções transformadas de Laplace na tabela 1. 


■ 

FUNÇÃO F(T) 

TRANSFORMADA F(S)=L{F(T)} 

1 

k 

k 

s 

2 

t 

1 

s 2 

3 

t n 

n\ 

s n+1 

4 

+-1/2 

£ 

5 

p/2 

W 

2s 3/2 

6 

sen(kt) 

k 

s 2 + k 2 

7 

cos (kt) 

k 

s 2 -k 2 

8 

senh(kt) 

k 

s 2 -k 2 

9 

cos (kt) 

k 

s 2 -k 2 

10 

sen 2 (kt) 

2 k 2 

s(s 2 +4k 2 ) 

1 1 

cos 2 (kt) 

s 2 + 2k 2 
s(s 2 +4k 2 ) 

12 

senh 2 (kt) 

2/c 2 

s(s 2 - 4/c 2 ) 
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13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 


cosh 2 (kt) 
tsen(kt) 
tcos (kt) 

sen(kt) + ktcos(kt) 
sen(kt) - ktcos(kt) 
tsenh(kt) 
t cosh(kt) 

1 - cos (kt) 
kt - sen (kt) 
cos (at) - cos(bt) 
sen(kt) senh(kt) 
sen(kt) cosh(kt) 
cos(kt) senh(kt) 

cos(kt) cosh(kt) 

2(l-cos(af)) 

t 

2(l-cosh(af)) 

t 

sen(af) 

t 


s 2 — 2/c 2 
s(s 2 - 4/c 2 ) 

2/cs 

(s 2 +/t 2 ) 2 

s 2 -/c 2 

(s 2 +/c 2 ) 2 

2/cs 2 

(s 2 +/c 2 ) 2 

2/c 3 

(s 2 +/t 2 ) 2 

2/cs 

(s 2 -/c 2 ) 2 

s 2 +/c 2 
(s 2 -/c 2 ) 2 

/c 2 

s(s 2 +/c 2 ) 
k 3 

s 2 (s 2 +/c 2 ) 

s(b 2 - a 2 ) 
s 2 +a 2 )(s 2 + ò 2 ) 

2s/c 2 
s 4 + 4/c 4 

/c(s 2 +2/c 2 ) 
s 4 +4/c 4 

/c(s 2 - 2/c 2 ) 
s 4 + 4/c 4 


s 4 + 4/c 4 


a 

arctg- 


CAPÍTULO 4 



30 

sen(at)cos(bt) 

t 

1 a+b 1 a-b 

— arctq + — arctq 

2 s 2 s 

31 

e kt 

1 

s-k 

32 

te kt 

{s-kf 

33 

t" e w 

n\ 

{■ s-k) n+ 1 

34 

e 4 f 

Q-kyfs 

35 

k 4 

e 4 í 

2/7^ 

Q-ky/s 

36 

e at sen(kt) 

k 

(s-a ) 2 +/c 2 

37 

e at cos(kt) 

s-a 

(s-a ) 2 +/c 2 

38 

e al senh(kt) 

k 

(s-a ) 2 - k 2 

39 

e at cosh(kt) 

s-a 

(s-a ) 2 -k 2 

40 

g af _ e bt 

1 

a-b 

(s-a)(s-ò) 

41 

ae at - be bt 

s 

a-b 

(s-a)(s-b) 

42 

e bt _ e af 

t 

. s-a 

In 

s-b 

43 

e at f(t) 

F (s - a) 

44 

f'(t) 

sF(s) - f(0) 

45 

f" (t) 

s 2 F(s) - sf(0) - f(0) 

46 

f n > (t) 

s" F(s)-s (n - 1 ) f(0)-...f<"- , )(0) 

47 

f , t , í°, paraO <t <a 
[l, para t > a 

e -as 

S 


Tabela 5.1 - Transformadas de Laplace 
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De posse da tabela e de todas as propriedades discutidas anteriormente, 
podemos iniciar as aplicações das transformadas de Laplace na resolução de 
equações diferenciais. 


4.6 Aplicações. 

Agora nesta seção, vamos estudar como as equações diferenciais são resolvidas 
com a ajuda das transformadas de Laplace, uma vez que já vimos todos os con- 
ceitos necessários para isso. 

1. Determine a solução da equação diferencial de segunda ordem, dados 
os valores iniciais, y(0)=0 e y'(0)=2, sendo y"+y'-2y=4t. 

Iniciamos aplicando a transformada de Laplace em ambos os lados da 
equação diferencial: 

Iniciamos aplicando a transformada de Laplace em ambos os lados da 
equação diferencial: 

C(y") + C(y')-2C(y) = AC{t) 

Com os valores da tabela da seção anterior temos de (45), (44), (2) e lembran- 
do que L(y) = F(s), que é transformada de interesse, podemos escrever: 

s 2 F(s)-sí(0)-F(0) + sF(s)-/(0)-2 F(s) = 4^- 

s 2 

Substituindo os valores iniciais como f(0) = 0 e f (0) = 2, temos: 

s 2 F(s)-s(0)-2 + sF(s)-0-2 F(s) = 4-^ 


Agrupando os termos em F(s), temos: 


(s 2 + s-2)F(s) = + 2 : 


4 + 2s 2 


Isolando F(s) na expressão, fatorando o denominador e separando a expres- 
são em frações parciais, temos: 


F(s) 


4 + 2s 2 
s 2 (s 2 +s-2) 


4 + 2s 2 

s 2 (s + 2)(s-1) 


ABC D 

1 — — H 1 

s s 2 s + 2 s - 1 
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Tirando o mínimo múltiplo comum e agrupando os termos em s3, s2, s e 
constantes, temos: 


Fazendo a associação entre os numeradores da expressão, o sistema linear 
a ser resolvido será: 

A+C+D = 0 
A + B -C + 2D = 2 


De onde você pode ver facilmente, partindo da última equação, determinan- 
do B = -2 e fazendo as devidas substituições chegamos aos valores de A = -1, C = 
-1 eD = 2. 

Com isso, a transformada de Laplace passa a ser vista como: 


Associando os resultados da tabela de transformadas de Laplace por 1, 2 e 
31, temos: 

/(f) = -1-2f -e -2 ' + 2e' 

É um excelente exercício fazer a verificação do resultado na equação dife- 
rencial. Tente e verifique a veracidade do resultado para ter maior confiança no 
método. 

2. A equação diferencial em que y é uma função de t (f(t)) 
y"+ 4y = t 

Está sujeita às condições iniciais y(0) = 2 e y'(0) = 3. 



s 3 (A + C + D) + s 2 (A + B-C + 2D) + s(B-2A)-2B 
s 2 (s + 2)(s-1) 


B-2A = 0 
-26 = 4 
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Aplicando a transformada de Laplace na equação diferencial, temos como 
resultado: 

C(y") + 4C(y) = C(t) 

Da tabela de transformadas de Laplace, 2 e 45, e ainda lembrando que 
L(y) = F(s), então: 

s 2 F(s)-sf(0)-f'(0) + 4F(s) = -1 

Ou ainda, 

s 2 F(s) + 4F(s)-2s-3 = \ 

Ou seja, 

(s 2 + 4)F(s) = 2s + 3 + -Jr- 
s 2 

Então, 

F( >_ 3 2s 1_ 

(s 2 +4) (s 2 +4) s 2 (s 2 +4) 


O último termo deve ser dividido em duas frações na forma de frações par- 
ciais, obtendo: 

1 _ 41 B _s 1 (A + B) + AA 

s 2 (s 2 +4) s 2 s 2 + 4 s 2 (s 2 + 4) 


Resolvendo o sistema linear em A e B, temos A = — e B = — ■ 

4 4 


Com isso, F(s) torna-se: 


^00 = 3 ^ 


1 


2s 


\( 1 


(s 2 +4)J (s 2 +4) 4 s 2 y ^ 4 y v (s 2 + 4) 


+ - - - 


1 


128 ■ CAPÍTULO 4 


Fazendo a operação entre o primeiro e o último termo da expressão e mani- 
pulando algebricamente o resultado para poder ser comparado à transformada 
de Laplace da tabela no. 6 e ainda utilizando as de no. 7 e 2, temos: 


Sendo essa a solução da equação diferencial que você também pode con- 
ferir substituindo o resultado, além de sua segunda derivada, na equação 
diferencial. 

3. [Adaptado de Simmons e Krantz] Determine a solução da equação dife- 
rencia abaixo sujeita às seguintes condições iniciais y(0) = 0 e y'(0) = 3. 

y "+ 2y 1 ' + 5y = 3 e^sent 

Mais uma vez devemos aplicar a transformada de Laplace na equação 
diferencial 

£(y") + 2£(y’) + 5£(y) = 3£(e~'senf) 

A partir dos resultados da tabela de transformadas de Laplace 36, 44 e 45, 
podemos escrever: 

í 1 1 

s 2 F(s) - sf( 0) - f'( 0) + 2(sF(s) -f( 0)) + 5F(s) = 3 ^ — 



Ou seja, 



^(s + 1) + 1, 


Considerando as condições iniciais, temos: 


r 


\ 


s 2 F(s) - s( 0) - 3 + 2(sF(s) - 0) + 5F(s) = 3 2 

U s + 1 ) +1 


y 


Isolando F(s), a expressão é escrita como: 
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E ainda, 


F(s) = 3 


(s 2 + 2s + 2)(s 2 + 2s + 5) 


s 2 + 2s + 5 


O primeiro termo da expressão deve ser escrito como frações parciais na 
forma: 


1 


1 


(s 2 +2s + 2)(s 2 + 2s + 5) (s 2 +2s + 2) (s 2 +2s + 5) 


Adaptada novamente à expressão de 

Fí \ ] 1 3 

F(s) = 


FO)F(s) 


s 2 + 2s + 2 s 2 + 2s + 5 s 2 + 2s + 5 


temos: 


Ou na forma: 

F(:] 1 | 2 

s 2 + 2s + 2 s 2 + 2s + 5 


Escrevendo os dois denominadores na forma de produtos notáveis, temos: 


F( s ) 


1 2 

( S -(-1)f+1 + (s-(-1)) 2 + 2 2 


Recorrendo mais uma vez a transformada 36 da tabela, a função f(t) passa 
a ser: 

/(f) = e^sent + e~'ser)(2f) 

Como mais uma verificação, faça a substituição da função e suas derivadas 
na equação diferencial para conferir se a função realmente é solução da mesma. 
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ET ATIVIDADES 

1 3. Para os exercícios abaixo obtenha a transformada de Laplace a partir da tabela: 

a) f(t) = 3t 5 

b) f(t) = 3t 2 - e- yt + 8 

c) f(t) = 5 (e* - e~') 2 

d) f(t) = 2 - 2cos (3t) 

e) f(t) = 2t 3 + 4e~ 3t - 7 sen 2 2t 

1 4. Para os exercícios abaixo obtenha a transformada inversa: 
a) 

b) ^'{(f 


15. Encontre a solução das equações diferenciais abaixo a partir das transformadas de 
Laplace: 

a) y' + 6y = e 31 , para y(0) = 1 

b) y" - 6y' + 9y = t 2 e 3t , para y(0) = 2 e y' (0) = 6 



REFLEXÃO 

Discutimos neste capítulo as transformadas de Laplace, que, assim como outros métodos de 
transformação, visa a melhorar o cenário de um problema ou modelo matemático a fim de 
torná-los mais simples. 

Em especial, as transformadas de Laplace servem para simplificar o cálculo de alguns 
tipos de equações diferenciais. 

No entanto, as transformadas não sevem apenas para simplificar cálculos, mas servem 
também para modificar a forma de representação de um modelo matemático. 

Por exemplo, em estudos de variação dos preços de um ativo na bolsa de valores, o que 
se observa é uma variação muito acentuada com um gráfico formado por ciclos com muitos 
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picos. Para se interpretar o que ocorre com o preço é comum adotar funções especiais, 
denominadas de filtros, que possam suavizar o comportamento, como por exemplo, médias 
móveis, modelos autoregressivos ou splines. Todos esses filtros são funções de transforma- 
ções, são transformadas, com o objetivo de produzir novas funções. 

Em situações de sinais elétricos ou de ondas, como as ondas quadradas associadas aos 
bits do computador, existe uma classe especial de transformada, denominada transformada 
de Fourier, que estudaremos no próximo capítulo. 


^ CONCEITO 

Alguns termos citados neste capítulo são de especial importância e sendo assim reforçamos 
seus conceitos neste glossário. 

1 1. Integral imprópria: integrais definidas de maneira que pelo menos um dos limites de 
integração seja igual ao infinito e para a resolução dependem do cálculo de limites; 

1 2. Transformada de Laplace: a transformada L{f(t)} de uma função f(t) integrável em um 
intervalo [0, °°) é definida por uma integral imprópria de uma função de ordem exponencial 
do tipo e stf (t) ; 

1 3. Função de ordem exponencial: é uma função f(t) cuja taxa de crescimento é menor do 
que uma função exponencial e é tal que |f(t)| < Me ct ; 


Ef LEITURA 

Para se aprofundar nos conceitos de transformadas de Laplace e suas aplicações, recomen- 
damos a leitura das obras listadas abaixo: 

ANDRADE, D. Transformada de Laplace. DM, Universidade Estadual de Maringá. 

www.dma.uem.br/kit/arquivos/arquivos_pdf/transforlaplace.pdf 

Acessado em 29/06/201 5. 

LIMA, P.C. Equações Diferenciais. DM, Universidade Federal de Minas Gerais. 

www.mat.ufmg.br/~lima/apostilas/apostila_eda.pdf 

Acessado em 29/06/201 5. 
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Séries de Fourier 



Neste quinto e último capítulo, vamos discutir os conceitos e propriedades das 
séries infinitas. Sobre as séries infinitas, ampliaremos a discussão conside- 
rando uma série muito especial, a série de Fourier, que possui muitos aplica- 
ções práticas associadas a transformação de funções em uma série de senos e 
cossenos. 


Hf OBJETIVOS 


■ Entender os conceitos de séries infinitas e suas propriedades. 

■ Entender os conceitos de séries periódicas e trigonométricas. 

■ Entender as conceituações de séries pares e ímpares. 

■ Desenvolver as séries de Fourier. 


136 


CAPÍTULO 1 



O ano era 1822, o país era a França, quando foi publicado um trabalho sobre a 
teoria analítica do calor, o autor, Jean Baptiste Joseph Fourier (1768-1830). Este 
trabalho, poucos anos antes de sua morte, marca definitivamente a presença 
de Fourier na Matemática e na Física e a sua transformada, a série de Fourier. 

@ 

I 


Jean Baptiste Joseph Fourier (1 768-1 830) 

A elaboração de seu trabalho, como apontado por MATOS (2002), foi estu- 
dar a propagação do calor em um corpo sólido, um trabalho focado na área de 
Termodinâmica, mas o que Fourier não esperava é que os resultados que obte- 
ve tivessem um impacto da magnitude que houve. 

Durante mais de 150 anos a evolução das ciências e tecnologias deve parte 
de sua história aos resultados obtidos por Fourier. 

Muitos são os exemplos práticos, e alguns talvez de seu cotidiano, que utili- 
zam a transformada de Fourier. 

Apenas para ilustrar a você, suponha que você selecione uma imagem ex- 
traída da internet, de um site que contenha informações de seu interesse. Ao 
gravar a imagem em seu computador, você pode decidir editá-la, melhorando a 
sua resolução e depois gravá-la novamente, por exemplo, no formato JPG (Joint 
Photographic Experts Group). Bem, esse procedimento lhe pareceu familiar? 
Alguma vez você já fez isso? Certamente sim, pois então, implicitamente, já fez 
uso da transformada de Fourier. Isso mesmo, toda vez que processamos ima- 
gens, mesmo que seja através de um aplicativo computacional, estamos na ver- 
dade executando uma transformação da imagem via séries de Fourier. É claro 
que estamos suprimindo muitos detalhes técnicos, mas o certo é que os resul- 
tados obtidos por Fourier são fundamentais para muitos eventos atuais. 

Outro exemplo, apenas para motivá-lo, está bem associado a um tipo de 
lazer bem conhecido, o uso das máquinas de karaokê. Nessas máquinas, as 
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músicas que conhecemos aparecem sem o vocal, apenas com o som dos instru- 
mentos, de maneira que o objetivo é reproduzir a letra da música. Bem, o isola- 
mento do vocal e sua retirada da música ocorrem pelo uso das séries de Fourier. 

Assim, situações de modulação e limpeza de sinais, como por exemplo, o 
sinal obtido pelo eletrocardiograma, processamento de imagens, como por 
exemplo, a limpeza de impressões digitais para facilitar o reconhecimento pe- 
los bancos de dados policiais, processamento de áudios, como por exemplo, o 
isolamento de instrumentos musicais para a mixagem de músicas, enfim, mui- 
tos casos utilizam a transformada desenvolvida por Fourier. 

Nessa altura você deve estar se perguntando, qual é a transformada de 
Fourier? Do que trata a série de Fourier? 

Bem, no trabalho publicado em 1822, a essência do resultado obtido por 
Fourier, e melhor concebido devido a contribuição de diversos matemáticos 
e físicos ao longo da História, é que inicialmente Fourier provou ser possível 
representar a onda de propagação através de uma combinação de funções peri- 
ódicas, que trataremos mais a frente, e o que é mais surpreendente é que essa 
combinação de funções periódicas pode, de fato, representar qualquer função 
ou comportamento variacional. Esse resultado é excelente, pois com o seu uso, 
podemos além dos exemplos já discutidos, transformar uma função matemá- 
tica mais complicada em outra função mais simples e fácil de ser manipulada. 

É claro que o trabalho envolvido na obtenção da representação correta da 
série de Fourier, que nada mais é do que essa combinação de funções periódi- 
cas, nem sempre é simples, mas tem muitas aplicações práticas. 

Não pretendemos neste livro levá-lo a um aprofundamento matemático se- 
vero e áspero, mas ao mesmo tempo nossa intenção é que você compreenda 
como obter a série de Fourier para depois poder aplicá-la em situações típicas 
de seu curso. 

Para isso, vamos caminhar de uma forma natural e detalhada, construindo 
o caminho até a série de Fourier pedaço a pedaço, de maneira que fique intuiti- 
vo compreender as ideias iniciais dos resultados de Fourier. 

Vamos lá. Preparado? 

De início, vamos entender o que é uma série infinita. 
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5.1 Conceitos sobre séries infinitas. 


Com certeza você já ouviu falar de progressão aritmética e progressão geomé- 
trica. Pois bem, uma progressão aritmética é uma sequência numérica de ter- 
mos que variam em função de uma razão aditiva, enquanto que a progressão 
geométrica é uma sequência de termos que variam em função de uma razão 
multiplicativa. Ambas podem ter seus termos somados, inclusive você apren- 
deu as fórmulas da soma de uma PA e de uma PG. Pois bem, quando decidimos 
calcular a soma de uma sequência numérica, temos uma série. 

Assim, uma série nada mais é do que a soma das parcelas de uma sequência. 

Não é muito difícil imaginar o que é uma série infinita. Isso mesmo, séries 
infinitas são somas de sequências que possuem infinitos termos. 

É claro que o fato de possuir infinitos termos gera uma pergunta inevitável: 
Se a sequência possuir infinitos termos como é possível calcular a soma dessa 
sequência? 

Bem, para avaliar a possibilidade de obter a soma, temos que recorrer a ou- 
tro conceito, o estudo da convergência da série. 

Como explicado em MATOS (2002), considerando uma sequência numéri- 
ca, é natural supor que a mesma será convergente se possuir um limite L, ou 
seja, supondo uma sequência {a n } então Mm {a n } = L . Sempre que houver um 
limitante para a sequência a mesma é dita ser convergente, caso contrário, se o 
limite for igual a infinito, a sequência é divergente. Por exemplo: 


XJ EXEMPLO 


Verifique se a sequência harmônica é ou nao convergente. 

Calculando o limite, temos: 



lim — = — = 0 


Assim, a sequência harmônica é convergente e o seu limite L = 0. 
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C7 EXEMPLO 

Verifique se a sequência \ a n } = é ou não convergente. 
Calculando o limite, temos: 


n 

lim — — = lim — — - 

n— >oo r) + 1 n— >oo n ^ 1 

n n 


lim 


J 1 + - 1 + - 
n oo 


1 + 0 


Assim, a sequência é convergente e o seu limite L = 1 . 


E T EXEMPLO 

Verifique se a sequência |a n j = 2.n é ou não convergente. 
Calculando o limite, temos: 

lim 2 .n = 2.oo = oo 

n — >oo 


Assim, a sequência é divergente. 

0 mesmo raciocínio vale para as séries infinitas, isto é, se uma série infinita possuir um 
limite diferente de infinito será uma série convergente, senão será divergente. 

Contudo, não é trivial avaliar a convergência de uma série infinita. Existem muitos resul- 
tados e teoremas relacionados com o assunto, mas que extrapolam o objetivo desse capítulo, 
que está focado em um tipo especial de série infinita, a série de Fourier. 

Contudo, é possível avaliar, não provar matematicamente, o comportamento apenas ob- 
servando o comportamento da sequência. 


XJ EXEMPLO 

Retomando a sequência harmônica, provamos que existe um limite e portanto a sequência é 
convergente. Contudo, se tivermos a seguinte série infinita: 


140 ■ CAPÍTULO 1 





00 


s = I 

n = 1 


2 

n 


A série harmônica também será convergente? 

Note que estamos somando infinítamente termos que convergem para zero, mas o valor 
da soma não para de crescer, por menores que fiquem as parcelas, o que significa que o valor 
de S converge para infinito, ou seja, a série harmônica é divergente, apesar da sequência ser 
convergente. 

Outra série muito conhecida é a destacada no exemplo 5. 


EXEMPLO 

Avaliar se a série abaixo é convergente ou não: 


Os termos da série são parcelas da conhecida expansão em série de Taylor, e em par- 
ticular a série retratada resulta no número de Euler. Assim, como a soma será diferente de 
infinito, a série é convergente. 

Não vamos nos estender sobre a convergência de séries infinitas, mas uma questão deve 
ser esclarecida neste momento. Se a série de Fourier é formada por uma sequência infinita 
de funções periódicas, ela será convergente? 

Sim, esse resultado já foi provado de diversas formas e por diversos matemáticos ao 
longo do tempo, como você poderá conferir na indicação do box conexão. 


CONEXÃO 

Para verificar a comprovação da convergência da série de Fourier, você pode acessar o 
artigo publicado pela Revista Matemática Universitária em 1 990: 
http://rmu.sbm.org.br/Conteudo/n1 1/nl 1_Artigo02.pdf 
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Do discutido nesta seção, o importante é você ter em mente que uma série 
infinita nada mais é do que a soma de infinitas parcelas resultantes de um ter- 
mo geral, que forma a série, e que essa soma pode resultar em um valor limitan- 
te, sendo a série convergente, como é o caso da série de Fourier. 

Contudo, para ficar mais claro esse comportamento convergente, vamos 
agora pensar um pouco sobre as parcelas que formam a série de Fourier, que 
dependem de funções periódicas. 

5.2 Séries periódicas e séries 
trigonométricas. 

Os resultados obtidos por Fourier estão associados a funções periódicas o que 
significa que todo o processo matemático das séries envolve uma combinação 
de funções reais que possam descrever ciclos e períodos. 

Nesse sentido, pense o seguinte. Quais funções você conhece que podem 
representar esses períodos cíclicos? Quais funções estudadas no passado são 
capazes de ilustrar um comportamento periódico? 

Bem, uma forma não muito complexa de avaliar essa questão é lembrar a re- 
presentação gráfica das principais funções estudadas até aqui e questionar se tal 
representação pode reproduzir um comportamento cíclico ou periódico. Você 
se lembra de algum caso? Se você pensou nas funções trigonométricas, acertou. 

As funções trigonométricas são uma família de funções que possuem um 
comportamento cíclico, em especial as funções seno e cosseno. Veja nas figu- 
ras 5.1 e 5.2 uma representação gráfica dessas funções. 



Figura 5.1 - Gráfico da função f(x)=sen(x). 

142 ■ CAPÍTULO 1 



Figura 5.2 - Gráfico da função f(x)=cos(x). 

É notável que as duas funções, seno e cosseno, possuem um comportamen- 
to periódico e podem, sem sombra de dúvidas, representar os resultados obti- 
dos por Fourier. Note também que tanto a função seno quanto a função cosse- 
no possuem um período, igual a 2n. O período representa o intervalo em x onde 
ocorre o ciclo das funções, que irá se repetir infinitamente, sempre no mesmo 
período. Nos dois casos, o seno e o cosseno repetem a ondulação, o ciclo, em 
intervalos de comprimento igual a 2tt. A única diferença entre as duas é um 
deslocamento de — no início do ciclo. 

Partindo das funções seno e cosseno, podemos agora pensar o seguinte: 
qualquer parametrização dessas funções gera um comportamento periódico 
constante? O que você acha, sim ou não? Infelizmente a resposta é não. Nem 
toda parametrização das funções gera um comportamento periódico constan- 
te, mas muitos possibilidades atendem ao critério. 

Vamos explanar alguns exemplos para avaliar alguns casos, o que com cer- 
teza facilitará seu entendimento da forma final da série de Fourier. 

Vejamos inicialmente o comportamento da função f(x)=sen(2x], figura 5.3, e 
da função f(x)=cos (3x), figura 5.4. 

Observando as figuras 5.3 e 5.4 notamos que ambas possuem um compor- 
tamento periódico constante, cada uma com um período específico. Assim, a 
mudança no argumento da função, de forma linear, mudou apenas o valor do 
período, mas não sua constância ao longo do domínio. 
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Figura 5.3 - Gráfico da função f(x)=sen(2x)t. 



- 1,5 


Figura 5.4 - Gráfico da função f(x)=cos(3x). 

Vejamos agora o que acontecerá se colocarmos um coeficiente multipli- 
cativo nas funções, por exemplo, funções do tipo f(x)=3.sen(2x), figura 5, e 
f(x)=4.cos(3x), figura 6. 

Observando as figuras 5 e 6 é perceptível que o comportamento é mantido, 
periodicidade constante, mas conforme mudamos o coeficiente, muda tam- 
bém o valor do período, formando ciclos mais abertos ou mais fechados, perío- 
dos mais longos ou mais curtos. 

Essa propriedade dos senos e cossenos foi fundamental para que Fourier 
pudesse concretizar os seus resultados obtendo uma série capaz de repre- 
sentar qualquer tipo de função através da combinação de períodos das séries 
trigonométricas. 
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Contudo, nem sempre modificações nos senos e cossenos serão bem 
sucedidas. 



Figura 5.5 - Gráfico da função f(x)=3.sen(2x). 



Figura 5.6 - Gráfico da função f(x)=4.cos(3x). 

Para exemplificar situações que fogem do esperado, vejamos as funções 
f(x)=sen(x 2 ) , figura 5.7, e a função f( x ) = cos [“ j , figura 8. 
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Figura 5.7 - Gráfico da função f(x)=sen(x2). 



Figura 5.8 - Gráfico da função f(x) = cos 

Desta forma, apenas argumentos ou coeficientes lineares produzem confi- 
gurações com períodos constantes. Agora, o que ocorre se fizermos uma com- 
binação de senos e cossenos? Vejamos. 

Vamos começar por uma combinação aditiva, senos somados a cossenos. 
Observe na figura 9 a função f(x)=sen(x)+cos (x) . 

É possível observar que o comportamento continuou cíclico, mas com uma 
diferença, o período é uma intermedia entre os anteriores. 
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Figura 5.9 - Gráfico da função f(x)=sen(x)+cos (x) 

Vejamos ainda um gráfico com as funções seno, cosseno e a soma das duas 
de maneira simultânea (figura 5.10). Note que três possuem diferentes perío- 
dos e ainda que a função da soma representada pela linha cheia corta o eixo 
horizontal em pontos onde o valor do seno é igual ao valor do cosseno, mas com 
sinais contrários. 

Da mesma forma, se tomarmos as funções seno e cosseno com coeficientes 
lineares e argumentos lineares, as somas também serão cíclicas, como pode ser 
confirmado nas figuras 5.11, 5.12 e 5.13. 



Figura 5.1 0 - Gráfico das funções seno, cosseno e soma de seno e cosseno. 
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Figura 5.1 1 - Gráfico da função f(x)=2.sen(3x)+3.cos(4x) 



Figura 5,1 2 - Gráfico da função f(x)=sen(2x)+cos(3x) 
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Figura 5.1 3 - Gráfico da função f{x) -sen(x) + — sen[3x) + — sen(5x) 

O o 


Você pode explorar muifas formas diferenciadas de curvas periódicas alterando os 
coeficientes e os argumentos em uma sequência de somas de senos e cossenos. Experi- 
mente utilizar o software livre Graph, o mesmo que utilizamos na confecção dos gráficos 
deste capítulo. 


Se somarmos os conceitos de séries infinitas com a visualização de que a 
soma de senos e cossenos pode gerar uma função periódica, não é difícil imagi- 
nar qual a forma geral da série de Fourier, não é? 

Relembrando, vimos que a resultado obtido por Fourier sugeria que era pos- 
sível representar uma função qualquer por uma série infinita de funções perió- 
dicas, trigonométricas. Podemos interpretar esse resultado da seguinte forma: 


f(x)=a_0 sen(0.x)+a_l sen(lx)+a_2 sen(2x)+-+b_l cos(lx)+b_2 cos(2x)+- 
Ou de forma mais resumida: 

co co 

f(x ) = ^a ; sen(/>) + ^b / cos(/>) 

/=o i=0 
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Já começamos a dar uma forma para a série de Fourier, uma soma infinita 
de funções periódicas formadas por senos e cossenos. 

O grande problema é, dada uma situação qualquer, encontrar os coeficien- 
tes a e b . 

I l 

Para entender o processo de cálculo dos coeficientes, vamos explorar inicial- 
mente dois conceitos: as séries pares e ímpares e o valor médio de uma função. 

5.3 Séries pares e ímpares. 

Antes de abordar o conceito de série par e ímpar, vejamos tais conceitos consi- 
derando funções matemáticas. Funções pares e funções ímpares estão associa- 
das a simetria do domínio e também ao comportamento simétrico, ou não da 
função. Por exemplo, você se lembra dos gráficos das funções de primeiro e de 
segundo graus, não é? 

Bem, se você considerar uma reta do tipo f(x)=ax, o resultado é uma função 
simétrica em relação a origem, mas não é simétrica em relação ao eixo hori- 
zontal, mas uma parábola do tipo f(x)=ax2, será simétrica em relação ao eixo 
horizontal, mas não em relação a origem. 

Assim, se tomarmos valores do domínio, simétricos em relação ao eixo ver- 
tical, por exemplo -2 e 2, vamos observar que para a reta, o valor da função será 
diferente nos dois pontos, mas para a parábola não, como pode ser observado 
nas figuras 5.14 e 5.15. 



Figura 5.1 4 - Gráfico da função linear - função ímpar. 
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Figura 5.1 5 - Gráfico da função quadrática - função par. 

O fato de a função acompanhar a simetria do domínio a caracteriza como 
uma função par, ou seja: 

Definição de Função Par: Seja uma função f(x) cujo domínio seja um con- 
junto simétrico. Assim, f(x) é uma função par se, e somente se, f(- x) = f(x), para 
todo x pertencente ao domínio. 

Apesar da função linear não ser simétrica em relação ao eixo horizontal, o é 
em relação a origem o que significa que o valor da função para -x e x são iguais, 
mas com sinais contrários, o que a caracteriza como uma função ímpar. Assim: 

Definição de Função ímpar: Seja uma função f(x) cujo domínio seja um 
conjunto simétrico. Assim, f(x) é uma função ímpar se, e somente se, f(- x) = - f 
(x), para todo x pertencente ao domínio. 

Voltando às figuras 5.1 e 5.2, você seria capaz de definir qual a classificação 
das funções seno e cosseno? 

Note que a função seno é uma função ímpar enquanto que a função cosseno 
é uma função par. Agora, o que você pode dizer dessas funções se elevarmos 
ambas ao quadrado? 

Bem, para responder a essa questão podemos recorrer a uma propriedade 
das funções pares e ímpares que diz: 

Propriedade 1: se f(x) for par, então f(x).f(x) será par. 

Propriedade 2: se f(x) for ímpar, então f(x).f(x) será par. 

Assim, tanto a função f{x) = sen 2 (x) quanto a função f(x) = cos 2 (x) serão 
pares. 
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Você deve estar se perguntando, qual a relevância em discutir funções pares 
e ímpares para o estudo das séries de Fourier? 

Bem, além dos aspectos relacionados com a simetria da função, temos ou- 
tros dois resultados, que apesar de serem bem intuitivos, são de grande ajuda 
para se trabalhar com a série de Fourier. Vejamos. 

Resultado 1. Considere uma função f(x) par variando no intervalo de -x a x, 
conforme indicado na figura 16. Considere ainda a área sob f(x) no intervalo de 
-x a 0 e a área sob f(x) no intervalo 0 a x, destacadas na figura 16. Note que, pela 
simetria, as áreas são iguais e com isso podemos afirmar que se uma função f(x) 
for par, vale a seguinte relação: 

0 x x x 

J f(x)dx = jf(x)dx => J f(x)dx = 2 -$f(x)dx 

-x 0 -x 0 



Figura 5.1 6 - Áreas simétricas em uma função par. 

Resultado 2. Considere uma função f(x) ímpar variando no intervalo de -x a 
x, conforme indicado na figura 5.17. Considere ainda a área sob f(x) no intervalo 
de -x a 0 e a área sob f(x) no intervalo 0 a x, destacadas na figura 5.17. Note que, 
pela simetria, as áreas são iguais, mas com sinal contrário e com isso podemos 
afirmar que se uma função f(x) for ímpar, vale a seguinte relação: 

0 x x 

j 'f(x)dx = -$f(x)dx^ $f(x)dx = 0 

—x 0 —x 
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Figura 5.1 7 - Áreas opostas em uma função ímpar. 

Note que tais resultados são importantes, pois a série de Fourier é formada 
por senos (funções ímpares) e cossenos (funções pares), mas que também será 
formada por funções do tipo seno ao quadrado e cosseno ao quadrado, que são 
funções pares e tais resultados abreviará a manipulação das integrais necessá- 
rias para a caracterização das séries de Fourier. 

Além das funções, e consequentemente séries pares e ímpares, temos ainda 
outro conceito a apresentar antes de começarmos nosso trabalho de dedução 
dos coeficientes da série de Fourier. Vamos relembrar o conceito de valor mé- 
dio de uma função. Supondo uma função f(x) contínua em um intervalo (0,x), 
sabemos que a integral definida de f(x) no intervalo equivale a área S sob a cur- 
va, ou seja: 

X 

S = J f(x)dx 
o 

Ainda, conforme SALAS et al (2005), sabemos que sempre é possível repro- 
duzir a área S através da área de um retângulo de base x com altura M, onde M 
é chamada de média da função f(x). Veja a representação desse resultado na 
figura 5.18. 
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Figura 5.18 - Área S sob a a curva F(x), de 0 a x e a área equivalente de um retângulo de 
altura M. 

O valor médio de uma função também pode ser representado com o uso dos sinais de 
maior e menor, da forma <f(x)>. 


E7 CONEXÃO 

Para maiores detalhes sobre o valor médio de uma função você pode consultar a seguin- 
te referência: 

http://www.icmc.usp.br/pessoas/andcarva/sma301/Aulas/Aula1 8.pdf 


Nessas condições, temos que a altura M, para que as áreas sejam equivalen- 
tes, deve ser igual a: 


M = 



x 


Agora, relembrado esse conceito, o que ocorre com o valor médio M se f(x) 
for uma função periódica como o seno e o cosseno? 

Observando as figuras 1 e 2, de 0 a 2n, notamos que em ambas a área S será 
uma soma de áreas contendo partes positivas e partes negativas, e devido à si- 
metria das curvas, a área positiva será igual a área negativa, o que faz com que 
M seja igual a zero. 

O mesmo comportamento ocorre se tomarmos funções que combinam o pro- 
duto de senos e cossenos, como podemos observar nas figuras 5.19, 5.20 e 5.21. 
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Figura 5.1 9 - Gráfico da função f(x)=sen(2x),sen(3x). 



Figura 5.20 - Gráfico da função f(x)=cos(2x).cos(3x). 



Figura 5.21 - Gráfico da função f(x)=sen(2x),cos(3x). 
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E sobre funções do tipo f(x)=sen2(x) e f(x)=cos2(x)? O que ocorrerá com a 
área entre 0 e 2n? Observe as figuras 5.22 e 5.23. 



Figura 5.22 - Gráfico da função f(x)=S6n 2 (x). 



Figura 5.23 - Gráfico da função f(x)=COS 2 (x). 


Note que neste caso a área S ainda será uma soma de áreas, mas todas posi- 
tivas. Assim, temos que: 

f seri 2 (x)dx 

M = í° — 

2n 

f cos' 2 (x)dx 

M = l -° K —— 

2k 
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Vamos calcular a integral definida do seno ao quadrado inicialmente. 
Temos: 


2 7i 2 n 

J sen 2 (x)dx = j 
o o 


2 ,sen 2 {x)dx 
2 


I 


[sen 2 (.*■) + sen 2 (x)] dx 
2 


2rc 

= í 


[sen 2 ( x) + sen 2 ( a-)]+[cos 2 (a-) - cos 2 (a - )] dx 


2 ? [sen 2 (a-) + cos 2 (;()]+[ sen 2 (a-) - cos 2 (a-)] dx 

J 2 

o 


^ [sen 2 (a) + cos 2 (a-)] dx ^ [sen 2 (a) -cos 2 (a)] dx 

J õ + J õ 


_ 2 r I [sen 2 (A') + cos 2 (^)] dx 2 ? [cos 2 (a) -sen 2 (a)] dx 

~ J 2 J 2 

o o 

Note que chegamos a duas relações bem conhecidas. A primeira relação en- 
tre senos e cossenos aponta que sen 2 (A) + c 
relação de interesse aponta que cos ( a ) = cos 2 
cos 2 (a-) - sen 2 (a-) = cos(2a) 

Substituindo nas integrais anteriores, temos então: 


os 2 . 

[a)=1 enq 

I a' 


- 

-sen^ — 

v2y 



Lo que quer dizer que 


2 7 o / \ , 2 ?dA 2 ?cos(2A)dA 
J sen 2 (a) g/a = J — - j 


2 

o o 


Utilizando uma substituição simples, u= 2x, verificamos então que: 


2u 


J sen 2 (x)dx = 


x sen(2x) 
2 4 


2it sen(47i) 
~2 4 


0 sen(O) 
2 4 


Verifique que o mesmo ocorrerá com a integral do cosseno ao quadrado. 
Assim, em resumo, temos os seguintes resultados: 


2n 

J sen 2 (x)dx = n 
o 
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2 71 

J cos 2 (x)dx = n 
o 

Substituindo no cálculo do valor médio M, temos: 



2n 2 7i 2 

M J%Q 5 2 (x)dx ji 1 
2n 2 k 2 

Considerando todos os resultados obtidos anteriormente, temos os seguin- 
tes valores para a média M (para cada um dos períodos de tamanho 2n): 


FUNÇÃO 

VALOR MÉDIO M 

f(x)=a.sen(ix) 

0 

f(x)=a.cos(ix) 

0 

f(x)=a.sen(ix).b.sen(jx) 

0 

f(x)=a.cos(ix).b.cos(jx) 

0 

f(x)=a.sen(ix).b.cos(jx) 

0 

f(x)=a.sen 2 (ix) 

1/2 

f(x)=a.cos 2 (ix) 

1/2 


Tabela 1 - Valores médios para as funções seno e cosseno. 


Nos tópicos discutidos até aqui, reunimos diversos resultados sobre a con- 
vergência da série de Fourier, sobre propriedades das séries pares e ímpares e 
também sobre as integrais definidas partindo do valor médio da função. Com 
todos esses resultados, estamos prontos para iniciar a obtenção dos coeficien- 
tes da série de Fourier. 
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A Série de Fourier 


Vimos anteriormente que a série de Fourier é na verdade uma expansão 
erando uma soma de senos e cossenos da forma: 

f[x) = a 0 + a^sen[x) + a^senfèx) h — + ò 1 cos(a') + ò 2 cos(2a")h — 

Ou de forma mais resumida: 


f{x ) = ^a ; sen(/>) + ^b ; cos(/>) 
/=o i=0 


Uma forma simples de procurar os valores dos coeficientes o. associados 
aos senos os valores b. associados ao cosseno é multiplicar toda a série pela fun- 
ção correspondente ao coeficiente que se deseja calcular. 

Ilustrando o processo, vamos determinar o valor do coeficiente a 2 que está 
associado a função sen(x). Assim, multiplicamos a série de Fourier por sen(x), 
obtendo o seguinte resultado: 


f(x).sen(x ) = a Q sen(x) + a^serí 1 ( x ) + a 2 sen(x)seo( 2x) — 
— b ò 1 seí?(A')cos(x) + b 2 sen(x)cos(2x) h — 


As funções resultantes lembram algum resultado obtido anteriormente? Se 
você pensou nos resultados colocados na tabela 1, acertou. Note que temos as 
mesmas funções colocadas na tabela, então se tomarmos os valores médios das 
mesmas, já sabemos quais os resultados esperados. Por uma questão de conve- 
niência, denotaremos o valor médio com a notação o. Assim: 


^f(x).sen(x)^ =< a 0 sen(x) > -b^sen 2 (x)\+ < a 2 sen(x)sen(2x) > h — 
■■■+ < b 1 sen(x)cos(x) > + < b 2 sen(x)cos(2x) > h — 
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Observando a tabela 5.1, vemos que os valores médios para os funções do 
tipo sen(x), sen(x)sen(ix) serão iguais a zero, bem como para as funções do tipo 
sen(x).cos(ix). Ainda, pela tabela 5.1 podemos observar que a média da função 
sen2(x) será igual a meio. 

Substituindo, temos: 

(f (x).sen(x)) = => a, = 2.(f(x).sen(x)) 


Como os resultados da tabela 5.1 valem para o seno e também para o cosse- 
no, podemos generalizar o resultado anterior da seguinte forma: 

a ; = 2 .(f(x).sen(i.x)^ 
b : = 2.{f(x).cos(i.x)'j 

É importante notar que os resultados anteriores não estão relacionados com 
o parâmetro aO. Neste caso, tomando a média da série, sem a multiplicação de 
qualquer função pela série, todos os termos serão mantidos como estão, senos 

e cossenos, e em todos os casos o valor médio será igual a zero. Assim a o = V M) 
Resumidamente, temos os seguintes coeficientes: 

a o=( f M) 

a ; = 2 .(f(x).sen(i.x)} 

b i = 2.(f(x).cos(i.x)^ 


Basta termos a função f(x) de interesse e determinar os valores médios para 
cada um dos coeficientes. Desde já é bom esclarecer que a série de Fourier é 
infinita e portanto não temos condições de calcular todos os coeficientes, mas 
apenas parte deles. Cabe ao usuário definir quantos coeficientes serão deter- 
minados, sempre lembrando que quanto maior o número de coeficientes mais 
próxima a série ficará da função objetivo. 
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XJ EXEMPLO 

Vejamos como trabalhar com a série de Fourier em uma caso clássico, muito frequente em 
circuitos de chaveamento digitais, a onda quadrada. Ela tem esse nome por estar associada 
ao conceito de dígitos binários, representando a ausência ou presença de um determinado 
sinal, por exemplo. Desta forma, é uma função não confínua que apresenta apenas dois re- 
sultados. A função abaixo é um exemplo de onda quadrada, que pode ser representada pela 
figura 5.24. 

f() A í 1,aji<x<(a + 1)7t 

[O, (a + 1) n < x < (a + 2) n 

Para a=0, 2, 4, 6, ... 


ji 2 jt 3it 4 ji 5ji 6jt 

Figura 5.24 - Gráfico da onda quadrada do exemplo 1 . 

Resolvendo o exemplo 1, temos que considerar que a função f(x) possui apenas dois 
valores, 0 ou 1 , em função do domínio. Faremos os cálculos para o primeiro período apenas, 
de 0 a 2 tt, pois para os demais os resultados serão os mesmos. 

Como a 0 é simplesmente a média da função f(x), e ela assume apenas dois valores, 0 e 
1 , é imediato deduzir que a 0 = 0,5. 

Observando os valores de b,-, temos que todos estão associados a cossenos, que no 
intervalo de 0 a 2 tt terão valor médio igual a zero, como visto anteriormente. Assim todos os 
coeficientes b,- são eliminados da série. 

Para calcular o valor de a lt lembramos que a 1 = 2.(f(x).sen(xf) . Para f(x)=0 temos que 
<f(x).sen(x)> será igual a zero. Para f(x)= 1 , temos que a média será igual a: 

M JÔ^n^dx -cos(x)g 2 

TC TC 71 
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Assim, temos para f(x)=0 média zero e para f(x)=1 média 2/n, então é natural que a 
média final de f(x).sen(x) seja igual a l/rr. 

Então, temos: 

a 1 = 2.(f (x).sen(x))2.— = — 

' n n 


Vamos agora calcular o valor de a v lembramos que a 2 = 2-(f(x)sen(2x Observando 
novamente a figura 3, colocada no início do capítulo, é notável observar que no intervalo de 0 
a 71 , onde f(x)= 1 , a função sen(2x) tem um comportamento ímpar, o que nos leva a concluir 
que o valor médio será igual a zero. Como para f(x)=0 o valor médio também será zero, con- 
cluímos que Ü 1 = 2.0 = 0. 

Ainda, se fizermos o gráfico das função sen(4x) ou sen(6x) esse comportamento será 
mantido. Assim, para todo a,-, com i par, o coeficiente será igual a zero. 

Resta observar o que ocorrerá com o parâmetro a 3 e os demais parâmetros com índice 
ímpar. Vejamos. 

Para calcular o valor de a 3 , lembramos que ^ = '2.(f(xysen(3x)'j p ara f( x )=0 temos 
que <f(x).sen(3x)> será igual a zero. Para f(x)= 1 , temos que a média será igual a: 

f 1.sen(3x)dx _ õ cos (3x) 0 q 

M = = = — 

71 71 371 


2 1 

Assim, temos para f(x)=0 média zero e p ^ J'<)=1 média igual a íttstt, então é natural 
que a média final de f(x).sen(3x) seja igual a i-r.in- 


Então temos: 

a 3 =2.( f (x).s en (3x))2.^ = AL 


Você consegue verificar uma lei de formação para os índices ímpares dos coeficientes? 
Não? Então façamos o cálculo do parâmetro a 3 . 

Para calcular o valor de a 5 , lembramos que 2.^/(x).sen( p ara f( x )=o temos 
que <f(x).sen(3x)> será igual a zero. Para f(x)= 1 , temos que a média será igual a: 

f\sen(5x)dx ~^ cos ( 5x )o 2 

M = 5 = — 

n n 5n 
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2 

Assim, temos para f(x)=0 média zero e para f(x)=1 média igual a — , então é natural 

07- 


que a média final de f(x).sen(3x) seja igual a 

Então temos: 

a 3 =2.(f(x).sen(3x))2.i- = |- 


E agora, será que podemos generalizar? Sim, certamente. É possível observar que para 
todos os parâmetros ímpares, vale a seguinte regra: 

2 

o; = — 

[JT 


Assim, finalmente temos a série de Fourier para o exemplo 1 , dada por: 

2 2 2 2 

f(x) = 0,5 + — sen(x) h sen(3x) h sen(5x) h sen(7x) + •■• 

71 3 71 Õ7T 7 7T 

Como já discutimos anteriormente, podemos representar, de forma finita, a série definin- 
do a quantidade de termos a serem utilizados, mas com o cuidado de lembrar que quanto 
maior o número de termos mais próxima da função original. E fácil verificar essa propriedade 
pelas figuras 5.25, 5.26 e 5.27. 



Figura 5.25 - Gráfico da série de Fourier do exemplo 1 com 3 termos. 
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Figura 5.26 - Gráfico da série de Fourier do exemplo 1 com 6 termos. 

iy 


l/WwW\/ 


\aaaw/v\/ 


Vwww\/ 



jt 


ji 5jt 3jV 7jt 

jr 9 ji 5 

' 1 í jt ' ^ 

Jt X 


2 2 2 2 2 2 


Figura 5.27 - Gráfico da série de Fourier do exemplo 1 com 1 0 termos. 

É nítido que na figura 5.27, com a série contendo 1 0 termos, o resultado está bem mais 
próximo da figura 5.24 que nos outros dois casos. É claro também que se colocarmos mais 
termos, a aproximação será cada vez melhor. 


t7 EXEMPLO 

Expandir a função f(x)=2x 2 em termos da série de Fourier, para -71 < x < n. 
Calculando o parâmetro a Q : 
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a 


o - 


f* 2x 2 dx 

J -n 

271 


2x 3 

n 

3 

-n 


2ti 


4 TI 3 
3 

2 ti 


2ti 2 

3 


Calculando o parâmetro a,: 


f 2x 2 sen(x)dx 

a, = 2 ■ — — — 

1 271 

Aplicando seus conhecimentos de cálculo integral, você pode verificar que a integral do 
numerador, sem os limites, será igual a: 

|2x 2 sen(x)dx = -2x 2 -cos(x) + 4x -sen(x) + 4-cos(x) 

Ao avaliarmos esse resultado no intervalo da função, devemos nos lembrar das proprie- 
dades das funções pares e ímpares, discutidas anteriormente. Como a função seno é ímpar, é 
natural associar que a área será igual à zero. Assim, podemos efetuar o cálculo considerando 
apenas os termos relacionados ao cosseno, da forma: 

n 

J 2x 2 sen(x)dx = [-2x 2 •cos(x) + 4-cos(x)]' 1 = 

—n 


= -2ji 2 cos(ti) + 4cos(7t) + 2ji 2 cos(ii) - 4cos(ti) = 

= 27i 2 - 4 - 2 ti 2 +4 = 0 
Então: 


= 2 - 


J 2x 2 sen(x)dx 


271 



Devido às propriedades do seno e cosseno em termos de funções pares e ímpares, o 
mesmo resultado ocorrerá com todos os termos sen(kx) da série de Fourier. Desta forma, 
temos que todos os a. serão iguais a zero. Neste caso dizemos que a série de Fourier será 
uma série de cossenos apenas. 
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Calculando agora o parâmetro b , temos: 
f 2x 2 cos(x)dx 
b, =2.— 

1 2it 

Aplicando seus conhecimentos de cálculo integral, você pode verificar que a integral do 
numerador, sem os limites, será igual a: 

j2x 2 cos(x)dx = 2x 2 -sen(x) + 4x-cos(x)-4-sen(x) 


Novamente, devido às características da função seno, temos que todos os termos de- 
pendente de sen(x) serão iguais a zero, reduzindo o cálculo da integral apenas em termos do 
cosseno. Ainda, pelo fato do cosseno ser uma função ímpar, podemos nos limitar a calcular a 
integral de zero a 71, e depois multiplicar 0 resultado final por 2 . 

Assim, temos: 

7T 71 

J 2x 2 cos(x)dx = 2- j"4x - cos(x)dx = 

-71 0 


= 2 -[4 ji -cos(7i)-4.0-cos(0)] = 2 [-471-0] = -871 


Então, temos: 

f 2x 2 cosíx)dx _o_ 

b, =2.— = 2- — = -8 

271 271 

Calculando agora o parâmetro b 2 , temos: 

f 2x 2 cos(2x)dx 
h - o jLis ' 1 

2 2 71 

Aplicando seus conhecimentos de cálculo integral, você pode verificar que a integral do 
numerador, sem os limites, será igual a: 

|2x 2 cos(2x)dx = x 2 • sen(2x) + x • cos(2x) - sen ^ x ) 

De maneira análoga ao caso anterior, a integral será reduzida a: 
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j 2x 2 cos(2x)dx = 2- jx -cos(2x)dx = 

— 71 0 

= 2 • [tc • cos(2ji) - 0 • cos(2.0)] = 2 • [71 - 0] = 2n 


Então, temos: 


b 2 = 2. 


J^2x 2 C0s(x)dx 27! 

2n 2n 


Comparando os dois parâmetros calculados, não fica evidente uma lei de formação. As- 
sim, vamos continuar calculando os parâmetros. 

Calculando agora o parâmetro b 3 , temos: 


b 3 = 2- 


J 2x 2 cos(3x)dx 


2n 


Aplicando seus conhecimentos de cálculo integral, você pode verificar que a integral do 
numerador, sem os limites, será igual a: 

2x 2 -sen(3x) 4x-cos(3x) 4-sen(3x) 


J2x 2 cos(3x)dx 


27 


De maneira análoga ao caso anterior, a integral será reduzida a: 


?_ „ . . ?4x-cos(3x) 

j 2x 2 cos(3x)dx = 2 ■ j — ^-dx = 


4n-cos(37t) 4-0-cos(3.0) 


-4it 0 

9 9 


9 9_ 


Então, temos: 


b 3 =2 


J 2x 2 cos(x)dx 


2n 


= 2 - 


8ít 

~~9~ 

2n 


8 

9 
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Calculando agora o parâmetro b 4 , temos: 


b 4 = 2, 


J 2x 2 cos(4x)dx 


2n 


Aplicando seus conhecimentos de cálculo integral, você pode verificar que a integral do 
numerador, sem os limites, será igual a: 


J2x 2 cos(4x)dx = 


•sen(4x) x-cos(4x) 


i(4x) 


16 


De maneira análoga ao caso anterior, a integral será reduzida a: 
J 2x 2 cos(4x)dx - 2- J * C ° S 1 4 *h x - 

-jt 0 4 


7i-cos(4tt) 0-cos(4.0) 


71 0 

4 4 


.4 4_ 


Então, temos: 


b 4 =2. 


J 2x 2 cos(x)dx 


2ti 


jt 



2 7i 2 


Vejamos agora se conseguimos encontrar uma lei de formaçao. Colocando os quatro 
parâmetros calculados e reescrevendo de uma forma conveniente: 


b 1 = 



_8 

I 2 


b 2 



^3 ~ 


8 

9 


8 

4 

_ 8 _ 

3 2 


_ 8 _ 

2 2 


b _ 1 _ 4 __ 8 ___ 8 _ 
2_ 2 _ 8 _ 16 _ 4 2 
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Então, podemos obter uma lei de formação dos coeficientes dos termos dos cossenos 
da série? Sim, sem dúvida. 


_8 

Para os coeficientes de índices ímpares, os coeficientes serão da forma p e para os 


8 

índices pares, 


Finalizando o exemplo 2, temos a seguinte série de Fourier: 
f W = - 3 - - ^ cos (*) + pCos(2x) - _cos(3x) + — 005(4^) 

Graficamente, considerando 9 termos da série, podemos observar o comportamento da 
expansão em relação ao gráfico original da função (figura 5.28). Note que a série, represen- 
tada pela linha pontilhada consegue fazer exatamenfe o mesmo fraçado da função f(x)= 2 x 2 . 



Figura 5.28 - Gráfico da série de Fourier e da função do exemplo2. 

Vejamos alguns exercícios para reforçar conceitos trabalhados neste capítulo. 
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El ATIVIDADES 


1 6. Verifique, utilizando limite, se as sequências abaixo são convergentes ou divergentes. 
4x-3 


a) f{x) = 

b) f{x) = 


5x 2 -4 

3x 2 +7x-6 
x + 3 


1 7. Determine a expansão em série de Fourier para a função dada abaixo: 

_ {. r , — jz s x < ü 
" i 2x, 0 < x < ít 


1 8. Determine a expansão em série de Fourier para a função dada abaixo: 

í-x,-1< x < 0 


f(x) = 


x,0 < x < 1 


REFLEXÃO 

Discutimos neste capítulo diversos conceitos sobre séries e funções periódicas que são utili- 
zadas para a composição da expansão via série de Fourier. Essas séries podem ser encontra- 
das utilizando-se a abordagem das médias das funções, mas será que temos outras maneiras 
de obter os coeficientes da série? 

Ainda, será que podemos gerar as séries sem trabalhar com as mesmas matematica- 
mente? 

Para melhor reflexão de sua parte, indicamos algumas referências para complementar 
o conteúdo deste capítulo. Em primeiro lugar, uma referência que além de se utilizar de um 
método alternativo para a obtenção dos parâmetros, associa a série de Fourier ao uso de 
osciloscópios em medidas de sinais: 

www.jdbte.com.br/wjrteleco/stfour.pdf [acessado em 22/06/201 5]. 

Ainda, a referência abaixo que apresenta outra maneira de obter os coeficientes da série: 
ZILL, D.G. e CULLEN, M.R. Equações Diferenciais. Volume 2. 3 a Ed. São Paulo: Makron 
Books, 2001 . Capítulo 1 1 , páginas de 206 a 209. 
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D CONCEITO 

Alguns termos citados neste capítulo são de especial importância e aparecerão nos capítulos 
futuros, sendo assim reforçamos seus conceitos neste glossário. 

1 4. Série: soma dos termos de uma dada sequência; 

1 5. Série infinita: soma dos termos de uma sequência infinita; 

1 6. Sequência convergente: sequência que possui um limite igual a um valor qualquer L, 
quando os termos tendem a infinito; 

1 7. Funções ou séries periódicas: ocorrem quando, a cada novo período (espaçamento 
no domínio) o comportamento da função ou da série se repete. 

18. Função par: é uma função que apresenta simetria ao longo do domínio, o que quer 
dizer que f(x) = f(-x). 

19. Função ímpar: é uma função que apresenta uma simetria em função da origem, ou 
seja, f(-x) = -f(x); 

20. Série de Fourier: é uma série infinita que contém a soma de funções periódicas tri- 
gonométricas, senos e cossenos, e serve para representar em forma expandida aproximada 
qualquer função matemática. 


ET leitura 

Para se aprofundar nos conceitos de séries de Fourier, recomendamos que você leia as 
obras listadas abaixo, relacionadas com aplicações práticas das séries de Fourier e o uso 
dessas séries em equações diferenciais: 

NETO, J.F. Aplicação da Transformada de Fourier no Processamento Digital de Imagens. 
UFPE, 1999. 

http://www.cin.uf pe.br/~ags/Sínais/Aplica 0 /oE7%E3o°/o20da%20Transformada%20 
de 0 /o20Fourier 0 /o20no%20processamento°/o20digital%20de%20imagens.pdf 
Acessado em 22/06/201 5. 

PUPIN, J.R. Introdução às Séries e Transformadas de Fourier e Aplicações no Processa- 
mento de Sinais e Imagens. DM, UFSCar, 201 1 

http://www.dm. ufscar.br/dm/attachments/article/6/TCCJosíana.pdf 
Acessado em 22/06/201 5. 

QUINTINO, LO. Séries de Fourier e Aplicações em Equações Diferenciais Parciais. Uni- 
versidade Federal de Alfenas, 201 3. 

http://www.dm. ufscar.br/dm/attachments/article/6/TCCJosíana.pdf 
Acessado em 22/06/201 5. 
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ET GABARITO 

Capítulo 1 
01 . 

a) Reescrevendo a equação diferencial, temos: 

— = sen2x -3 
dx 

dy = [sen 2x - 3 )dx 
Integrando ambos os lados 


Jdy = ^[sen2x -3)dx 
y + C, =-^-cos2a'-3a' + C 2 


y 

y 


= -^cos2a' - 3* + (C 2 

= — - cos 2 x - 3 a ' + C 
2 


c,) 


Que é a solução geral da equação diferencial. Para testar, basta derivar a solução geral, 
da forma: 


d\ cos 2 a , - 3 a ' + C I 

dy l 2 J 

dx dx 

* < 'B cosgJr ] , ^(-3A , d(Ç) 

dx dx dx dx 

^- = -~2(-sen2x) + (-3) + 0 

— = sen2x -3 
dx 
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Utilizando a condição y(0) = 3 na solução geral, temos: 


y(0) = -^cos(2.0)-3.0 + C = 3 


— cos0-0 + C = 3 
2 


— --1 + C = 3 
2 

1 1 7 

— - + C = 3=>C = 3 + — =>C = — 

2 2 2 


Assim, a solução particular dada a condição é 

1 7 

V = — cos2^-3a" + — 

/ 2 2 


b) Reescrevendo a equaçao diferencial, temos: 

d y x 
dx 

dy = xe x dx 

Integrando ambos os lados 
jdy = ^xe x dx 

Resolvendo a integral do lado direito por partes, temos: 

u = x => du = dx 

d v = e x dx => v = e x 

^udv =uv- Ji /du 

J xe x dx = xe x - j e x dx 

jxe x dx = xe x - (e x + C ^ 

Portanto, temos: 

Jdy = xe x - e x - C, 
y + C 2 = xe x - e x - C, 
y = e x (x-1)-C 1 -C 2 
y = e x (x-1) + C 


CAPÍTULO 1 


173 


Que é a solução geral da equação diferencial. Para testar, basta derivar a solução geral, 
da forma: 


dy _ c/(e JÍ (x-1) + C) 

dx dx 

dy _ d [xe x ) d ( e *) , d (C) 

dx dx dx dx 

— = e x + xe x - e x + 0 
dx 



dx 

Utilizando a condição y(0) = 1 na solução geral, temos: 
y(0) = e°(0 - 1) + C = 1 
1 • (0 - 1) + C = 1 
-1 + C = 1=>C = 2 

Assim, a solução particular dada a condição é 
y = e x (x- 1) + 2 

c) Reescrevendo a equação diferencial, temos: 

1 dy sen x 
cosx dx 4 
dy sen x 
cos xdx 4 
4dy = sen x ■ cos x dx 

Integrando ambos os lados 
J4 dy = Jsen x ■ cos xdx 

Resolvendo a integral do lado direito por substituição, temos: 

u = sen x => du = cos xdx 

r , u 2 sen 2 x 

l/ du = h C, = r C, 

J 2 1 2 1 
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Portanto, temos: 


J4 dy = - 


- + q 


2 

+ (q-q>) 


4 y + c 2 - 

9 

. serr 

7 2 

. sen 2 x _ 

4 y = + C 

2 

1 f sen 2 x _ 

y = -\ — - — i-c 


+ q 


sen 2 x C 

y = + — 

8 4 


y- 


8 


- + C 


Que é a solução geral da equaçao diferencial. Para testar, basta derivar a solução geral, 
da forma: 


dy_ 

dx 

± 

dx 


dx 


/ sen 2 x 

+c ) 

1 8 

dx 

,( sen 2 x 

1 

( 8 

J . d(C) 

dx 

dx 

— •( 2sen x-cosa') 
8 V ’ 

sen a', cos a' 


4dy = sen x-cosx dx 

Utilizando a condição y(0) = 2 na solução geral, temos: 

y(o) = — +C = 2 

— + C = 2 
8 

0 + C = 2=>C = 2 

Assim, a solução particular dada a condição é 


y- 


8 


- + 2 
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02 . 


a) é de segunda ordem, de grau 2, pois a derivada de grau mais elevado é de ordem 2. 

b) é de terceira ordem, de grau 3, pois a derivada de grau mais elevado é de ordem 3. 

c) é de primeira ordem, de grau 1 , pois a derivada de grau mais elevado é de ordem 1 . 

03. 

a) A equação diferencial é não linear. 

b) A equação diferencial é linear. 

c) A equação diferencial é linear. 

04. Reescrevendo a equação diferencial, temos: 


2 dy 

2dy 

±__ 

y 


+ y = 0 
= -ydx 



Integrando ambos os lados 


In y + C, = -^x + C 2 

In y=-^x + (C 2 -C 1 ) 

In y = -— x + C 
2 

e ln y = J 

u 

y = e'- 2 

X 

y = C- e ~2 

Representando graficamente a solução geral para C=2, C=5 e para C=10, temos o 
gráfico reproduzido na figura G4. 



176 ■ CAPÍTULO 1 



Figura G4 - Soluçoes para o exercício 4, com C=2, C=5 e C= 1 0 


Capítulo 2 

01 . 

a) y 1 (x 2 + 9) = xy 
Reescrevendo, temos: 

^(x 2 +9) = xy 

dy(x 2 +9) = xydx 
dy dx 
Y “ x 2 +9 

Como já está separada, podemos integrar: 
r dy r dx 
' y ■* x 2 + 9 

A integral do lado direito pode ser resolvida por substituição, da forma: 

u = x 2 + 9 => du = 2x.dx => dx = — 

2 
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Temos então a seguinte solução geral: 
rdy rdu 

Jy “J 

Lny + C, = -^l_nu + C 2 

Lny = LnVü + C 
Lny = Ln >/x 2 +9 + C 
y _ gl_n Vx 2 +9+C 

y = C-Vx 2 + 9 

Aplicando a condição inicial y(4) = 1 0, temos: 

1 0 = C ■ V4 2 + 9 => 1 0 = C ■ Vi 6 + 9 => 1 0 = C ■ 5 => C = 2 

Então a solução particular será: 
y = 2>/x 2 +9 

b) y' = y -3 sen x 
Reescrevendo, temos: 

— = y _3 sen x 
dx 

y 3 dy = senxdx 

Como já está separada, podemos integrar: 

Jy 3 dy = J senxdx 

y 4 

— + C, =-cosx + C 9 
4 

y 4 

— = -cosx + C 
4 

y 4 = -4-cosx + 4C 

2 

y = (-4-cosx + C)4 
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Aplicando a condição inicial y(0) = 2, temos: 

i i 

2 = (-4.cosO + C)4 =>2 = (C-4)4 =>2 4 =C-4=>C = 16 + 4 = 20 
Então a solução particular será: 


(2y 4 + x 4 )dx-xy 3 dy = 0 

Vamos verificar a homogeneidade da equação. Temos: 

g(x,y) = 2y? 4 + x 4 = 2x°y 4 + x 4 y° 
h(x, y) = -xy 3 = ?- x'y 3 

Assim, tanto a função g quanto a função h são homogêneas de grau 4. 

Fazendo a substituição y=ux e escrevendo dy como dy=u.dx+x.du, temos: 

g(x, y)dx - h(x, y)dy = 0 

(2u 4 x 4 + x 4 ) • dx - (u 3 x 4 )(u.dx + x.du) = 0 

(2u 4 x 4 + x 4 ) • dx - (u 3 x 4 )u.dx - (u 3 x 4 ) • x ■ du = 0 

(2u 4 x 4 + x 4 ) • dx - u 4 x 4 • dx - u 3 x 5 • du = 0 

(2u 4 x 4 + x 4 -u 4 x 4 ) -dx-(u 3 x 5 )-du = 0 

(u 4 x 4 + x 4 ) ■ dx - (u 3 x 5 ) • du = 0 

x 4 (u 4 + 1) ■ dx - x 5 (u 3 )vdu = 0 

Dividindo a expressão por x 4 , temos: 

(u 4 + l)-dx-x(u 3 )-du = 0 

Dividindo a expressão por u 2 , temos: 


y = (-4.cosx + 20)4 


02 . 


a) xy 3 | = 2y 4 + x 4 


Reescrevendo a equação diferencial, temos: 



dx 
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Multiplicando a expressão por 


> 4 + 1 ) 


-dx- 

x 


u 4 +1 


■du = 0 


temos: 


Como está separável, basta utilizar a integral dos dois lados. 


í-dx-j 

f f \ 

J x J 

lu 4 +lj 


Resolvendo a integral para du com substituição, temos: 

Chamando de u 4 + 1 de z, du será igual a — g e assim: 

r 1 r 1 
f-dx-f— dz = 0 
J x J 4z 

Lnx + C. — Lnz + C„ =0 
i 4 4 

Lnx--j-Ln(u 4 +l) = C 


Como fizemos y = ux, então u = — : 

x 

Lnx-— Lní ^- + 1 | = C 


Lnx — Ln 
4 


y 4 +x 4 


= C 
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Multiplicando a expressão por 4, temos: 


4.l_nx-l_n| y4+ 4 x4 1 = C 
Lnx 4 -Ln|L±*L] = C 


Ln 


Ln 


Ln 


Ln 


Ln 


y 4 + x 4 


x 4 


x 8 

= C 

y 4 + x 4 

X 8 

= L 

y 4 + x 4 

x 8 

-L 

y 4 +x 4 

X 8 



C.(y 4 + x 4 ) 


= 0 


C.(y 4 + x 4 ) 


Temos finalmente uma solução geral para a equação diferencial dada por: 

C.(y 4 +x 4 ) _1 = 0 

dy _ y 2 - x 2 
; dx 2xy 

Reescrevendo a equação diferencial, temos: 

2xy • dy = (y 2 ~ x 2 ) • dx => (y 2 - x 2 ) ■ dx - 2xy ■ dy = 0 


Multiplicando por (-1), temos: 
(x 2 -y 2 )-dx + 2xy - dy = 0 
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Vamos verificar a homogeneidade da equaçao. Temos: 


g(x,y) = x 2 - y 2 = x 2 y° - x°y 2 
h(x,y) = 2xy = 2xV 

Assim, tanto a função g quanto a função h são homogêneas de grau 2. 

Fazendo a substituição y=ux e escrevendo dy como dy=u.dx+x.du, temos: 

g(x,y)dx + h(x,y)dy = 0 

(x 2 - u 2 x 2 ) • dx + (2ux 2 )(u ■ dx + x ■ du) = 0 

(x 2 -u 2 x 2 ) - dx + (2ux 2 )u - dx + (2ux 2 ) - x - du = 0 

(x 2 - u 2 x 2 ) • dx + 2u 2 x 2 • dx + (2ux 3 ) ■ du = 0 

(x 2 - u 2 x 2 + 2u 2 x 2 ) • dx + (2ux 3 ) • du = 0 

(x 2 + u 2 x 2 ) • dx + (2ux 3 ) • du = 0 

x 2 (1 + u 2 ) ■ dx + x 2 (2ux) • du = 0 

Dividindo a expressão por x 2 , temos: 

(l + u 2 )-dx + 2ux-du = 0 

Dividindo a expressão por 2u, temos: 



Multiplicando a expressão por — v 

x(1 + u 2 J 


temos: 



Como está separável, basta utilizar a integral dos dois lados. 



Resolvendo a integral para du com substituição, temos: 


Chamando dei + u 2 de z, du será iqual a — e assim: 

2u 
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f-dx + f -dz = 0 
J x 1 z 

Lnx + C, + Lnz + C 0 = 0 
Lnx + Ln (1 + u 2 ) = C 
Ln [x-(1 + u 2 )] = C 
x-(1 + u 2 ) = e c 


Como fizemos y = ux, então 


x 2 +y 2 


= C 



Temos fínalmente uma solução geral para a equação diferencial dada por: 
-C = 0 


x 2 + y 2 


03. 

a) y' 


2x -3y^| 
2y-3x J 


Reescrevendo a equação diferencial, temos: 


s-(irl)="< 2s '- 3x » d> =-< 2x - 3s '»* 

(2y-3x)dy + (2x-3y)dx = 0 

Vejamos se a equação diferencial é ou não exata. Temos: 

g( *, y) .2x-3,=.5!M.-3 

õy 

h(x iy ).2y-3x=.?2CV.-3 
v ' õx 

Assim a equação diferencial é exata e portanto: 

= g(x,y) =i> = 2x - 3y => õf (x,y) = (2x - 3y)dx 

Calculando a integral em ambos os lados, em relação a x, temos: 

Jôf(x,y)dx = J (2x -3y)dx 
f(x,y) = x 2 -3xy + u(y) 
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Assim, temos: 


df(x,y) ,, . õ(x 2 -3xy + u(y)) 

a - h ( x -y)^— T = 2y - 3x 

õy õy 

-3x + u (y) = 2y - 3x => u'(y ) = 2y - 3x + 3x => u'(y) = 2y 

Integrando u’(y) em função de y, temos: 
u-(y) = 2y => Ju-(y)dy = J2ydy => u(y) = y 2 + C, 

Desta forma, f(x,y) será igual a : 
f(x,y) = x 2 -3xy + y 2 + C, 

Finalmente temos: 

C = x 2 - 3xy + y 2 + C, 
x 2 - 3xy + y 2 - C = 0 

ht ydy _ ~ xdx 

; x 2 + y 2 x 2 + y 2 

Reescrevendo a equação diferencial, temos: 

y(x 2 + y 2 )dy = -x(x 2 + y 2 )dx 
y(x 2 + y 2 )dy + x(x 2 + y 2 )dx = 0 
(yx 2 + y 3 )dy + (x 3 + ?xy ? 2 )dx = 0 

Vejamos se a equação diferencial é ou não exata. Temos: 
g(x,y) = x 3 + xy 2 ^^^ = 2y 

h(x,y) = yx 2 + y 3 => = 2y 

õx 

Assim a equação diferencial é exata e portanto: 

Mpl - g( x ,y) => = X 3 + xy 2 => 3f(x,y) = (x 3 + xy 2 
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Calculando a integral em ambos os lados, em relação a x, temos: 
J<3f(x,y)dx = j (x 3 + xy 2 )dx 

f M=T + í +u(y) 


Assim, temos 
ôf(x,y) 


= h(x,y): 


5 lT + T + U(y) 


= yx 2 + y 3 


õy õy 

yx 2 + u (y) = yx 2 + y 3 => u'(y) = yx 2 + y 3 - yx 2 => u'(y) = y 3 


Integrando u’(y) em função de y, temos: 


u-(y) = y 3 Ju-(y)dy = Jy 3 dy => u(y) = y — + C, 


Desta forma, f(x,y) será igual a : 


, x 4 x 2 y 2 y 4 

f(x,y) = — + — — + — + C, 
V 4 2 4 1 

, X x 4 x z y z y 4 

f ( x -y)=T + 2 4 +Ci 

f(x,y) = l(x 2 + y 2 ) 2 + C 1 

, x (x 2 + y 2 ) 2 
f(x,y) = l / ; +C, 


Finalmente temos: 


(x 2 + y 2 ) 2 
C = ^ ^- + C, 


(x 2 +y 2 ) 


-C = 0 


04 . 


a) y 1 = e 3x + 5y, com y(0) = 5 


Reescrevendo a equação diferencial, temos: 


dy 

dx 


5y = e 3x 
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A equação já está linearizada e a função P(x) = - 5 
Então: 

p(x) = e-f~ 5dx = e~ 5 í dx = e -5x 
Multiplicando na equação diferencial: 

e -5x ^ _ 5yg-5x _ e 3xg-5x 

dx 


Note que o lado esquerdo pode ser escrito como a derivada 1 e portanto: 

d (ye- 5x ) g_ 2x dx 

dx 

Integrando dos dois lados, temos: 

y e -5x _ J e -2x c | x 

y e -5x _ _2e _2x + C 
Dividindo a expressão por e~ 5x : 

y = -2e 3x + Ce 5x 

Que é a solução geral da equação diferencial. Utilizando a condição inicial y(0) 
temos: 

5 = -2e 30 + Ce 50 =>5 = -2 + C=>C = 7 

Assim, a solução particular é: 
y = -2e 3x +7e 5x 

b) dy = x 2 (3y + 1 ) dx 

Reescrevendo a equação diferencial, temos: 

= x 2 (3y + 1) => — = 3x 2 y + x 2 

^-3x 2 y = x 2 
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A equaçao já está linearizada e a função P(x) = - 3x 2 
Então: 


/ \ f-3x 2 dx 

|i(x) = eJ 



= e 


3 


Multiplicando na equação diferencial: 

e -x 3 _ 3x 2 ye~ x3 = x 2 e“ x3 
dx 


Note que o lado esquerdo pode ser escrito como a derivada 



d(ye x3 ) 
dx 


3 


Integrando dos dois lados, temos: 
ye ~ x3 = Jx 2 e~ x3 dx 


Resolvendo o lado direito por substituição, temos: 

u = e^ x3 => du = -3x 2 e _x3 dx => dx = ^ — T 

-3x 2 e~ x 


Assim, 

fx 2 e~ x3 dx = — — fdu = -— u + C = — -e ~ x3 + C 
J 3 J 3 3 


Substituindo na resolução da equação diferencial, temos: 

ye~ x3 = — — e~ x3 +C 
3 


(-X 3 ) 

Dividindo a expressão por e 



Que é a solução geral da equação diferencial. 


) e, portanto: 
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Capítulo 3 


01 . Fazendo a assocíaçao da equação diferencial com a equaçao auxiliar teremos: 
m 2 - 4m +13 = 0 

Onde determinamos A = -36, ou seja, serão duas raízes complexas e distintas, dadas por: 
m, = 2 + 3i 
m 2 = 2 - 3i 

Em associação com a solução complexa, temos s = 2 e t = 3, resultando na solução 
geral da equação diferencial: 

y = e 2x (C 1 cos (3x) + C 2 sen(3x)) 

Submetendo a solução geral às condições iniciais teremos: 

6 = e <2(0)) (C, cos (3(0)) + C 2 sen(3(0))) 

Ou seja, 

6 = e° (C, cos (0) + C 2 sen(0)) 

Ou ainda 

6= 1(C, (1) + C 2 (0)) 

Então, 

C, = 6 

E substituindo y' = 1 0 quando x = 0 na derivada da solução geral, dada por: 
y' _ 2 e 2x (C 1 cos (3x) + C 2 sen(3x)) + e 2x (3C 2 cos(3x) - 3C, sen(3x)) 

Teremos: 

1 o = 2e (2(0,) (C, cos (3(0)) + C 2 sen (3(0))) + e' 2 <°» (3C 2 cos(3(0)) - 3C, sen(3(0))) 

Ou seja, 

1 0 = 2e° (C, cos 0 + C 2 senO) + e° (3C 2 cosO - 3C t senO) 

Resultando em: 

10 = 2(1X0, (1) + C 2 (0)) + (1)(3C 2 (1) - 3C, (0)) 
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Sendo 

1 0 = 2C, + 3C, 


'2 



E a solução particular da equação diferencial será: 



02. Se associarmos a equação diferencial 
y" - 2y' + 2y = 0 

Com sua equação auxiliar dada por: 
m 2 - 2m + 2 = 0 

Encontramos as raízes complexas da equação, uma vez que A = -4 < 0. 

A solução dada por nrq = 1 + i e m 1 = 1 - i 
Chegando à solução geral dada por: 
y = e x (C, cos (x) + C 2 sen(x)) 

03. Para iniciar a resolução devemos passar a equação diferencial de segunda para o modo 
padrão, uma vez que deveremos encontrar a solução da equação diferencial a partir de uma 
solução já conhecida. Então, dividindo a equação por x, temos: 



x x 


De forma que a equaçao toma a forma: 


y"+— y' = 0 
x 


Entendemos que para a resolução, P(x) = — e y, = Inx. Então, substituindo na expressão 
para o cálculo da outra solução, temos: 



■ju(x)dx 
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Ou seja, 


y * (x) - lnx -%^ dx 


-f-dx 

e J x 


Como a integral do expoente do numerador é dada por: 
e -0 = e~ lnx = e lnx ~ =x _1 


Então, 



Resolvendo a integral por substituição simples, fazendo u = Inx, chegamos à integral 


Dessa forma a outra solução será da forma 

y 2 (x) = - 1 

A solução geral será da forma: 
y = C, Inx - C 2 

04. A equação auxiliar associada à equação linear de segunda ordem é dada por: 
y"+y' + ly = 0 

m 2 + m + — = 0 
4 

Cujas soluções são duas raízes iguais a - 0,5. Portanto, m = - 0,5. 

E como a solução geral para a equação diferencial com duas raízes iguais assume a 
forma, temos: 

y = C, e mx + C 2 xe mx 

Então, a solução será escrita como: 



y = C t e ( ~°’ 6x) + C 2 xe ( ~ 0,6x) 
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Se substituirmos os valores y = 2, quando x = 0 na solução geral ey'= 1, quando x = 
0 na primeira derivada dessa solução, então encontraremos duas equações. Resolvendo o 
sistema linear, teremos: 

í q =2 

{-0,5c, + C 2 = 1 

Encontramos, então C t = 2 e C 2 = 2, sendo a solução particular da equação diferencial 
y + 1y’+ — y = 0 , da forma: 

y = 2e ( ’ 0i6x) + 2xe w,5x) 


Capítulo 4 

01 . Aplicando a transformada de Laplace na função temos: 
a) £(t(f)) = 3£(f 5 ) 

^(0) = 3-!t = ^ 


b) L(f(t)) = 3L(t 2 )-L(e(- 7t >) + L(8) 

£ (,(,)) = 3^__L + i = 4__L + 8 

' 7 s 2+ s + 7 s s 3 s + 7 s 


c) C(f(t)) = 5£((e' -e-'f j = 5£(e 2 ' -2e'e-' + e- 2 ') = 5£(e 2í -2 + e - 2 ') 
£(/(?)) = 5£(e 2f ) - £(2) + £(e- 2! ) = — — + - 5 


s-2 s s+2 


d) 


£(/(f )) = £(2 -2cos(3f)) = 2£(l-cos(3f)) 


£(/(?)) = 2£(l-cos(3f)) = 2 


f 9 ' 


18 


s(s 2 +9)J s(s 2 + 9) 
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e) 


£(f (t)) = £(2t 3 + 4e“ 3t - 7sen 2 2t) = 2£(t 3 ) + 4£(e~ 3t ) - 7£(sen 2 2t) 
£(f(t)) = M( t 3) + 4«e-0-7£(sa„%), 5 |l + 4^-7J^ 


56 


s + 3 s(s 2 + 16) 


02. A partir da tabela de transformada de Laplace, associamos as transformadas inversas: 

a) 


CrHX 


£“i j— l = = — t 2 

1s 3 J 2!s 3 2 


b) 


Cr 


2 1 


s s 1 - 


2 1 


4 4 1 


4 3! 5! 'l 

4 1 

s 2 3!s 4 5!s 6 J 


2 1 V ,4 ,47, t 5 

= 4t - 4t 3 + 

s s 3 J 120 


c) 

£ _1 

Cr 1 


— 1 
s 2 + 49 [ 

5 


! +49 I 


W4* 


7 1 s 2 + 7 2 


= -sen{7t) 


d) 
Cr 1 


s + 1 
s 2 +2 

s + 1 
s 2 +2 


= Cr 


s 1 

/ y 

V2 


bK)J 


V 2 t + — =senV 2 t 

V2 
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03. Resolvendo a equação diferencial a partir da aplicação da transformada de Laplace e 
das condições iniciais: 

a) y' + 6y = e 3t , para y(0) = 1 

£(/') + 6£(y) = £(e 3í ) 

sF(s)-/(0) + 6F(s) = -L 
s-o 

(s + 6)F(s)-1 = — ^ 

c , x 1 , 1 s-2 

(s-3)(s + 6) (s + 6) (s-3)(s + 6) 


Por frações parciais, temos: 



Portanto, a solução da equação diferencial é: 

fít) = — e 3í + — e -6 ' 

W 9 9 

b) y" - 6y’ + 9y = t 2 e 3t , para y(0) = 2 e y'(0) = 6 


£(y ") -6£(y') + 9£(y ) = £(t 2 e 3t ) 

s 2 F(s) - sf (0) - f (0) - 6(sF(s) - f (0)) + 9F(s) = -*L 

(s-3) 

(s 2 - 6s + 9)F(s) - 2s + 6 = 21 

(s-3) 


(s-3) 2 F(s)-2(s-3) = 


F(s) 


2 

(s-3) 3 (s-3) 2 


(s-3) 3 

2(s-3) 

(s-3) 2 


2 2 
(s-3) 5 + s-3 


F(s) = 2 


(s-3) 5 



2 _í 4! 

4i[(s-3) 4+1 


A 


/ 



Portanto, a solução da equação diferencial é: 

+4p3í 

f(t) = + 2e 3í 

v ; 12 
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Capítulo 5 


01 . 


4x_3_ 4 _ J3^ ^_J3_ 

4x-3 x 2 x 2 _i irv ,x x 2 _oo oo 2 _0-0_0_q 


ai lim — g = lim = lim 

J n->co 5x 2 — 4 n->co 5x 2 4 

^“x 2 

Então a sequência é convergente pois tem limite, 
b) lim 3* 2 +7 x- 6 . n m (x + 3).(3x-2) 


n_>cc 5 — O-A 5-0 5 

X 2 oo 2 


= lim (x + 3) = oo + 3 = oo 


n->oo X + 3 n->oo x + 3 
Então a sequência é divergente pois não tem limite. 

02 . 

Utilizando o conceito de média da função, como feito anteriormente e considerando que 
a função é dividida em duas de acordo com a região do domínio, chegamos aos seguintes 
resultados: 


a 0 - a , - a 3 - a 5 - . . . - % ~ a A ~ a 6 : 


b i =b 3 = b 5 ò 2 =ò 4 =ò 6 = ... = 0 


E a série de Fourier será igual a 

f (x) = ^ + jsen(x) - ^sen(2x) + ^sen(3x) - ^sen(4x) h 


_ 2 _ 

Til 2 


2 2 

cos(x) " ^ 2 " cos ( 3x ) - ^2 ■ c°s(5 x ) - • ■ 


7t5 2 


03. 

Utilizando o conceito de média da função, como feito anteriormente e considerando que 
a função é dividida em duas de acordo com a região do domínio, chegamos aos seguintes 
resultados: 


_ 1 

a 0 “ 2 ’ a l = a 2 = a 3 = ■ ■ ■ = 0 ; 

4 

5l = ^3 = bg = . . . = -- “2 ; b 2 = Ò 4 = Ò 6 = ... = 0 

H 


E a série de Fourier será igual a 
f 1 ( x ) = \ - ^sen(x) - A-cos(3x) 


25ji 2 


3 (5x) 


4 

49 ji 2 


cos(7x) + 
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